Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fragmentation functions”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

B -Alkyl-borabicyclo[3.3.1]nonane Reagents Promote Closed-Shell Nickel-Catalyzed Alkylarylation Toward Encoded Cyclooctene Monomers

Access to new tailored monomers is essential to explore unprecedented polymer structures and modulate their properties. In previous work, we disclosed an iteroselective diarylation of 1,5-cyclooctadiene as an attractive platform for preparing diverse cyclooctene monomers; however, it was limited to C(sp 2 ) coupling partners. Established difunctionalization methods with C(sp 3 ) fragments rely on carbon-based radical formation, which is incompatible with 1,5-cyclooctadiene. This work describes a two-electron redox manifold for the nickel-catalyzed alkylarylation of 1,5-cyclooctadiene and discloses B -alkyl-borabicyclo[3.3.1]nonane (alkyl-9-BBN) as an effective transmetalating reagent for maintaining a polar reaction mechanism. The method provides 5,6-alkylarylated cyclooctenes suitable for ring-opening metathesis polymerization to obtain new materials. The properties of these polymers are benchmarked and fine-tuned by variation of coupling partners in the nickel catalysis. Density functional theory calculations revealed that destabilization of the pretransmetalation complex promotes the reactivity of alkyl-9-BBN in transmetalation compared to alkylboronic esters.

alkyls↗
DOE OSTI2025

Unimodal Imaging of Monovalent Metal-Chelator Complexes and Lipids by MALDI Imaging Mass Spectrometry

Careful regulation of monovalent metal ions (M + ) is necessary to maintain a functional cellular system. Of these ions, appropriate sodium (Na + ) and potassium (K + ) concentrations are particularly integral for electrochemical signaling, as well as the secondary transport of nutrients and waste. Dysregulation of M + homeostasis can disrupt these mechanisms, potentially influencing the metabolism of downstream biomolecules such as lipids. Thus, the relationship between M + abundances and related biomolecular distributions must be elucidated to better understand the physiology of healthy and disordered tissues. Traditional techniques for imaging biological metal distributions include SIMS, LA-ICP-MS, and XRF; however, these capabilities are limited to elemental analysis or the analysis of molecular fragments and must be paired with other modalities to visualize distributions of more complex biomolecules within the same or similar samples. Conversely, matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) is a powerful tool often used for mapping such biomolecular distributions, but current methods are unable to detect metals within tissue. This study illustrates a novel methodology that adds metal detection to the MALDI IMS repertoire through which the simultaneous detection of M + metals and lipids is achievable. Using a robotic sprayer for homogeneous application, on-tissue deposition of the chelator deferiprone (DEF) enables subsequent detection of the ionizable metal-chelator complex by MALDI without hindering lipid detection. Our work provides proof-of-concept data for the simultaneous detection of K + , Na + , and intact lipids using MALDI IMS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

In Vitro Selection of Antibodies Targeting Yersinia pestis Membrane Lipids Using Nanodisc-Based Antigen Presentation

Proteins are the most common targets for antibody discovery and vaccine development, but their sequence variability can limit the breadth of resulting antigens. Lipids represent an alternative class of antigens due to their structural conservation and roles in host–pathogen interactions. Here, we describe the development and optimization of an in vitro antibody selection workflow using lipid-containing nanodiscs as antigen presentation platforms to enable phage and yeast display selections under conditions adapted for these non-protein targets. Lipopolysaccharide (LPS) nanodiscs were first used as a model system to evaluate selection strategies, including competitive and subtractive approaches to reduce non-specific binders, yielding peptide and single-chain variable fragment (scFv) binders that were affinity matured to improve binding signals. The same approach was subsequently used to select scFv antibodies that recognize lipid nanodiscs prepared from Yersinia pestis membrane lipid extracts. These antibodies show binding to lipid nanodiscs derived from Y. pestis, with evidence of selectivity relative to control nanodiscs. Overall, this work establishes a workflow for antibody selection against lipid-containing nanodisc antigens and highlights practical considerations associated with these targets. The approach may be useful for generating affinity reagents to membrane-associated lipids, although further characterization is required to define antigen specificity and functional activity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A solid-state extrusion pathway extensible to upcycling aluminum scrap and end-of-life permanent magnets into ductile Al-6061/SmCo 5 structural magnetic composites

Permanent magnets offer strong potential for multifunctional structural components, but their inherent brittleness limits direct integration. Here, SolidStir® Extrusion (SSE) is used to embed SmCo5 powder into AA6061 in the solid state, and the resulting microstructure, interfacial chemistry, mechanical properties, and magnetic response are systematically characterized. SSE fragments the initial SmCo 5 powder into predominantly fine particles dispersed within an ultrafine Al matrix, with a volume fraction of ∼ 4–6 % in this proof-of-concept study. Interfacial reactions form a continuous Al–Co reaction layer, with local atomic-resolution indexing consistent with Al 13 Co 4 , together with Sm–Si-rich reaction products, while the matrix contains Sm–Co–Si-rich precipitates and Mg 2 Si, producing features from nanometers to micrometers. Compared with an Al extrudate, this architecture increases yield strength by ∼ 42 % and tensile strength by ∼ 31 %, while preserving uniform elongation, and it sustains work hardening in compression from 150 to 400 °C. The composite retained hard-magnetic behavior, with axial magnetization of 12.7 A·m 2 kg −1 at 6 T and intrinsic coercivity of 9.6 kOe. These findings demonstrate that SSE can produce Al-based load-bearing magnetic composites while preserving useful SmCo5 functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rewiring Aromatic Compound Consumption: Chromosomal Amplification and Evolution of a Foreign Pathway in Acinetobacter baylyi ADP1

Rational engineering strategies that seek to harness the remarkable diversity of microbial metabolism can be limited by incomplete biological knowledge. As described here, a novel approach to address this challenge involved replacing a native pathway for degrading lignin-derived aromatic compounds via ortho cleavage of protocatechuate in Acinetobacter baylyi ADP1 with a foreign meta-cleavage pathway that uses different enzymes, metabolites, and redox carriers. This alteration may improve lignin valorization and coordinate catabolism with bioproduction strategies. When a 14-kbp region of foreign DNA was inserted in the chromosome, the heterologous genes failed to confer growth on target substrates. Regional gene dosage was increased using a synthetic DNA fragment to promote recombination, and higher copy number enabled growth. During adaptive laboratory evolution, compensatory mutations arose that permit growth with one copy of the foreign genes. This complex metabolic remodeling was accomplished without assumptions about the impediments that initially prevented growth. To understand the changes that emerged, a novel transformation assay identified a combination of mutations sufficient for the new phenotype. Three unexpected changes were revealed: loss of one foreign enzyme, loss of one native enzyme, and loss of a two-component transcriptional regulatory system. This study establishes that large multicopy tandem arrays of poorly adapted pathway genes can confer new functions and improve understanding of metabolism.

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

Shining light on nanoscale ‘vine-on-stick’ eutectic structures in the Al-Ce-Ni system

Multiphase eutectics often comprise entangled solid phases with nanoscale periodicity, making it difficult to unravel their 3D connectivity using conventional 2D techniques. Among such systems, the three-phase eutectic Al-Al 11 Ce 3 -Al 3 Ni stands out for its ultrafine (∼100 nm) interphase spacing and promising creep resistance, yet its microstructure remains relatively unexplored despite its potential for high-temperature applications. Here, we use scanning hard X-ray microscopy with an unprecedented ∼10 nm pixel size to resolve its 3D morphology. Reconstructions reveal a novel “vine-on-stick” motif, wherein Al 11 Ce 3 wraps around Al 3 Ni pillars. Analysis of phase tortuosities confirms that Al 11 Ce 3 exhibits more convoluted morphologies than Al 3 Ni. No orientation relationship was observed between the two intermetallics, suggesting that growth is controlled by local solute gradients rather than epitaxy. Fragmentation of intermetallics along their longitudinal axis suggests a Rayleigh-type breakup mechanism in solid state. The “vine-on-stick” pattern may generalize to other alloy systems with variable interfacial anisotropies and low mutual solubilities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sensitive Detection of Structural Differences using a Statistical Framework for Comparative Crystallography

Chemical and conformational changes underlie the functional cycles of proteins. Comparative crystallography can reveal these changes over time, over ligands, and over chemical and physical perturbations in atomic detail. A key difficulty, however, is that the resulting observations must be placed on the same scale by correcting for experimental factors. We recently introduced a Bayesian framework for correcting (scaling) X-ray diffraction data by combining deep learning with statistical priors informed by crystallographic theory. To scale comparative crystallography data, we here combine this framework with a multivariate statistical theory of comparative crystallography. By doing so, we find strong improvements in the detection of protein dynamics, element-specific anomalous signal, and the binding of drug fragments.

Hekstra, Doeke R. [Harvard Univ., Cambridge, MA (U↗
DOE OSTI2025

Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series

Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.

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

Production of neutron-rich isotopes for 𝑍 ≥ 98 in the 238 U + 248 Cm reaction

Multinucleon transfer (MNT) reactions in actinide systems are a promising method to synthesize transuranium neutron-rich elements. Appropriate theoretical approaches are needed to understand the mechanisms behind MNT. We employ a microscopic approach to calculate neutron-rich isotope production in the reaction 238 U + 248 Cm system. Here, the stochastic mean-field (SMF) approach is used to calculate the primary cross sections in MNT reactions based on the quasifission and inverse quasifission processes, and a statistical de-excitation model with GEMINI ++ code to calculate the secondary fragment cross sections. The calculated cross sections using SMF and GEMINI ++ explain available experimental results for the 238 U + 248 Cm system at 𝐸 c.m. = 898.7 MeV energy. This shows the effectiveness and applicability of the quantal diffusion approach, based on the SMF theory, in heavy-ion collisions. Production of transuranium neutron-rich elements with a proton number up to 𝑍=101 is obtained with sizable cross sections. Theoretical results calculated for the 𝑍=102–105 region, for which there are no experimental data, show that the cross-section values would be lower than the microbarn level. SMF theory does not contain any adjustable parameters other than the standard parameters of the energy density functional used in the TDHF theory and is an important approach for the microscopic understanding of reaction mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dodecane Radiolysis Yields by Time-Resolved and Steady-State Methods

Liquid organic molecules are present as solvents, complexing ligands, and additives in both used nuclear fuel reprocessing solvent systems and in their subsequent nuclear waste streams. Under these extreme environments, these organic molecules are constantly exposed to ionizing radiation which promotes their radiolysis, forming a variety of short-lived, highly energetic, excited state and radical species.1-4 Here, we demonstrate new experimental results for the steady-state and time-resolved irradiations of dodecane (C12H26), a long chain, liquid, aliphatic hydrocarbon that is the prototypical solvent used for benchtop studies of aqueous-organic solvent extraction systems. When ionizing radiation interacts with neat dodecane, the energy transfer can result in molecular ionization, to give the dodecane radical cation (C12H26+•) and the solvated electron (eS–), and molecular electronic excitation (C12H26*), which rapidly produces transient carbon-centered radical fragments (CxHy•) and hydrogen atoms (H•).1-4 Studies on the initial yields of the ionization and excitation products were performed using time-resolved picosecond electron pulse radiolysis with the use of molecular probes. Using steady-state cobalt-60 gamma irradiations, the suite of products formed by dodecane radiolysis in aerated and deaerated solutions was determined. Then, using iodine as an alkyl radical scavenger, the loss of molecular iodine with dose was quantified, and by correlating with the molecular hydrogen yields of the system, the initial yields of the various carbon-centered radicals were also determined. Finally, the rates of reactions of the C12H26+• and eS– with ligands proposed for use in spent nuclear fuel reprocessing were studied as a function of temperature from 10 – 40 °C.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2025

Deer Vigilance and Movement Behavior Are Affected by Edge Density and Connectivity

ABSTRACT Animal behavior is an important component of individual, population, and community responses to anthropogenic habitat alteration. For example, antipredator behavior (e.g., vigilance) and animal movement behavior may both be important behavioral responses to the increased density of habitat edges and changes in patch connectivity that characterize highly modified habitats. Importantly, edge density and connectivity might interact, and this interaction is likely to mediate animal behavior: linear, edge‐rich landscape features often provide structural connectivity between patches, but the functional connectedness of patches for animal use could depend upon how edge density modifies animal vigilance and movement. Using remote cameras in large‐scale experimental landscapes that manipulate edge density (high‐ vs. low‐density edges) and patch connectivity (isolated or connected patches), we examined the effects of edge density and connectivity on the antipredator behavior and movement behavior of white‐tailed deer ( Odocoileus virginianus ). Deer vigilance was 1.38 times greater near high‐density edges compared to low‐density edges, regardless of whether patches were connected or isolated. Deer were also more likely to move parallel to connected high‐density edges than all other edge types, suggesting that connectivity promotes movement along high‐density edges. These results suggest that increases in edge density that accompany human fragmentation of existing habitats may give rise to large‐scale changes in the antipredator behavior of deer. These results also suggest that conservation strategies that simultaneously manipulate edge density and connectivity (i.e., habitat corridors) may have multiple effects on different aspects of deer behavior: linear habitat corridors were areas of high vigilance, but also areas where deer movement behavior implied increased movement along the habitat edge.

Bartel, Savannah L. [University of Wisconsin‐Madis↗
DOE OSTI2024

A DFT Comparison of C–C Reductive Coupling from Terminal Cyanido and Cyaphido Complexes of Nickel

The density functional theory study of the thermal C–C reductive coupling from terminal cyanido and hypothetical cyaphido complexes of [Ni(dmpe)] (dmpe = 1,2-bis(dimethylphosphino)ethane) revealed the key reaction intermediate in the reductive C–CP coupling being a σ-CC complex unlike an η 2 -aryl complex in the Ni C–CN system, as already observed in our previous studies. The reaction in THF is endothermic by 4.9 kcal/mol for cyanido with a 32.0 kcal/mol activation barrier and exothermic by 28.5 kcal/mol for cyaphido with an 11.3 kcal/mol activation barrier. To compare our results with the existing experimental data, we chose mesityl as the aryl group and also studied the CP reaction with [Pt(dmpe)] and [Pt(dmpm)] (dmpe = 1,2-bis(dimethylphosphino)methane) fragments. Our findings are consistent with the thermodynamically uphill photolytic C–CP bond activation in phosphaalkynes with Pt and a faster thermal back-reaction with [Pt(dmpe)] compared to that of [Pt(dmpm)]. Furthermore, based on the natural population analysis, when the polarity of the C–C bond is inverted, the sign of ΔG° is also inverted.

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

Intelligent Experiments through Real-Time AI: Fast Data Processing and Autonomous Detector Control for High-Energy Nuclear Experiments

The aim of this project is to develop software and hardware for fast real-time data processing and autonomous detector control and calibration for the sPHENIX and the future EIC experiments. Below summarizes Georgia Tech team efforts in the past year: 1. We developed a real-time clustering algorithm and FPGA-based pipeline architecture for processing fired pixel data from ALPIDE sensors in sPHENIX experiments. Our Columnar Clustering Co-Design introduces a hardware-aware, stream-friendly approach that segments pixel data by column pairs using a Column Pair Clustering (CPC) strategy, followed by Cluster Stitching to merge adjacent subclusters. Implemented in Vitis HLS, the pipeline comprises five stages—read-in, subclustering, stitching, analysis, and write-out—connected by tagged HLS streams with custom end-of-event signaling for robust synchronization. We designed a pipelined dataflow model optimized for throughput, low latency, and minimal buffering, enabling scalable clustering across events of arbitrary size. Our system maintains spatial precision via center-of-mass and shape key extraction and efficiently handles edge cases such as fragmented or nested clusters. Compared against DBSCAN in both software and hardware, our approach demonstrates competitive performance under FPGA constraints. 2. We also conducted a comprehensive algorithm-to-hardware co-design of connected component analysis tailored for sPHENIX experiments, focusing on real-time, low-latency processing using FPGAs and High-Level Synthesis (HLS). Starting from a Python-based particle tracking pipeline, the team translated the core logic—graph traversal via DFS and Union-Find—into an HLS-compatible C++ model, replacing dynamic memory and recursion with static arrays and pipelined control flow. The final design includes a fully streamed and dataflow-compatible Union-Find kernel optimized across five iterations, incorporating loop pipelining, array partitioning, AXI/FIFO interface tuning, and function flattening. Experimental results show up to 14.8× speedup over the CPU baseline, reducing per-graph latency to 1.58 μs and demonstrating strong resource efficiency with only ~7k LUTs and zero BRAM usage. The design maintains functional correctness against the Python reference using a Python-based C-simulation framework and Mean Squared Error metrics. This work validates the potential of HLS-driven FPGA designs for edge-level HEP data acquisition, laying a scalable foundation for future integration with real-time detector pipelines and multi-graph processing systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Bridging semantics, control specifications and assessment: A library for scalable demand flexibility controls

There is growing recognition that Demand Flexibility (DF) can play a major role in enhancing grid reliability, with building control applications emerging as key enablers for DF. However, the traditional approach to deploying new control applications in buildings, including those for DF, remains largely manual and tailored to individual buildings, making it difficult to scale. While research efforts have explored semantics-driven portability, DF controls specification, and assessment approaches, these initiatives are fragmented and limited in scope. This paper proposes a novel methodology, grounded in design science research, to integrate these elements and create a comprehensive DF controls library for both industry and academia. This approach is applied to develop the Demand FLEXibility controls LIBrary using Semantics (DFLEXLIBS), an extensible open-source library that provides DF controls for HVAC systems in Python. DFLEXLIBS enables portable, easy-to-deploy controls that abstract building-specific data points, facilitating assessment across diverse buildings. DFLEXLIBS features nine different control applications, and it is successfully implemented and tested across four virtual and two real buildings, bridging the gap between semantics-driven portability, DF controls specification, and rigorous performance assessment. Its benefits are measured by a reusability ratio greater than 90% and a functional overlap ratio of around 70% for the most common functions used in the library, significantly reducing time for deploying new controls.

Controls library↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗