Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-T-196: Heat Diffusion Modeling for Digital Holographic Interferometry Dosimetry

Purpose: We have previously demonstrated that with Digital Holographic Interferometry (DHI) 2D spatial calorimetric measurements of high dose rate radiation sources can be obtained. The impact of heat transfer must be considered when undertaking any form of calorimetric measurement, as the radiation induced temperature distributions are subject to degradation due to heat diffusion. Unaccounted for, this limits the accuracy of the approach especially for long delivery times. Methods: 3D modelling of the heat diffusion in water was undertaken, and two different approaches developed to account for this effect. The mathematical framework to describe heat diffusion in 3D was applied, with the differential equations solved numerically using an implicit method. The first approach involved the comparison of the DHI measurements to an independent dose model of the source. The model was forward modeled to account for the heat diffusion during irradiation, allowing a direct comparison to validate the measured results. The second approach involved the correction of the measured data directly, by comparing the temperature distribution of two instances and subtracting the effects of heat diffusion of the first distribution from the second instance. This required the use of the Abel transform to approximate the 3D dose distribution from the 2D DHI results, thus limiting the approach to radiation applications possessing cylindrical symmetry. Results: The first approach resulted in higher accuracy and was more straightforward, but has a major limitation in that the measured results are only able to be utilized in comparison with an independent dose model. The applicability of the second approach is affected by noise in the measurement data and introduces higher uncertainties, but results in higher usability of the final data. Conclusion: Both approaches were implemented, and if used in conjunction would provide the most utility for the interpretation and use of DHI measurements.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

The Optical–Mid-infrared Extinction Law of the l = 165° Sightline in the Galactic Plane: Diversity of the Extinction Law in the Diffuse Interstellar Medium

Understanding the effects of dust extinction is important to properly interpret observations. The optical total-to-selective extinction ratio, R{sub V}=A{sub V}/E(B−V), is widely used to describe extinction variations in ultraviolet and optical bands. Since the R{sub V}=3.1 extinction curve adequately represents the average extinction law of diffuse regions in the Milky Way, it is commonly used to correct observational measurements along sightlines toward diffuse regions in the interstellar medium. However, the R{sub V} value may vary even along different diffuse interstellar medium sightlines. In this paper, we investigate the optical–mid-infrared (mid-IR) extinction law toward a very diffuse region at l=165{sup ∘} in the Galactic plane, which was selected based on a CO emission map. Adopting red clump stars as extinction tracers, we determine the optical–mid-IR extinction law for our diffuse region in two APASS bands (B,V), three XSTPS-GAC bands (g,r,i), three 2MASS bands (J,H,K{sub s}), and two WISE bands (W1,W2). Specifically, 18 red clump stars were selected from the APOGEE–RC catalog based on spectroscopic data in order to explore the diversity of the extinction law. We find that the optical extinction curves exhibit appreciable diversity. The corresponding R{sub V} ranges from 1.7 to 3.8, while the mean R{sub V} value of 2.8 is consistent with the widely adopted average value of 3.1 for Galactic diffuse clouds. There is no apparent correlation between R{sub V} value and color excess E(B−V) in the range of interest, from 0.2 to 0.6 mag, or with specific visual extinction per kiloparsec, A{sub V}/d.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Auxiliary field diffusion Monte Carlo calculations of magnetic moments of light nuclei with chiral effective field theory interactions

We calculate the magnetic moments of light nuclei $( A < 20 )$ using the auxiliary field diffusion Monte Carlo method and local two- and three-nucleon forces with electromagnetic currents from chiral effective field theory. For all nuclei under consideration, we also calculate the ground-state energies and charge radii. We generally find a good agreement with experimental values for all of these observables. For the electromagnetic currents, we explore the impact of employing two different power counting schemes, and study theoretical uncertainties stemming from the truncation of the chiral expansion order by order for select nuclei within these two approaches. In conclusion, we find that it is crucial to employ consistent power counting schemes for interactions and currents to achieve a systematic order-by-order convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Status Report on Photothermal Radiometry Development

Thermal transport properties such as thermal conductivity and diffusivity, are critical parameters in the development of improved nuclear fuels and components. It is well known that thermal transport properties change temporally and spatially due to dynamic microstructural evolution as a result of neutron irradiation. Laser based photothermal techniques are an ideal candidate for thermal properties measurement of irradiated materials because they are non contact, non destructive and can be performed remotely. In this report, we present the progress in the development of a photothermal radiometry (PTR) system to measure the thermal diffusivity of solid samples. The PTR technique measures the local temperature response from a modulated heat source by detecting the sample radiation heat flux. Initially, a well developed thermal wave analytical model was fit to the data to obtain the thermal diffusivity. However, systematic errors from the diffraction effect and non linearity issues resulted in thermal diffusivity measurements higher than literature values. These effects have been studied, an improved analytical model has been developed, and the physics behind the phenomena are still under investigation. A procedure for optimizing measurements from the benchtop system has been provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Extending the capillary wave model to include the effect of bending rigidity: X-ray reflectivity and diffuse scattering

The surface roughness of a thin film at a liquid interface exhibits contributions of thermally excited fluctuations. This thermal roughness depends on temperature (𝑇), surface tension (𝛾), and elastic material properties, specifically the bending modulus (𝜅) of the film. A nonzero 𝜅 suppresses the thermal roughness at small length scales compared to an interface with zero 𝜅, as expressed by the power spectral density of the thermal roughness. The description of the x-ray scattering of the standard capillary wave model (CWM), which is valid for zero 𝜅, is extended to include the effect of 𝜅. The extended CWM (eCWM) provides a single analytical form for both the specular x-ray reflectivity (XRR) and the diffuse scattering around the specular reflection, and recovers the expression of the CWM at its zero 𝜅 limit. This theoretical approach enables the use of single-shot grazing incidence x-ray off-specular scattering (GIXOS) measurements for characterizing the structure of thin films on a liquid surface. The eCWM analysis approach decouples the thermal roughness factor from the surface scattering signal, providing direct access to the intrinsic surface-normal structure of the film and its bending modulus. Moreover, the eCWM facilitates the calculation of reflectivity at any desired resolution (pseudo-XRR approach). The transformation into pseudo-XRR provides the benefit of using widely available XRR software to perform GIXOS analysis. The extended range of the vertical scattering vector (𝑄 𝑧 ) available with the GIXOS pseudo-XRR approach allows for a higher spatial resolution than with conventional XRR. Experimental results are presented for various lipid systems, showing strong agreement between conventional specular XRR and pseudo-XRR methods. This agreement validates the proposed approach and highlights its utility for analyzing soft, thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach can be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of Framework Composition and NH 3 on the Diffusion of Cu + in Cu-CHA Catalysts Predicted by Machine-Learning Accelerated Molecular Dynamics

Cu-exchanged zeolites rely on mobile solvated Cu + cations for their catalytic activity, but the role of the framework composition in transport is not fully understood. Ab initio molecular dynamics simulations can provide quantitative atomistic insight but are too computationally expensive to explore large length and time scales or diverse compositions. We report a machine-learning interatomic potential that accurately reproduces ab initio results and effectively generalizes to allow multinanosecond simulations of large supercells and diverse chemical compositions. Biased and unbiased simulations of [Cu(NH 3 ) 2 ] + mobility show that aluminum pairing in eight-membered rings accelerates local hopping and demonstrate that increased NH 3 concentration enhances long-range diffusion. The probability of finding two [Cu(NH 3 ) 2 ] + complexes in the same cage, which is key for SCR-NOx reaction, increases with Cu content and Al content but does not correlate with the long-range mobility of Cu + . Supporting experimental evidence was obtained from reactivity tests of Cu-CHA catalysts with a controlled chemical composition.

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

Theory of Nonlinear Diffusion with a Physical Gapped Mode

In a system with one conserved charge the charge diffusion is modified by nonlinear self-interactions within an effective field theory (EFT) of diffusive fluctuations. We include the slowest ultraviolet (UV) mode, constructing a UV-regulated EFT. The relaxation time of this UV mode is protected from renormalization, as supported by experimental data in a bad metal system. Furthermore, the retarded density-density Green’s function acquires four branch points, eventually increasing the range of applicability. We discuss the fate of long time tails as well as implications for the quark gluon plasma. Published by the American Physical Society 2024

Physics↗
DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗
DOE OSTI2016

Quantitatively in Situ Imaging Silver Nanowire Hollowing Kinetics

We report the in-situ investigation of the morphological evolution of silver nanowires to hollow silver oxide nanotubes using transmission x-ray microscopy (TXM). Complex silver diffusion kinetics and hollowing process via the Kirkendall effect have been captured in real time. Further quantitative x-ray absorption analysis reveals the difference between the longitudinal and radial diffusions. In conclusion, the diffusion coefficient of silver in its oxide nanoshell is, for the first time, calculated to be 1.2 × 10 -13 cm 2 /s from the geometrical parameters extracted from the TXM images.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Finite Larmor radius effects on E x B weak turbulence transport

Transport of test particles in two-dimensional weak turbulence with waves propagating along the poloidal direction is studied using a reduced model. Finite Larmor radius (FLR) effects are included by gyroaveraging over one particle orbit. For low wave amplitudes the motion is mostly regular with particles trapped in the potential wells. As the amplitude increases the trajectories become chaotic and the Larmor radius modifies the orbits. For a thermal distribution of Finite Larmor radii the particle distribution function (PDF) is Gaussian for small $\unicode[STIX]{x1D70C}_{th}$ (thermal gyroradius) but becomes non-Gaussian for large $\unicode[STIX]{x1D70C}_{th}$ . Furthermore, the time scaling of transport is diffusive, as characterized by a linear dependence of the variance of the PDF with time. An explanation for this behaviour is presented that provides an expression for an effective diffusion coefficient and reproduces the numerical results for large wave amplitudes which implies generalized chaos. When a shear flow is added in the direction of wave propagation, a modified model is obtained that produces free-streaming particle trajectories in addition to trapped ones; these contribute to ballistic transport for low wave amplitude but produce super-ballistic transport in the chaotic regime. As in the previous case, the PDF is Gaussian for low $\unicode[STIX]{x1D70C}_{th}$ becoming non-Gaussian as it increases. The perpendicular transport presents the same behaviour as in the case with no flow but the diffusion is faster in the presence of the flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Final Project Report (Oct 2014-Dec 2017): Zeolite Membranes for Krypton/Xenon Separation from Spent Nuclear Fuel Reprocessing Off-Gas

The overall focus of this project is to develop and understand SAPO-34 zeolitic membranes that can separate mixtures of radioisotope krypton-85 and xenon released as off-gases during used nuclear fuel recycling. The primary advantage of separating 85 Kr from Xe is to reduce the volume of radioactive waste for storage. The second advantage is the revenue generated from the sale of high-purity Xe. Zeolite membranes are attractive because of their much lower energy requirements relative to cryogenic distillation, and their high resistance to radiation degradation. We report the detailed study of silicoaluminophosphate zeolite SAPO-34 materials and membranes for this application, due to hypothesized favorable molecular sieving properties. In the 3-year Mission Support project, we developed a novel, high-performance, low-energy intensity, lower-cost zeolite membrane process for Kr/Xe separation during SNF processing; and investigated the underlying molecular adsorption and transport processes in both ‘idealized’ and ‘realistic’ operating conditions to develop reliable synthesis-structure-property relationships for such membranes. Adsorption and diffusion measurements on SAPO-34 crystals indicate their potential for use in Kr-Xe separation membranes, but also highlight competing effects of adsorption and diffusion selectivity. SAPO-34 membranes synthesized on α-alumina substrates via steam-assisted conversion seeding and hydrothermal growth are characterized in detail, with Kr permeances 26 GPU and ideal Kr/Xe selectivities >20 at 298 K after thickness reduction. Post-synthesis cation exchange shows large (>50%) increases in selectivity at ambient or slight sub-ambient conditions. In addition, we confirm that SAPO-34 membrane is stable under radiation exposure and the impact of radiation exposure on membrane performance would not be substantial. We also successfully synthesized hollow-fiber SAPO-34 membranes with the same performance levels as the disk-type and tubular membranes. This important development will allow a very compact and low-cost Kr/Xe separation system. Finally, a detailed process calculation for techno-economic analysis was performed by integrating Maxwell-Stefan model into cross-flow membrane system, in order to estimate the required number of membrane stages and the total cost.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Physics of forced magnetic reconnection in coaxial helicity injection experiments in National Spherical Torus Experiment

We numerically examine the physics of fast flux closure in transient coaxial helicity injection (CHI) experiments in National Spherical Torus Experiment (NSTX). By performing resistive Magnetohydrodynamics (MHD) simulations with poloidal injector coil currents held constant in time, we find that closed flux surfaces are formed through forced magnetic reconnection. Through a local Sweet-Parker type reconnection with an elongated current sheet in the injector region, closed flux surfaces expand in the NSTX global domain. Simulations demonstrate outflows approaching poloidally Alfvénic flows and reconnection times consistent with the Sweet-Parker model. Critical requirements for magnetic reconnection and flux closure are studied in detail. These primary effects, which are magnetic diffusivity, injector flux, injector flux footprint width, and rate of injector voltage reduction, are simulated for transient CHI experiments. The relevant time scales for effective reconnection are τ{sub V}<τ{sub rec}≈τ{sub A}√(S)(1+Pm){sup 1/4}<τ{sub R}, where τ{sub V} is the time for the injector voltage reduction, τ{sub A} is the poloidal Alfvén transit time, τ{sub R} is the global resistive diffusion time, and Pm and S are Prandtl and Lundquist numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2015

Evaluation of Technetium Getters to Improve the Performance of Cast Stone

Cast Stone has been selected as the preferred waste form for solidification of aqueous secondary liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) process condensates and low-activity waste (LAW) melter off-gas caustic scrubber effluents. Cast Stone is also being evaluated as a supplemental immobilization technology to provide the necessary LAW treatment capacity to complete the Hanford tank waste cleanup mission in a timely and cost effective manner. One of the major radionuclides that Cast Stone has the potential to immobilize is technetium (Tc). The mechanism for immobilization is through the reduction of the highly mobile Tc(VII) species to the less mobile Tc(IV) species by the blast furnace slag (BFS) used in the Cast Stone formulation. Technetium immobilization through this method would be beneficial because Tc is one of the most difficult contaminants to address at the U.S. Department of Energy (DOE) Hanford Site due to its complex chemical behavior in tank waste, limited incorporation in mid- to high-temperature immobilization processes (vitrification, steam reformation, etc.), and high mobility in subsurface environments. In fact, the Tank Closure and Waste Management Environmental Impact Statement for the Hanford Site, Richland, Washington (TC&WM EIS) identifies technetium-99 ( 99 Tc) as one of the radioactive tank waste components contributing the most to the environmental impact associated with the cleanup of the Hanford Site. The TC&WM EIS, along with an earlier supplemental waste-form risk assessment, used a diffusion-limited release model to estimate the release of different contaminants from the WTP process waste forms. In both of these predictive modeling exercises, where effective diffusivities based on grout performance data available at the time, groundwater at the 100-m down-gradient well exceeded the allowable maximum permissible concentrations for 99 Tc. (900 pCi/L). Recent relatively short-term (63 day) leach tests conducted on both LAW and secondary waste Cast Stone monoliths indicated that 99 Tc diffusivities were at or near diffusivities where the groundwater at the 100-m down-gradient well would exceed the allowable maximum permissible 99 Tc concentrations. There is, therefore, a need and an opportunity to improve the retention of Tc in the Cast Stone waste form. One method to improve the performance of the Cast Stone waste form is through the addition of “getters” that selectively sequester Tc inside Cast Stone.

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

Governing Constraints of Chromium(VI) Formation from Chromium(III)-Bearing Minerals in Soils and Sediments

The potential for geogenic Cr(VI) contamination is vast, yet it is difficult to predict susceptible environments as groundwater Cr(VI) concentrations vary significantly, even within a given aquifer, due to physical and hydrologic heterogeneity. The physical structure of soils and sediments exerts a dominant control on Cr(VI) production by dictating the separation distance of reactive phases, the diffusion distance from Cr(VI) generation sites to advecting groundwater, and by influencing infiltration rates and porewater velocity. Here, we used a dual-pore domain model to investigate the relative control of these parameters on Cr(VI) production. The reaction distance between Cr(III)-bearing minerals and Mn oxides predominantly controls Cr(VI) export to advecting groundwater, while changes in diffusion distance between sites of Cr(VI) generation and advective flow channels generally have little impact on steady-state Cr(VI) concentrations. Changes in Cr(VI) diffusion distance can, however, increase the time required for groundwater Cr(VI) concentrations to reach a steady-state; thus, under fluctuating hydrologic and biogeochemical conditions, long diffusion distances still have the potential to suppress Cr(VI) supply to advecting water. Furthermore, we show that high porewater flow velocities effectively dilute Cr(VI) diffusing from soil/sediment aggregates, thus minimizing Cr(VI) concentrations relative to lower porewater velocities. The strong control that the physical/hydrologic parameters exert on Cr(VI) production appears to overwhelm the impact of Cr(III)-mineral solubility within soils and sediments.

Hausladen, Debra (ORCID:0000000342356294)↗