Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion Couple Experiments”

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

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2020

Cross-field electron diffusion due to the coupling of drift-driven microinstabilities

In this paper, the nonlinear interaction between kinetic instabilities driven by multiple ion beams and magnetized electrons is investigated. Electron diffusion across magnetic field lines is enhanced by the coupling of plasma instabilities. Here, a two-dimensional collisionless particle-in-cell simulation is performed accounting for singly and doubly charged ions in a cross-field configuration. Consistent with prior linear kinetic theory analysis and observations from coherent Thomson scattering experiments, the present simulations identify an ion-ion two-stream instability due to multiply charged ions (flowing in the direction parallel to the applied electric field) which coexists with the electron cyclotron drift instability (propagating perpendicular to the applied electric field and parallel to the ExB drift). Small-scale fluctuations due to the coupling of these naturally driven kinetic modes are found to be a mechanism that can enhance cross-field electron transport and contribute to the broadening of the ion velocity distribution functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The Scattered Light Time-history Diagnostic suite at the National Ignition Facility

The Scattered Light Time-history Diagnostic (SLTD) is being implemented at the National Ignition Facility (NIF) to greatly expand the angular coverage of absolute scattered-light measurements for direct- and indirect-drive inertial confinement fusion (ICF) experiments. The SLTD array will ultimately consist of 15 units mounted at a variety of polar and azimuthal angles on the NIF target chamber, complementing the existing NIF backscatter suite. Each SLTD unit collects and diffuses scattered light onto a set of three optical fibers, which transport the light to filtered photodiodes to measure scattered light in different wavelength bands: stimulated Brillouin scattering (350 nm–352 nm), stimulated Raman scattering (430 nm–760 nm), and ω/2 (695 nm–745 nm). SLTD measures scattered light with a time resolution of ~1 ns and a signal-to-noise ratio of up to 500. Currently, six units are operational and recording data. Measurements of the angular dependence of scattered light will strongly constrain models of laser energy coupling in ICF experiments and allow for a more robust inference of the total laser energy coupled to implosions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Evaluation of subsurface transport processes of delayed gas signatures applicable to underground nuclear explosions

Radioactive gas signatures from underground nuclear explosions (UNEs) result from gas-migration processes occurring in the subsurface. The processes considered in this study either drive or retard upward migration of gases from the detonation cavity. The relative importance of these processes is evaluated by simulating subsurface transport in a dual-permeability medium for the multi-tracer Noble Gas Migration Experiment (NGME) originally intended to study some aspects of transport from a UNE. For this experiment, relevant driving processes include weak two-phase convection driven by the geothermal gradient, over pressuring of the detonation cavity, and barometric pumping while gas sorption, dissolution, radioactive decay, and usually diffusion represent retarding processes. From deterministic simulations we found that over-pressuring of the post-detonation chimney coupled with barometric pumping produced a synergistic effect amplifying the tracer-gas reaching the surface. Bounding simulations indicated that the sorption and dissolution of gases, tending to retard transport, were much smaller than anticipated by earlier laboratory studies. The NGME observations themselves show that differences in gas diffusivity have a larger effect on influencing upward transport than do the combined effects of tracer-gas sorption and dissolution, which is consistent with a Sobol’ sensitivity analysis. Both deterministic simulations and those considering parametric uncertainties of transport-related properties predict that the excess in concentration of SF 6 compared to 127 Xe as might be captured in small volumetric samples should be much smaller than the order-of-magnitude contrast found in the large-volume gas samples taken at the site. While extraction of large-volume subsurface gas samples is shown to be capable of distorting in situ gas compositions, the highly variable injection rate of SF 6 into the detonation cavity relative to that of 127 Xe at the start of the field experiment is the most likely explanation for the large difference in observed concentrations.

58 GEOSCIENCES↗
DOE OSTI2025

Selective Recovery of Critical Minerals from Simulated Electronic Wastes Via Reaction‐Diffusion Coupling

Abstract Atom‐ and energy‐efficient chemical separations are urgently needed to meet the surging demand for critical materials that has strained supply chains and threatened environmental damage. In this study, we used reaction‐diffusion coupling to separate iron, neodymium, and dysprosium ions from model feedstocks of permanent magnets, which are typically found in electronic wastes. Feedstock solutions were placed in contact with a hydrogel loaded with potassium hydroxide and/or dibutyl phosphate, resulting in complex precipitation patterns as the various metal ions diffused into the reaction medium. Specifically, we observed the precipitation of up to 40 mM of iron from the feedstock, followed by the enrichment of 73 % dysprosium, and the extraction of >95 % neodymium product at a further distance from the solution‐gel interface. We designed a series of experiments and simulations to determine the relevant ion diffusivities, D Nd =5.4×10 −10 and D Dy =5.1×10 −10 m 2 /s, and precipitation rates, k Nd =1.0×10 −5 and k Dy =5.0×10 −3 m 9 mol −3 s −1 , which enabled a numerical model to be established for predicting the distribution of products in the reaction medium. Our proof‐of‐concept study validates reaction‐diffusion coupling as an effective and versatile approach for critical materials separations, without relying on ligands, membranes, resins, or other specialty chemicals.

Wang, Qingpu [Physical and Computational Sciences ↗
DOE OSTI2020

In Situ Identification of Reaction Intermediates and Mechanistic Understandings of Methane Oxidation over Hematite: A Combined Experimental and Theoretical Study

Effective methane utilization for either clean power generation or value-added chemical production has been a subject of growing attention worldwide for decades, yet challenges persist mostly in relation to methane activation under mild conditions. Here, we report hematite, an earth-abundant material, to be highly effective and thermally stable to catalyze methane combustion at low temperatures (<500 °C) with a low light-off temperature of 230 °C and 100% selectivity to CO 2 . The reported performance is impressive and comparable to those of precious-metal-based catalysts, with a low apparent activation energy of 17.60 kcal·mol –1 . Our theoretical analysis shows that the excellent performance stems from a tetra-iron center with an antiferromagnetically coupled iron dimer on the hematite (110) surface, analogous to that of the methanotroph enzyme methane monooxygenase that activates methane at ambient conditions in nature. Isotopic oxygen tracer experiments support a Mars van Krevelen redox mechanism where CH 4 is activated by reaction with a hematite surface oxygen first, followed by a catalytic cycle through a molecular-dioxygen-assisted pathway. Surface studies with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and density functional theory (DFT) calculations reveal the evolution of reaction intermediates from a methoxy CH 3 –O–Fe, to a bridging bidentate formate b-HCOO–Fe, to a monodentate formate m-HCOO–Fe, before CO 2 is eventually formed via a combination of thermal hydrogen-atom transfer (HAT) and proton-coupled electron transfer (PCET) processes. Finally, the elucidation of the reaction mechanism and the intermediate evolutionary profile may allow future development of catalytic syntheses of oxygenated products from CH 4 in gas-phase heterogeneous catalysis.

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

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2023

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Static and dynamic properties of multi-ionic plasma mixtures

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "Pseudo-Ions in Jellium" (PIJ), based on an iso-electronic assumption (iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Further, applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Calibration of the sub-aperture backscatter system on OMEGA EP

The sub-aperture backscatter (SABS) diagnostic on the OMEGA EP Laser System is a diagnostic that is used to measure the backscattered and sidescattered light during laser–plasma interaction experiments that are relevant to high-energy-density physics and inertial confinement fusion. The diagnostic collects stimulated Brillouin scattering (SBS) UV light at around 351 nm and stimulated Raman scattering (SRS) in the visible-light regime in the 420–720-nm-wavelength range and provides spectrally and temporally resolved information. Five 1-in. light collectors, composed of a lens, ground glass diffuser, and coupling into a 300-μm fiber, are positioned behind the last steering mirror on one of the four beamlines to catch a portion of the beam cross section (~1.5%) of the emission that is scattered into the beamline. The SRS light is collected in two light collectors, combined, and transported via graded index fibers to a streaked spectrometer. The SABS-SRS streak spectrometer has a temporal and spectral resolution of 100 ps and 1 nm, respectively. Three other light collectors collect, combine, and transport the SBS signal to a Hamamatsu high-voltage photodiode, where an oscilloscope digitizes the data, providing a time resolution of better than 1 ns. To obtain an absolute energy calibration of SRS measurements, light signals of known energy and wavelength were injected into the light collectors one at a time. The resulting counts on the streak camera charge-coupled device for SRS are then correlated with the incident fluence of scattered light at the light collector in order to allow a quantitative assessment of streak camera sensitivity to determine the energy of the scattered light during experiments. The measurements were performed in situ from the light collectors to the detectors. Other offline measurements provided the transmission of the optics between the target chamber center and the light collectors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗
DOE OSTI2023

Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination

Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Non-galvanic mass transport in molten fluoride salt isothermal corrosion cells

Mass transport of corrosion products between Ni-201 test samples and 316 L capsules in a 500-hour FLiNaK corrosion experiment has been investigated. In one experiment, the capsule and samples were electrically connected for galvanic coupling, while in the other they were electrically isolated using boron-nitride. Post-corrosion characterization of Ni-201 from the two experiments did not reveal any remarkable differences. Substantial transport of Fe and Cr from the 316 L to the Ni-201 via the salt was observed, with notable diffusion of these elements in the near-surface regions of the Ni-201 samples, while a Cr-rich film was also observed on their surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Towards performing high-resolution inelastic X-ray scattering measurements at hard X-ray free-electron lasers coupled with energetic laser drivers

High-resolution inelastic X-ray scattering is an established technique in the synchrotron community, used to investigate collective low-frequency responses of materials. When fielded at hard X-ray free-electron lasers (XFELs) and combined with high-intensity laser drivers, it becomes a promising technique for investigating matter at high temperatures and high pressures. This technique gives access to important thermodynamic properties of matter at extreme conditions, such as temperature, material sound speed, and viscosity. The successful realization of this method requires the acquisition of many identical laser-pump/X-ray-probe shots, allowing the collection of a sufficient number of photons necessary to perform quantitative analyses. Here, a 2.5-fold improvement in the energy resolution of the instrument relative to previous works at the Matter in Extreme Conditions (MEC) endstation, Linac Coherent Light Source (LCLS), and the High Energy Density (HED) instrument, European XFEL, is presented. Some aspects of the experimental design that are essential for improving the number of photons detected in each X-ray shot, making such measurements feasible, are discussed. A careful choice of the energy resolution, the X-ray beam mode provided by the XFEL, and the position of the analysers used in such experiments can provide a more than ten-fold improvement in the photometrics. The discussion is supported by experimental data on 10 µm-thick iron and 50 nm-thick gold samples collected at the MEC endstation at the LCLS, and by complementary ray-tracing simulations coupled with thermal diffuse scattering calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗