Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Core Analysis”

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 1,477 records · Page 82

Multiphysics Coupling Methods for Molten Salt Reactor Modeling and Simulation in VERA

New multiphysics coupling capabilities for molten salt reactor (MSR) analysis have been developed in the Virtual Environment for Reactor Applications (VERA). This development consisted of two main efforts. First, a generic species transport module was added in the CTF code, which is the thermal-hydraulics (TH) code for VERA. This module uses the velocity fields for which CTF solves during the TH calculation to transport species through the core and around the primary loop. Additionally, a gas sparging model has been added to CTF to model the movement of certain species, namely, fission products such as xenon gas, to transport between the molten salt and gas bubbles present in the salt. The second effort in this development was coupling this capability to VERA’s neutron transport code MPACT. This effort focused on coupling the detailed TH transport models in CTF to MPACT to account for feedback effects in the neutron transport calculations. Finally, the thermochemistry code Thermochimica has also been coupled to VERA. Thermochimica performs pointwise calculations for chemical potential and Gibbs free energy and determines what phases are produced by the temperature, pressure, and elemental concentrations at different locations in the primary loop. These capabilities are demonstrated using a model of the Molten Salt Reactor Experiment (MSRE). Furthermore, this reactor operated at Oak Ridge National Laboratory in the 1960s, providing sources of experimental data that were used to develop the model. Various combinations of species were modeled using VERA’s new multiphysics coupling capabilities. Species distributions and reactivity effects behaved as expected for the MSRE model, demonstrating that the coupling is behaving correctly and causing appropriate feedback. The results of these calculations show the potential for VERA to be used for a wide variety of MSR analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Protocols for x-ray transient grating pump/optical probe experiments at x-ray free electron lasers

Abstract Transient grating spectroscopy is a specialized application of the four-wave-mixing methodology and constitutes a versatile technique for investigating the dynamics of vibrational, magnetic and electronic degrees of freedom of matter in a background-free fashion. Recent developments in free-electron laser sources have enabled the extension of this technique into the extreme ultraviolet range. Ongoing efforts to expand transient grating spectroscopy into the x-ray regime promise numerous advantages: (1) substantial penetration depths that allow for probing bulk material properties, (2) element specificity via specific core-excited states, and (3) short wavelengths that allow for excitation gratings with higher momentum transfer and improved spatial resolution. In this study, we comprehensively outline the procedures for conducting x-ray transient grating pump/optical probe experiment. The process encompasses the design and alignment of the experimental setup, as well as the subsequent steps involved in data acquisition and analysis. This paper is intended as a comprehensive guide for researchers interested in implementing x-ray transient grating spectroscopy, providing valuable insights into the intricacies of the experimental workflow required for this novel technique. Furthermore, we discuss the potential for extending this methodology to an x-ray pump/x-ray probe scheme, envisioning a future direction that holds promise for enhancing the capabilities and scope of x-ray transient grating spectroscopy, opening new opportunities for studying ultrafast processes with unprecedented temporal and spatial resolutions.

Optics↗

Oscillations and random walk of the soliton core in a fuzzy dark matter halo

We report that a fuzzy dark matter halo consists of a soliton core close to the center and a Navarro-Frenk-White–like density profile in the outer region. Previous investigations found that the soliton core exhibits temporal oscillations and random walk excursions around the halo center. Analyzing a set of numerical simulations, we show that both phenomena can be understood as the results of wave interference - a suitable superposition of the ground (solitonic) state and excited states in a fixed potential suffices to account for the main features of these phenomena. Such an eigenmode analysis can shed light on the evolution of a satellite halo undergoing tidal disruption. As the outer halo is stripped away, reducing the amplitudes of the excited states, the ground state evolves adiabatically. This suggests diminished soliton oscillations and random walk excursions, an effect to consider in deducing constraints from stellar heating.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The AXEAP2 program for K β X-ray emission spectra analysis using artificial intelligence

The processing and analysis of synchrotron data can be a complex task, requiring specialized expertise and knowledge. Our previous work addressed the challenge of X-ray emission spectrum (XES) data processing by developing a standalone application using unsupervised machine learning. However, the task of analyzing the processed spectra remains another challenge. Although the non-resonant K β XES of 3 d transition metals are known to provide electronic structure information such as oxidation and spin state, finding appropriate parameters to match experimental data is a time-consuming and labor-intensive process. Here, a new XES data analysis method based on the genetic algorithm is demonstrated, applying it to Mn, Co and Ni oxides. This approach is also implemented as a standalone application, Argonne X-ray Emission Analysis 2 ( AXEAP2 ), which finds a set of parameters that result in a high-quality fit of the experimental spectrum with minimal intervention. AXEAP2 is able to find a set of parameters that reproduce the experimental spectrum, and provide insights into the 3 d electron spin state, 3 d –3 p electron exchange force and K β emission core-hole lifetime.

36 MATERIALS SCIENCE↗

Key Learnings from Hydraulic Fracturing Test Site - 2 (HFTS-2), Delaware Basin

Hydraulic fractures (HF) play a dominant role in oil and gas production from unconventional reservoirs, fracture characteristics, geometry, and extent are generally not well known and are typically inferred from various monitoring methods or subsurface models. The Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp formation of the Permian (Delaware) Basin. HFTS-2, a unique joint industry program (JIP) that brings together government, academia, and industry, focused on field data acquisition and interpretation. The project has built a comprehensive HF research dataset that includes microseismic (MS), 1,500 ft of core, diagnostic formation injection tests (DFIT), fluid characterization, advanced logs, geochemistry, downhole gauges, varying completions, and permanent fiber optic (FO) cable in vertical and horizontal wells. These data were integrated and used to characterize HF geometry/dimensions and resulting depletion profiles, gain insights for completion design, and achieve a better understanding of parent well depletion effects on child well performance. The analysis and integration of the comprehensive HFTS-2 diagnostics dataset provide insights for optimal well placement and completions and improves understanding about the effects of parent well depletion. Furthermore, this dataset can be used to blind test and develop subsurface workflows for unconventional resources. Key learnings from HFTS-2 and the new technologies tested can easily be transferred to other unconventional plays.

58 GEOSCIENCES↗

SCALE Analyses of Scenarios in the TRISO-based Heat Pipe Microreactor Fuel Cycle

This report documents the application of the SCALE code to the analysis of a TRistructural-ISOtropic (TRISO)-based heat pipe microreactor (HPMR) within the context of its nuclear fuel cycle stages. The evaluation was conducted in support of the US Nuclear Regulatory Commission’s ongoing efforts to assess modeling capabilities for advanced non–light-water reactor technologies. The generic HPMR selected as a representative microreactor concept features a compact core design that incorporates TRISO fuel compacts, passive heat removal via heat pipes, and a transportable configuration intended for deployment in remote environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Revisiting the Architecture of the KOI-89 System

While high stellar obliquities observed in exoplanetary systems may be attributed to processes that tilt the planetary orbits, it is also possible that they reflect misalignments between protoplanetary disks and stellar spins. This latter hypothesis predicts the presence of coplanar multiplanetary systems misaligned with their central stars. Here we reevaluate the evidence of such an architecture that has been claimed for the KOI-89 system. An early-type star, KOI-89 has one validated transiting planet, KOI-89.01/Kepler-462b (period 84.7 days, radius 3.0 R {sub ⊕}), and one transiting planet candidate, KOI-89.02 (period 207.6 days, radius 4.0 R {sub ⊕}), where the latter exhibits transit timing variations (TTVs). A previous modeling of the stellar gravity-darkening effect in the transit light curves inferred a high stellar obliquity of ≈70°. We perform photodynamical modeling of the Kepler transit light curves and use the resulting constraints on the orbital configuration and transit times to update the gravity-darkened transit model. As a result, we find no firm evidence for the gravity-darkening effect in the transit shapes and conclude that stellar obliquity is not constrained by the data. Given the evidence for low orbital eccentricities from the dynamical analysis, the system architecture can thus be consistent with many other multitransiting systems with flat, near-circular orbits aligned with the stellar spin. We find that the TTVs imparted on its neighbor imply that KOI-89.01 has a mass ≳20 M {sub ⊕}. This would render it one of the densest known sub-Neptunes, mostly composed of a solid core. Lower masses are possible if the TTVs are instead due to an unseen third planet.

79 ASTRONOMY AND ASTROPHYSICS↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integrative analysis of CAM photosynthesis reveals its impact on primary metabolism in Yucca

Crassulacean Acid Metabolism (CAM) is an adaptation that temporally separates carbon uptake at night from photosynthesis during the day. CAM has evolved repeatedly across vascular plants, as its emergence may depend on simple regulatory changes to deeply conserved metabolic pathways. Modern CAM research relies heavily on interpretation of transcriptomic data, though regulation occurs at multiple levels following transcription. Additionally, while most research to date has focused on a handful of genes and metabolites in the core CAM pathway, the co-option of conserved regulatory and functional genes is bound to have wide-ranging effects on other aspects of primary metabolism. In this study, we integrate transcriptomic, proteomic, and metabolomic data to compare primary metabolism between the CAM species Yucca aloifolia and closely related C 3 species, Y. filamentosa. We observe minimal correlation between protein abundance and mRNA expression, suggesting significant post-transcriptional regulation in CAM species. We also find evidence of shifts in gene expression and metabolite accumulation outside of the central CAM pathway, suggesting that the shift to CAM has cascading effects across primary metabolism, especially nitrogen metabolism. Our findings provide insights into the metabolic shifts associated with CAM evolution, and highlight the complexity of its regulation at multiple biological levels.

59 BASIC BIOLOGICAL SCIENCES↗

A missing layer in COVID-19 studies: Transmission of enveloped viruses in mucus-rich droplets

Here we evaluate the influence of mucus layers on the evaporation time and transport of enveloped viruses, including SARS-CoV-2. Enveloped viruses must remain moist to be fully infective. Yet, the Wells model based on water droplets divides respiratory droplets into either quickly evaporated aerosolized particles termed droplet nuclei (<10 s) or liquid droplets that fall to the nearest surface, leaving no physical mechanism for airborne transmission of fully infective enveloped viruses over large distances (greater than a few meters). Yet, the role of mucus layers on evaporation times has not been considered even though the formation of mucus shells around liquid cores of respiratory droplets has been shown experimentally. Here we show that mucus shells increase the drying time by orders of magnitude so that enveloped virions may remain well hydrated and, thus, fully infective at substantial distances. Further, this provides a mechanism by which infective enveloped virus particles can transmit as aerosols within buildings and between buildings over extended distances. This analysis is important because public health agencies typically follow the Wells model to establish health policies including social/physical distancing guidelines.

60 APPLIED LIFE SCIENCES↗

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗

Magnetars and axion-like particles: probes with the hard X-ray spectrum

Quiescent hard X-ray and soft gamma-ray emission from neutron stars constitute a promising frontier to explore axion-like-particles (ALPs). ALP production in the core peaks at energies of a few keV to a few hundreds of keV; subsequently, the ALPs escape and convert to photons in the magnetosphere. The emissivity goes as ~T 6 while the conversion probability is enhanced for large magnetic fields, making magnetars, with their high core temperatures and strong magnetic fields, ideal targets for probing ALPs. We compute the energy spectrum of photons resulting from conversion of ALPs in the magnetosphere and then compare it against hard X-ray data from NuSTAR, INTEGRAL, and XMM-Newton, for a set of eight magnetars for which such data exists. Upper limits are placed on the product of the ALP-nucleon and ALP-photon couplings. For the production in the core, we perform a calculation of the ALP emissivity in degenerate nuclear matter modeled by a relativistic mean field theory. The reduction of the emissivity due to improvements to the one-pion exchange approximation is incorporated, as is the suppression of the emissivity due to proton superfluidity in the neutron star core. A range of core temperatures is considered, corresponding to different models of the steady heat transfer from the core to the stellar surface. For the subsequent conversion, we solve the coupled differential equations mixing ALPs and photons in the magnetosphere. The conversion occurs due to a competition between the dipolar magnetic field and the photon refractive index induced by the external magnetic field. Semi-analytic expressions are provided alongside the full numerical results. Furthermore, we also present an analysis of the uncertainty on the axion limits we derive due to the uncertainties in the magnetar masses, nuclear matter equation of state, and the proton superfluid critical temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurements of radioxenon activities during periods of gaseous release from an advanced gas-cooled reactor

Activities of radioxenon isotopes are reported from an Advanced Gas-cooled Reactor (AGR) during periods of gaseous release. Xenon isotopes are relevant to the monitoring regime of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These releases may affect detections from the International Monitoring System (IMS) by influencing the radioxenon background. Time series activity plots have been produced from both in-core monitoring and direct measurement at the point of release using a stack monitor system. Ratio plots using both data sets have been produced with results compared with a commonly chosen nuclear explosion “discrimination line.” The reported results improve the understanding of radionuclide emissions from civil nuclear power plants. This work is being performed as part of the Xenon Environmental Nuclide Analysis at Hartlepool (XENAH) collaboration between the Atomic Weapons Establishment (AWE, UK), EDF Energy (UK), Pacific Northwest National Laboratory (PNNL, USA), and the Swedish Defence Research Agency (FOI).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SCALE Analysis of a Fluoride Salt-Cooled High-Temperature Reactor in Support of Severe Accident Analysis

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a fluoride salt-cooled high-temperature reactor (FHR). Based on the preconceptual design for a small modular 236 MWth FHR developed by the University of California, Berkeley (PB-FHR-Mk1), a SCALE model of the PB-FHR-Mk1 reactor core was developed. The reactor was modeled at equilibrium state with different fuel compositions in different regions of the reactor. An iterative approach was used to interpolate and mix the burnup-dependent fuel compositions obtained through the depletion calculation of a core slice model. After demonstrating the applicability of SCALE’s multigroup (MG) approach for the simulation of the PB-FHR-Mk1, the resulting equilibrium core was studied in terms of the power profile, the flux profile, temperature reactivity coefficients, and the xenon reactivity. Furthermore, the tritium production rate in the salt coolant was determined, and the dependence of the one-group cross sections on the burnup and location in the reactor core was studied. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the equilibrium core, a zone-wise power profile, temperature feedback coefficients, the tritium production rate, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers (Final Technical Report)

This report summarizes the carbon capture research and development conducted by The State University of New York at Buffalo (UB) and GTI Energy (GTI) for award “DE-FE0031969: Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers” sponsored by the U.S. Department of Energy (DOE). The objective of this project is to develop an innovative sorbent structure of trapped small amines in HNC embedded in PEF for DAC. This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An innovative sorbent structure of trapped small amines in hierarchical nanoporous capsules (HNC) embedded in porous electrospun fibers (PEF) was developed for direct air capture (DAC). This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An interfacial polymerization process was developed, which utilized loaded amines inside mesoporous silica and trimesoyl chloride (TMC) dissolved in organic solvents as the precursors, to generate a polyamide (PA) coating layer on mesoporous silica and thus trap amines. Reaction conditions, including TMC concentration, organic solvents, reaction time, etc., for interfacial polymerization were optimized to effectively trap loaded amines, and cyclic heating-cooling operation was conducted to evaluate the coating quality. Larger pore volume mesoporous silica was also synthesized to increase amine loading and thus increase CO 2 capacity. The optimized sorbent material exhibited CO 2 capacity as high as 4.88 mmol/g under humid DAC conditions and negligible loss (<1%) during 10 cyclic heating-cooling operations. The optimized PA-coated sorbent also showed fast adsorption and desorption kinetics, with <20% t1/2 increase compared to uncoated sorbent. PEF fabrication conditions, including organic solvents for dissolving core and shell polymers, voltage, distance from the nozzle to the collection panel, etc. were adjusted to better incorporate HNC. After incorporating the optimized sorbent material into PEF, the structured sorbent had a CO 2 capacity of approximately 4.0 mmol/g under humid DAC conditions, with capacity loss of 0.17% per cycle and t1/2 increase less than 10%. A techno-economic analysis (TEA) for the process design for a DAC system based on our developed sorbent structure of trapped small amines in HNC embedded in PEF was conducted. The process design included process description and major equipment sizing and energy and mass balances in addition to scale-up research results and estimated capture cost. Aspen Adsorption Simulator was used to fit the experimentally measured breakthrough curves and extract equilibrium and kinetic data of the optimized sorbent. Our results indicated that for a DAC plant with CO 2 productivity of 3,000 tonne/year, the levelized cost of CO 2 capture was $\$$612/tonne, with the largest contribution of 44.33% from the fixed operation cost. Increasing CO 2 productivity, while maintaining similar fixed operation cost, is expected to significantly reduce the CO 2 capture cost. A sensitivity study was also conducted to understand the influence of total plant cost, sorbent cost, CO 2 concentration in the feed, sorbent mat lifetime, sorbent regeneration electricity, and adsorption blower pressure drop on the levelized cost of CO 2 capture, revealing a capture cost range of $\$$520-870/tonne.

36 MATERIALS SCIENCE↗

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗

A mechanistic model of a PWR-based nuclear power plant in response to external hazard-induced station blackout accidents

Natural hazard-induced nuclear accidents, such as the Fukushima Daiichi Accident that occurred in Japan in 2011, have significantly increased reactor safety studies in understanding nuclear power plant (NPP) responses to external hazard events such as earthquakes and floods. Natural hazards could cause the loss of offsite power in nuclear power plants, potentially leading to a Station Blackout (SBO) accident that significantly contributes to the overall risk of nuclear power plant accidents. Despite the fact that extensive research has been conducted on the station blackout accident for nuclear power plant, further understanding of these events is needed, particularly in the context of the dynamic nature of external hazards such as external flooding. This paper estimates the progression of station blackout events for a generic pressurized water reactor (PWR) in response to external flooding events. The original RELAP5-3D model of the Westinghouse four-loop design pressurized water reactor was adopted and modified to simulate the external flood-induced station blackout accident, including the short-term and long-term station blackout scenarios. A sensitivity analysis of long-term station blackout, examining reactor operation times and analyzing key parameters over time, was also conducted in this work. The results of the analyses, especially the critical timing parameters of key event sequences, provide useful insights about the time during the external flooding event, which is important for plant operators to make timely decisions to prevent potential core damage. This paper represents significant progress toward developing an integrated risk assessment framework for further identifying and assessing the effects of the critical sources of uncertainties of nuclear power plant under external hazard-induced events.

Liu, Tao↗

Temporal evolution of Pu and Cs sediment contamination in a seasonally stratified pond

There remains a lack of knowledge regarding ecosystem transfer, transport processes, and mechanisms, which influence the long-term mobility of Pu-239 and Cs-137 in natural environments. Monitoring the distribution and migration of trace radioisotopes as ecosystem tracers has the potential to provide insight into the underlying mechanisms of geochemical cycles. This study investigated the distribution of anthropogenic radionuclides Pu-239 and Cs-137 along with total organic carbon, iron, and trace element in contaminated sediments of Pond B at the Savannah River Site (SRS). Pond B received reactor cooling water from 1961 to 1964, and trace amounts of Pu-239 and Cs-137 during operations. Our study collected sediment cores to determine concentrations of Pu-239, Cs-137, and major and minor elements in solid phase, pore water and an electrochemical method was used on wet cores to determine dissolved elemental concentrations. More than 50 years after deposition, Pu-239 and Cs-137 in sediments are primarily located in the upper 5 cm in area where deposition of particulate-bound contaminants was prevalent and located between 5 and 10 cm in areas of high sedimentation, showing a limited migration of Pu-239 and Cs-137. A Factor analysis demonstrated different sediment facies across the pond resulting in a range of geochemical processes controlling accumulation of Pu and Cs. Highest concentrations appear to be controlled by particulate input from the influent canal, dominated by clay, silt, and sand minerals bearing Fe. Elevated Pu-239 in the sediments were observed in areas with high organic matter and higher deposition rate relative to the Pond B system near the outlet indicating strong association of Pu with OM and particulates. Therefore, organic matter cycling likely plays a role in Pu redistribution between sediment and overlying pond water, and deposition in organic rich sediments accumulating near the outlet. Though Pu appears to have been distributed throughout the pond, Cs-137 concentrations remained the highest near the influent canal.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗