Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance bound”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Safety-Related Instrumentation & Control Pilot Upgrade Initiation Phase Implementation Report

This research report (1) describes the process followed and products developed during the SR I&C Pilot Project Initial Scoping Phase, and (2) captures lessons learned. Exelon Generation and LWRS collaborated to develop a Digital Transformation Strategy as part of a larger Advanced Concept of Operations. The proposed SR I&C Pilot Upgrade provides a foundation stone for this Digital Transformation that will improve plant safety, reliability, and operational performance while lowering plant Total Cost of Ownership (TCO). Initial Scoping Phase activities for this Pilot Project have been performed in accordance with industry processes that have been adapted to better support digital upgrades. These processes include IP-ENG-001, Standard Design Process (SDP) [Reference 2], NISP-EN-04, Standard Digital Engineering Process (SDEP) [Reference 3], and Electric Power Research Institute (EPRI) Report 3002011816, Digital Engineering Guide (DEG). Completing Initial Scoping Phase Engineering and Operations, Licensing, and Project Management Activities was necessary to sufficiently bound the scope, schedule, and estimated cost of the Project to enable utility management to authorize moving into the Conceptual Design Phase. A significant finding of the Business Case Analysis (BCA) methodology developed and applied as part of this effort was that the growth rate of material costs for sustaining the operation of obsolete SR I&C equipment is accelerating. This directly contributed to the Project Economic Analysis created to justify continuing the Project. Project Initial Scoping Phase lessons learned have also been captured to assist the larger industry in understanding the Digital Transformation Strategy and SR I&C Pilot Project Initial Scoping Phase efforts. This is in keeping with the public/private partnership that has been established between the Department of Energy (DOE) and Exelon for this effort with engagement from the NRC. By addressing first-of-a-kind (FOAK) risks and capturing lessons learned, the SR I&C Pilot Upgrade Project addresses technical, regulatory, and business risks to enable subsequent implementers of similar upgrades.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Pellet cladding mechanical interaction as a potential failure mechanism during a control rod drop accident in a boiling water reactor

Boiling water reactors (BWRs) represent approximately one-third of the operating fleet in the United States, contributing significantly towards the global effort in reducing carbon emissions. Even though most of the operating fleet has been in operation for quite some time, continued advancements in new nuclear fuel (such as accident tolerant fuel) or operating regimes (such as power up-rates and higher burnup operation) necessitates similar advancements in modeling and simulation capabilities. Bison, a component of the Virtual Environment for Reactor Applications (VERA), is a high-fidelity fuel performance code able to explore the fuel performance of a wide variety of fuel types in one-, two-, and three-dimensions. Until recently, the code had not been used for analyses of BWRs. Modeling capabilities have been added for Gd-bearing UO{sub 2} and pure zirconium liners. New models have been added based upon the U.S. Nuclear Regulatory Commission (NRC) guide-lines for hydrogen pickup in Zircaloy-2 claddings and pellet-clad mechanical interaction (PCMI) failure during a reactivity insertion accident (RIA), known as a rod drop accident (CRDA) in BWRs. Implementing and/or improving Bison modeling capabilities extended its analytical reach to areas beyond its original intended purpose. This paper demonstrates one of these capabilities as a proof of concept. Recently developed models enable Bison to provide an alternate approach for cladding integrity determination in CRDA evaluations, which currently use bounding conservative estimates. As part of this demonstration, NRC guidance on hydrogen-pickup and PCMI failure models were utilized in this research. Even though more research is needed in establishing right inputs and process in this area, this paper demonstrates Bison's ability to determine cladding integrity in a CRDA evaluation. This first of a kind demonstration is a proof of concept in this area, which could potentially be extended to a number of other areas where a more accurate cladding integrity determination would be needed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Structures and function of a tailoring oxidase in complex with a nonribosomal peptide synthetase module

Nonribosomal peptide synthetases (NRPSs) are large modular enzymes that synthesize secondary metabolites and natural product therapeutics. Most NRPS biosynthetic pathways include an NRPS and additional proteins that introduce chemical modifications before, during or after assembly-line synthesis. The bacillamide biosynthetic pathway is a common, three-protein system, with a decarboxylase that prepares an NRPS substrate, an NRPS, and an oxidase. Here, the pathway is reconstituted in vitro. The oxidase is shown to perform dehydrogenation of the thiazoline in the peptide intermediate while it is covalently attached to the NRPS, as the penultimate step in bacillamide D synthesis. Structural analysis of the oxidase reveals a dimeric, two-lobed architecture with a remnant RiPP recognition element and a dramatic wrapping loop. The oxidase forms a stable complex with the NRPS and dimerizes it. We visualized co-complexes of the oxidase bound to the elongation module of the NRPS using X-ray crystallography and cryo-EM. The three active sites (for adenylation, condensation/cyclization, and oxidation) form an elegant arc to facilitate substrate delivery. The structures enabled a proof-of-principle bioengineering experiment in which the BmdC oxidase domain is embedded into the NRPS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Supercharging simulation-based inference for Bayesian optimal experimental design

Abstract Bayesian optimal experimental design (BOED) seeks to maximize the expected information gain (EIG) of experiments. This requires a likelihood estimate, which in many settings is intractable. Simulation-based inference (SBI) provides powerful tools for this regime. However, existing work explicitly connecting SBI and BOED is restricted to a single contrastive EIG bound. We show that the EIG admits multiple formulations which can directly leverage modern SBI density estimators, encompassing neural posterior, likelihood, and ratio estimation. Building on this perspective, we define a novel EIG estimator using neural likelihood estimation. Further, we identify optimization as a key bottleneck of gradient based EIG maximization and show that a simple multi-start parallel gradient ascent procedure can substantially improve reliability and performance. With these innovations, our SBI-based BOED methods are able to match or outperform by up to 22% existing state-of-the-art approaches across standard BOED benchmarks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Temperature-Time Modeling of Spent Fuel Cladding in Dry Storage Casks

This report uses thermal and decay heat modeling to investigate spent fuel performance and the potential for cladding to anneal in dry storage conditions. Annealing is an important feature to investigate in SNF cladding because it has a direct relationship to cladding response in storage, transportation and disposal conditions. Annealed cladding may have lower strength than unannealed cladding, however its increased ductility would provide better protection against rupture in high strain rate situations such as severe accidents. These consequences are the driver for this work however are outside of the scope of this report. The modeling focused on three representative storage systems, the TN-32B, MAGNASTOR with TSC-37 canister, and NUHOMS AHSM with a 32PTH2 canister. This covers the vertical dual-purpose, vertical ventilated and horizontal ventilated casks respectively. Decay heat modeling using high and low enrichment assemblies was used to bound the decay heat curves that might be expected in dry storage. To bound the temperature relationship, the storage casks were modeled starting at the design basis heat loads with heat decaying through time. Although the results are bounding there is not an attempt to maximize conservatism, rather the intent to form a reasonable upper limit on temperature that will be broadly applicable to the U.S. cask fleet. This will allow materials testing to focus on relevant conditions for annealing that may affect the U.S. spent fuel inventory. The results show a clear dependence on heat load pattern in time in Figure S-1 and Figure S-2. This dependence is due to the different assembly decay heat curves for different assemblies in preferentially loaded casks. It shows the need for careful decay heat modeling when examining in service fuel temperatures that are less than the cask design basis heat load. The body of the report shows percent cladding cutoffs of 300°C and 350 °C as well. These results can be used to inform testing and conclusions about cladding performance through time and the potential for cladding annealing during dry storage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Search for central exclusive production of top quark pairs in proton-proton collisions at $\sqrt{s}$ = 13 TeV with tagged protons

A search for the central exclusive production of top quark-antiquark pairs ($t\overline{t}$) is performed for the first time using proton-tagged events in proton-proton collisions at the LHC at a centre-of-mass energy of 13 TeV. The data correspond to an integrated luminosity of 29.4 fb -1 . The $t\overline{t}$ decay products are reconstructed using the central CMS detector, while forward protons are measured in the CMS-TOTEM precision proton spectrometer. An observed (expected) upper bound on the production cross section of 0.59 (1.14) pb is set at 95% confdence level, for collisions of protons with fractional momentum losses between 2 and 20%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Diaminopurine in Nonenzymatic RNA Template Copying

In the RNA World before the emergence of an RNA polymerase, nonenzymatic template copying would have been essential for the transmission of genetic information. However, the products of chemical copying with the canonical nucleotides (A, U, C, and G) are heavily biased toward the incorporation of G and C, which form a more stable base pair than A and U. We therefore asked whether replacing adenine (A) with diaminopurine (D) might lead to more efficient and less biased nonenzymatic template copying by making a stronger version of the A:U pair. As expected, primer extension substrates containing D bound to U in the template more tightly than substrates containing A. However, primer extension with D exhibited elevated reaction rates on a C template, leading to concerns about fidelity. Our crystallographic studies revealed the nature of the D:C mismatch by showing that D can form a wobble-type base pair with C. We then asked whether competition with G would decrease the mismatched primer extension. We performed nonenzymatic primer extension with all four activated nucleotides on randomized RNA templates containing all four letters and used deep sequencing to analyze the products. We found that the DUCG genetic system exhibited a more even product distribution and a lower mismatch frequency than the canonical AUCG system. Furthermore, primer extension is greatly reduced following all mismatches, including the D:C mismatch. Our study suggests that D deserves further attention for its possible role in the RNA World and as a potentially useful component of artificial nonenzymatic RNA replication systems.

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

Entanglement Capacity Estimates and Throughput Measurements of Quantum Channels

The throughput is an important performance metric of entangled qubit distribution quantum networks, and may be characterized by the number of distributed entangled qubit pairs per second (ebps). It is measured over physical quantum network connections using specialized instruments, including photonic entanglement sources and single photon detectors. Extensive theory has been developed to estimate the entangled qubit capacity of quantum channels using abstractions of physical connections. These two quantities both characterize the throughput performance but in different ways, and typically have been hard to relate to each other in concrete terms, in part due to the lack of precise measurements with matching analytical models and derivations. We describe measurements on a physical testbed with fiber connections of lengths 0-75 kilometers. We obtain the normalized analytic capacity estimates using the transmissivity approximations derived using single photon coincidence measurements, and convert them to bounds on throughput (measured in ebps) using a multiplier derived from co-located detector measurements. The results indicate consistent throughput measurements upper-bounded by their analytical capacity estimates across all connections. We show that previous capacity estimates using light measurements are below ebps measurements for some connections, due to the inclusion of non-representative decrease of light levels outside C-band with distance.

Rao, Nageswara↗
DOE OSTI2026

Capacities of Entanglement Distribution From a Central Source

Distribution of entanglement is an essential task in quantum information processing and the realization of quantum networks. In our work, we theoretically investigate the scenario where a central source prepares an N -partite entangled state and transmits each entangled subsystem to one of N receivers through noisy quantum channels. The receivers are then able to perform local operations assisted by unlimited classical communication to distill target entangled states from the noisy channel output. In this operational context, we define the EPR distribution capacity and the GHZ distribution capacity of a quantum channel as the largest rates at which Einstein-Podolsky-Rosen (EPR) states and Greenberger-Horne-Zeilinger (GHZ) states can be faithfully distributed through the channel, respectively. We establish lower and upper bounds on the EPR distribution capacity by connecting it with the task of assisted entanglement distillation. We also construct an explicit protocol consisting of a combination of a quantum communication code and a classical-post-processing-assisted entanglement generation code, which yields a simple achievable lower bound for generic channels. As applications of these results, we give an exact expression for the EPR distribution capacity over two erasure channels and bounds on the EPR distribution capacity over two generalized amplitude damping channels. We also bound the GHZ distribution capacity, which results in an exact characterization of the GHZ distribution capacity when the most noisy channel is a dephasing channel.

42 ENGINEERING↗
DOE OSTI2021

Organic Iodide Sorption from Dilute Gas Streams

Reprocessing used nuclear fuel releases volatile radionuclides, including 129 iodine (I), into the off-gas of a processing plant. Volatile radioiodine could be present in several forms, depending on the chemistry of the process used and the off-gas stream. Inorganic I 2 is expected to be the predominant I species in the dissolver off-gas (DOG), with minor organic iodides present. The bulk of the I is expected to volatilize into the DOG in parts-per-million-level (ppm) concentrations. In contrast, in the vessel off-gas (VOG), most of the volatile I is expected to be found as organic iodides, such as CH 3 I, C 4 H 9 I, and C 12 H 25 I. These species are expected to be present in parts-per-billion-level (ppb) concentrations but require abatement, even at low expected concentrations, to meet regulatory emissions limits in the United States. Historically, studies of I abatement by Ag-functionalized sorbents have focused on inorganic I in the DOG, but in the last few years, more research attention has been given to organic iodides, especially longer chain species, such as C 4 H 9 I, and C 12 H 25 I. This report has three main goals: (1) to present new data generated at Oak Ridge National Laboratory (ORNL) in FY21 on the sorption behavior of organic iodides on AgZ, (2) to summarize and synthesize organic iodide data produced by ORNL and Idaho National Laboratory (INL) over the last 4 years to answer questions on organic iodides behavior outlined in the 2018 joint test plan (Jubin et al. 2018), and (3) to propose a VOG abatement system design that can provide the capture efficiencies required to meet I emission limits. The 2021 ORNL experimental campaign tested the effects of organic iodide speciation and concentration in the off-gas, superficial velocity of the off-gas, and effects of aging on AgZ sorbent capacity. These studies found that the sorption rate of organic iodides by AgZ depends on the hydrocarbon chain length and the concentration in the off-gas. Higher molecular weight organic iodides adsorb to AgZ more slowly than I. At a concentration of 50 ppm concentration in the off-gas, CH 3 I loads 8% slower, C 4 H 9 I loads 20% slower, and C 12 H 25 I loads 40% slower than I. The lowest concentration loading rates calculated were in 5 ppm organic iodide gas streams in which AgZ gained on average 0.14 mg I/g sorbent/hour in the bench scale test system. Thus, longer sorbent beds might be needed to accommodate slower loading rates onto AgZ in lower concentration gas streams. Although sorption rate varies as a function of hydrocarbon chain length, the saturation concentration of the sorbent for these I-bearing species does not vary. Aging AgZ in a humid air stream for 9 months drops the overall sorbent capacity by ~35% for CH 3 I, ~50% for C 4 H 9 I, and ~40% for C 12 H 25 I. This results in a saturation capacity between 35 and 70 mg I/g sorbent for the aged AgZ. In conjunction with recent data produced by INL, these data are used to estimate the mass transfer zone (MTZ) and decontamination factor (DF) for sorbent beds of AgZ. Sorption tests performed with iodide gas concentrations of about 1 ppm and higher at a superficial gas velocity of 10 m/min, indicate that MTZ depths for these conditions tend to range between about 8-20 cm. Tests performed at lower concentrations between 50-90 ppb and at gas superficial velocities of 1, 10, and 20 m/min indicate that the MTZ depth increases with increasing superficial gas velocity. The 20 m/min test indicates that the MTZ for those conditions was at least 22 cm, and doubling the superficial gas velocity from 10 to 20 m/min could roughly double or triple the MTZ depth. Doubling and tripling the bounding MTZ depth of 20 cm for the body of MTZ estimates made at with 10 m/min superficial velocity would extend the MTZ for a superficial velocity of 20 m/min to 40-60 cm. This bounding limit applies to all of the organic iodides that have been tested. These results also indicate that the sorption rate-limiting step is not sensitive to the superficial gas velocity; otherwise the MTZ depth would not have increased approximately in proportion to the increase in the gas superficial velocity. This further suggests that the rate limiting step is not associated with mass transfer of the sorbate to the sorbent surface, or mass transfer of the reaction byproducts from the sorbent surface, but is associated with sorption or chemical reactions on the sorbent surface or in sorbent pores. Deep bed testing at INL has established DFs of >2,000 for I, CH 3 I, and C 4 H 9 I under a range of conditions (Soelberg et al. 2021, Bruffey et al. 2019). DF does not seem to be affected by the concentration of the organic iodide in the gas stream over the range of 1 to 50 ppm. Thus, if the MTZ is accommodated in sorbent bed design for the DOG and VOG, then regulatory DFs will be met. To meet the third objective outlined in this report, these experimental data were used to update an engineering evaluation of the VOG first completed in 2016. The updated VOG design can be found in an accompanying document (Welty et al., 2021; INL- LTD-21-64587). This report finds that the VOG will decrease in both size and complexity, relative to previous designs, and will still meet regulatory requirements for all iodine forms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

PROTEUS: Machine Learning Driven Resilience for Extreme-scale Systems

The objective of this project is to design, develop, and evaluate scalable software to enhance resilience, data checkpointing, program restart, and analysis. The proposed tasks are to 1) develop scalable machine learning techniques to learn temporal change patterns in a scalable and in-situ manner, and to minimize data movement and maximize learning locally closest to data; 2) design a concise data representation and indexing mechanism to capture the distribution of changes in data that can guarantee point-wise user-defined tolerable errors while reducing the data storage requirements by an order of magnitude or more; 3) develop data reduction techniques as library modules; 4) exploit local SSD for minimizing data movement in storage hierarchy; 5) develop anomaly detection algorithms that can predict corruptions based on learning of emerging patterns; 6) develop software libraries to be incorporated within widely used data formats and APIs; and 7) evaluate the proposed software using DOE scientific applications. The outcomes of the proposed work are to satisfy many synergistic data reduction and resilience requirements for large-scale data intensive applications executed on extreme-scale computing systems. The developed mechanism for error-bound data approximation is directly applicable to existing scientific applications. Through machine learning from historical events and change distribution, this work will enable anomaly detection for DOE computer facility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nucleon off-shell structure and the free neutron valence structure from A=3 inclusive electron scattering measurements

Understanding the differences between the distribution of quarks bound in protons and neutrons is key for constraining the mechanisms of SU(6) spin-flavor symmetry breaking in Quantum Chromodynamics (QCD). While vast amounts of proton structure measurements were done, data on the structure of the neutron is much more spars as experiments typically extract the structure of neutrons from measurements of light atomic nuclei using model-dependent corrections for nuclear effects. Recently the MARATHON collaboration performed such an extraction by measuring inclusive deep-inelastic electron-scattering on helium-3 and tritium mirror nuclei where nuclear effects are expected to be similar and thus be suppressed in the helium-3 to tritium ratio. Here we evaluate the model dependence of this extraction by examining a wide range of models including the effect of using instant-form and light-cone nuclear wave functions and several different parameterizations of nucleon modification effects, including those with and without isospin dependence. We find that, while the data cannot differentiate among the different models of nuclear structure and nucleon modification, they consistently prefer a neutron-to-proton structure function ratio of at $x_B \rightarrow 1$ of $\sim 0.4$ with a typical uncertainty ($1\sigma$) of $\sim0.05$ and $\sim0.10$ for isospin-independent and isospin-dependent modification models, respectively. While strongly favoring SU(6) symmetry breaking models based on perturbative QCD and the Schwinger-Dyson equation calculation, the MARATHON data do not completely rule out the scalar di-quark models if an isospin-dependent modification exist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2023

The influence of subhaloes on host halo properties

Within the ΛCDM cosmology, dark matter haloes are composed of both a smooth component and a population of smaller gravitationally bound subhaloes. These components are often treated as a single halo when properties, such as density profiles, are extracted from simulations. Recent work has shown that density profiles change substantially when subhalo mass is excluded. In this paper, we expand on this result by analysing three specific host halo properties – concentration (cNFW), spin (λB), and shape (c/a) – when calculated only from the smooth component of the halo. This analysis is performed on both Milky Way-mass haloes and cluster-mass haloes in high-resolution zoom-in N-body simulations. We find that when subhaloes are excluded, the median value of (1) cNFW is enhanced by $\approx 30\pm 11$ and $\approx 77\pm 8.1~{{\ \rm per\ cent}}$ for Milky Way-mass ($10^{12.1}\, \text{M}_\odot$) and cluster-mass ($10^{14.8}\, \text{M}_\odot$) haloes, respectively, (2) λB is reduced for Milky Way-mass by $\approx 11\pm 4.9~{{\ \rm per\ cent}}$ and cluster-mass haloes by $\approx 27\pm 3.5~{{\ \rm per\ cent}}$. Additionally, with the removal of subhaloes, cluster-mass haloes tend to become more spherical as the ratio of minor-to-major axis, c/a, increases by $\approx 11\pm 3.6~{{\ \rm per\ cent}}$, whereas Milky Way-mass haloes remain approximately the same shape with c/a changed by $\approx 1.0\pm 5.8~{{\ \rm per\ cent}}$. Fractional changes of each of these properties depend primarily on the amount of mass in subhaloes and, to a lesser extent, mass accretion history. Our findings demonstrate that the properties of the smooth components of dark matter haloes are biased relative to the total halo mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Transversity Distributions and Tensor Charges of the Nucleon: Extraction from Dihadron Production and Their Universal Nature

We perform the first global quantum chromodynamics (QCD) analysis of dihadron production for a comprehensive set of data in electron-positron annihilation, semi-inclusive deep-inelastic scattering, and proton-proton collisions, from which we extract simultaneously the transversity distributions of the nucleon and π + π − dihadron fragmentation functions. We incorporate in our fits known theoretical constraints on transversity, namely, its small- x asymptotic behavior and the Soffer bound. We furthermore show that lattice-QCD results for the tensor charges can be successfully included in the analysis. This resolves the previously reported incompatibility between the tensor charges extracted from dihadron production data and lattice QCD. We also find agreement with results for the transversity and tensor charges obtained from measurements on single-hadron production. Overall, our work demonstrates for the first time the universal nature of all available information for the transversity distributions and the tensor charges of the nucleon. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Explicit model predictive control through robust optimization

A strategy that calculates an explicit state feedback policy to regulate constrained uncertain discrete‐time uncertain linear systems is presented. We consider uncertain processes, affected by box‐bounded multiplicative uncertainty as well as bounded additive uncertainty with linear state and inputs constraints. The proposed method includes (i) the calculation of a terminal set constraint and (ii) the robust reformulation of state constraints in the prediction horizon. These features allow the derivation of the desired policy by solving a single multiparametric quadratic programming problem that guarantees feasible operation in the presence of uncertainty. Additionally, we employ variable and constraint elimination approaches to enhance the computational performance of the strategy. We demonstrate the steps and benefits of these developments with a numerical example and a chemical engineering case study.

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

Capturing Ultrafast Electron Driven Chemical Reactions in Molecules

The goals of this project are to develop experimental tools to investigate the dynamics of dissociative electron attachment (DEA) reactions in real time. In a DEA reaction, the attachment of a low energy electron to a molecule results in bond breaking and the fragmentation of the molecule. A new instrument has been designed and is under construction to perform time-resolved photoelectron and photo-fragment spectroscopy following electron transfer from an iodide anion to a molecule within a mass-selected molecular dimer. Iodide binds only weakly through the dipole interaction with many polar molecules, and such molecules forming a dimer with iodide undergo only small changes from the equilibrium geometry in isolation. When the electron is transferred from the iodide, the iodine atom is no longer bound, so the electron transfer step can be considered analogous to electron attachment and brings the key advantage of ultrafast timing precision. The experiments proposed here will establish the experimental tools to interrogate the dynamics of dissociative electron attachment on the fundamental timescales of the nuclear motion, by simultaneously probing the valence electronic structure and the fragmentation pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗