Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Breakthrough Curves”

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 109 records · Page 6

Adsorption characteristics of water vapor on gear-pellet and honeycomb-pellet types of adsorbents containing A-type zeolite

It is necessary to recover or process tritiated species that are extensively coexistent in nuclear fusion installations. A conventional way to recover tritium release to atmosphere is catalytic oxidation of tritiated species and adsorption of tritiated water vapor on adsorbents with high surface areas. Therefore, new adsorbents with low pressure loss and high surface areas need to be developed and utilized for such large-scale adsorption systems. In this study, attention was focused on new adsorbents, which are gear-type pellet MS5A adsorbent, gear-type pellet MS4A adsorbent and honeycomb-type pellet MS5A adsorbent. The adsorption characteristics of the new adsorbent were comparatively studied with conventional type of adsorbents (pellet-type MS5A adsorbent and pebble-type MS5A adsorbent), in terms of adsorption capacity, pressure loss and adsorption rate. It was found that the adsorption capacity of water vapor on the gear-type adsorbents is higher than that on a honeycomb-type adsorbent. The experimental breakthrough curves indicate that the adsorption rates of water vapor on gear-type and honeycomb-type adsorbents are smaller than that on conventional type adsorbents. Various adsorption models were also tested to correlate the experimental isotherms. It was found that the Langmuir-Freundlich model could properly correlate the experimental adsorption isotherms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of water saturation using gas phase partitioning tracers and time-lapse electrical conductivity measurements

Water saturation is an important indicator of contaminant distribution and plays a governing role in contaminant transport within the vadose zone. Understanding the water saturation distribution is critical for both remediation and contaminant flux monitoring in unsaturated environments. In this work we propose and demonstrate a method of remotely determining water saturation levels using gas phase partitioning tracers and time-lapse bulk electrical conductivity measurements. The theoretical development includes the partitioning chemistry for the tracers we demonstrate (ammonia and carbon dioxide), as well as a review of the petrophysical relationship governing how these tracers influence bulk conductivity. We also investigate methods of utilizing secondary information provided by electrical conductivity breakthrough magnitudes induced by the tracers. We test the method on clean, well characterized, intermediate-scale sand columns under controlled conditions. Results demonstrate the capability to predict partitioning coefficients and accurately monitor gas breakthrough curves along the length of the column according to the corresponding electrical conductivity response, leading to accurate water saturation estimates. This work is motivated by the need to develop effective characterization and monitoring techniques for contaminated deep vadose zone environments, and provides a proof-of-concept toward uniquely characterizing and monitoring water saturation levels at the field scale and in three-dimensions using electrical resistivity tomography.

30 DIRECT ENERGY CONVERSION↗

TRIPLICATE SODIUM IODIDE GAMMA RAY MONITORS FOR THE SMALL COLUMN ION EXCHANGE PROGRAM

This technical report contains recommendations from the Analytical Development (AD) organization of the Savannah River National Laboratory (SRNL) for a system of triplicate Sodium Iodide (NaI) detectors to be used to monitor Cesium-137 ({sup 137}Cs) content of the Decontaminated Salt Solution (DSS) output of the Small Column Ion Exchange (SCIX) process. These detectors need to be gain stabilized with respect to temperature shifts since they will be installed on top of Tank 41 at the Savannah River Site (SRS). This will be accomplished using NaI crystals doped with the alpha-emitting isotope, Americium-241({sup 241}Am). Two energy regions of the detector output will be monitored using single-channel analyzers (SCAs), the {sup 137}Cs full-energy {gamma}-ray peak and the {sup 241}Am alpha peak. The count rate in the gamma peak region will be proportional to the {sup 137}Cs content in the DSS output. The constant rate of alpha decay in the NaI crystal will be monitored and used as feedback to adjust the high voltage supply to the detector in response to temperature variation. An analysis of theoretical {sup 137}Cs breakthrough curves was used to estimate the gamma activity expected in the DSS output during a single iteration of the process. Count rates arising from the DSS and background sources were predicted using Microshield modeling software. The current plan for shielding the detectors within an enclosure with four-inch thick steel walls should allow the detectors to operate with the sensitivity required to perform these measurements. Calibration, testing, and maintenance requirements for the detector system are outlined as well. The purpose of SCIX is to remove and concentrate high-level radioisotopes from SRS salt waste resulting in two waste streams. The concentrated high-level waste containing {sup 137}Cs will be sent to the Defense Waste Processing Facility (DWPF) for vitrification and the low-level DSS will be sent to the Saltstone Production Facility (SPF) to be incorporated into grout.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigation into the Relationship between Heterogeneity and Heavy-Tailed Solute Transport (Final Report)

This project combines outcrop-scale heterogeneity characterization, laboratory experiments, and numerical simulations. The study is designed to test whether established dispersion theory accurately predicts the behavior of solute transport through heterogeneous media and to investigate the relationship between heterogeneity and the parameters that populate these models. The dispersion theory tested by this work is based upon the fractional advection-dispersion equation (fADE) model. Unlike most dispersion studies that develop a solute transport model by fitting the solute transport breakthrough curve, this project will explore the nature of the heterogeneous media to better understand the connection between the model parameters and the aquifer heterogeneity. Our work at the Colorado School of Mines was focused on the following questions: 1) What are the effects of multi-scale geologic variability on transport of conservative and reactive solutes? 2) Can those transport effects be accounted for by classical methods, and if not, can the nonlocal fractional-order equations provide better predictions? 3) Can the fractional-order equations be parameterized through a link to some simple observable geologic features? 4) Are the classical equations of transport and reaction sufficient? 5) What is the effect of anomalous transport on chemical reaction in groundwater systems? The work is predicated on the observation that upscaled transport is defined by loss of information, or spatio-temporal averaging. This averaging tends to make the transport laws such as Fick's 2nd-order diffusion equation similar to central limit theory. The fractional-order advection-dispersion equations rely on limit theory for heavy-tailed random motion that has some diverging moments. The equations predict larger tails of a plume in space and/or time than those predicted by the classical 2nd-order advection-dispersion equation. The heavy tails are often seen in plumes at field sites.

58 GEOSCIENCES↗

OSPREY Model

The absence of industrial scale nuclear fuel reprocessing in the U.S. has precluded the necessary driver for developing the advanced simulation capability now prevalent in so many other countries. Thus, it is essential to model complex series of unit operations to simulate, understand, and predict inherent transient behavior and feedback loops. A capability of accurately simulating the dynamic behavior of advanced fuel cycle separation processes will provide substantial cost savings and many technical benefits. The specific fuel cycle separation process discussed in this report is the off-gas treatment system. The off-gas separation consists of a series of scrubbers and adsorption beds to capture constituents of interest. Dynamic models are being developed to simulate each unit operation involved so each unit operation can be used as a stand-alone model and in series with multiple others. Currently, an adsorption model has been developed within Multi-physics Object Oriented Simulation Environment (MOOSE) developed at the Idaho National Laboratory (INL). Off-gas Separation and REcoverY (OSPREY) models the adsorption of off-gas constituents for dispersed plug flow in a packed bed under non-isothermal and non-isobaric conditions. Inputs to the model include gas, sorbent, and column properties, equilibrium and kinetic data, and inlet conditions. The simulation outputs component concentrations along the column length as a function of time from which breakthrough data is obtained. The breakthrough data can be used to determine bed capacity, which in turn can be used to size columns. It also outputs temperature along the column length as a function of time and pressure drop along the column length. Experimental data and parameters were input into the adsorption model to develop models specific for krypton adsorption. The same can be done for iodine, xenon, and tritium. The model will be validated with experimental breakthrough curves. Customers will be given access to OSPREY to used and evaluate the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of Microbially Induced Calcite Precipitation for the Immobilization of Strontium-90 Using a Variable Velocity Streamtube Ensemble (Final Technical Report)

Biogeochemical modeling using PHREEQC2 and a streamtube ensemble approach is utilized to understand a well-to-well subsurface treatment system at the Vadose Zone Research Park (VZRP) near Idaho Falls, Idaho. Treatment involves in situ microbially-mediated ureolysis to induce calcite precipitation for the immobilization of strontium-90. PHREEQC2 is utilized to model the kinetically-controlled ureolysis and consequent calcite precipitation. Reaction kinetics, equilibrium phases, and cation exchange are used within PHREEQC2 to track pH and levels of calcium, ammonium, urea, and calcite precipitation over time, within a series of one-dimensional advective-dispersive transport paths creating a streamtube ensemble representation of the well-to-well transport. An understanding of the impact of physical heterogeneities within this radial flowfield is critical for remediation design; we address this via the streamtube approach: instead of depicting spatial extents of solutes in the subsurface we focus on their arrival distribution at the control well(s). Traditionally, each streamtube maintains uniform velocity; however in radial flow in homogeneous media, the velocity within any given streamtube is spatially-variable in a common way, being highest at the input and output wells and approaching a minimum at the midpoint between the wells. This idealized velocity variability is of significance in the case of ureolytically driven calcite precipitation. Streamtube velocity patterns for any particular configuration of injection and withdrawal wells are available as explicit calculations from potential theory, and also from particle tracking programs. To approximate the actual spatial distribution of velocity along streamtubes, we assume idealized radial non-uniform velocity associated with homogeneous media. This is implemented in PHREEQC2 via a non-uniform spatial discretization within each streamtube that honors both the streamtube’s travel time and the idealized “fast-slow-fast” pattern of non-uniform velocity along the streamline. Breakthrough curves produced by each simulation are weighted by the path-respective flux fractions (obtained by deconvolution of tracer tests conducted at the VZRP) to obtain the flux-average of flow contributions to the observation well.

54 ENVIRONMENTAL SCIENCES↗

Investigation into the Relationship between Heterogeneity and Heavy-Tailed Solute Transport (Final Technical Report)

The objective of this project was to characterize the influence that naturally complex geologic media has on anomalous dispersion and to determine if the nature of dispersion can be estimated from the underlying heterogeneous media. The UNM portion of this project was to provide detailed representations of aquifer heterogeneity through producing highly-resolved models of outcrop analogs to aquifer materials. This project combined outcrop-scale heterogeneity characterization (conducted at the University of New Mexico), laboratory experiments (conducted at Sandia National Laboratory), and numerical simulations (conducted at Sandia National Laboratory and Colorado School of Mines). The study was designed to test whether established dispersion theory accurately predicts the behavior of solute transport through heterogeneous media and to investigate the relationship between heterogeneity and the parameters that populate these models. The dispersion theory tested by this work was based upon the fractional advection-dispersion equation (fADE) model. Unlike most dispersion studies that develop a solute transport model by fitting the solute transport breakthrough curve, this project explored the nature of the heterogeneous media to better understand the connection between the model parameters and the aquifer heterogeneity. We also evaluated methods for simulating the heterogeneity to see whether these approaches (e.g., geostatistical) could reasonably replicate realistic heterogeneity. The UNM portion of this study focused on capturing realistic geologic heterogeneity of aquifer analogs using advanced outcrop mapping methods.

54 ENVIRONMENTAL SCIENCES↗

Scale-Dependent Fracture-Matrix Interactions And Their Impact on Radionuclide Transport - Final Report

Matrix diffusion and adsorption within a rock matrix are widely regarded as important mechanisms for retarding the transport of radionuclides and other solutes in fractured rock (e.g., Neretnieks, 1980; Tang et al., 1981; Maloszewski and Zuber, 1985; Novakowski and Lapcevic, 1994; Jardine et al., 1999; Zhou and Xie, 2003; Reimus et al., 2003a,b). When remediation options are being evaluated for old sources of contamination, where a large fraction of contaminants reside within the rock matrix, slow diffusion out of the matrix greatly increases the difficulty and timeframe of remediation. Estimating the rates of solute exchange between fractures and the adjacent rock matrix is a critical factor in quantifying immobilization and/or remobilization of DOE-relevant contaminants within the subsurface. In principle, the most rigorous approach to modeling solute transport with fracture-matrix interaction would be based on local-scale coupled advection-diffusion/dispersion equations for the rock matrix and in discrete fractures that comprise the fracture network (Discrete Fracture Network and Matrix approach, hereinafter referred to as DFNM approach), fully resolving aperture variability in fractures and matrix property heterogeneity. However, such approaches are computationally demanding, and thus, many predictive models rely upon simplified models. These models typically idealize fracture rock masses as a single fracture or system of parallel fractures interacting with slabs of porous matrix or as a mobile-immobile or multi-rate mass transfer system. These idealizations provide tractable approaches for interpreting tracer tests and predicting contaminant mobility, but rely upon a fitted effective matrix diffusivity or mass-transfer coefficients. However, because these fitted parameters are based upon simplified conceptual models, their effectiveness at predicting long-term transport processes remains uncertain. Evidence of scale dependence of effective matrix diffusion coefficients obtained from tracer tests highlights this point and suggests that the underlying mechanisms and relationship between rock and fracture properties are not fully understood in large complex fracture networks. In this project, we developed a high-resolution DFN model of solute transport in fracture networks to explore and quantify the mechanisms that control transport in complex fracture networks and how these may give rise to observed scale-dependent matrix diffusion coefficients. Results demonstrate that small scale heterogeneity in the flow field caused by local aperture variability within individual fractures can lead to long-tailed breakthrough curves indicative of matrix diffusion, even in the absence of interactions with the fracture matrix. Furthermore, the temporal and spatial scale dependence of these processes highlights the inability of short-term tracer tests to estimate transport parameters that will control long-term fate and transport of contaminants in fractured aquifers.

58 GEOSCIENCES↗

Adsorptive Separation and Sequestration of Krypton, I and C14 on Diamond Nanoparticles

The objective of this research proposal was to address the separation and sequestration of Kr and I from each other using nano-sized diamond particles and retaining these in diamond until they decay to the background level or can be used as a byproduct. Following removal of Kr and I, an adsorbent will be used to adsorb and store CO 2 from the CO 2 rich stream. A Field Enhanced Diffusion with Optical Activation (FEDOA-a large scale process that takes advantage of thermal, electrical, and optical activation to enhance the diffusion of an element into diamond structure) was used to load Kr and I on micron or nano sized particles having a larger relative surface area. The diamond particles can be further increased by doping it with boron followed by irradiation in a neutron flux. Previous studies showed that the hydrogen storage capacity could be increased significantly by using boron-doped irradiated diamond particles. Diamond powders were irradiated for a longer time by placing them in a quartz tube. The surface area was measured using a Quantachrome Autosorb system. No significant increase in the surface area was observed. Total surface area was about 1.7 m 2 /g. This suggests the existence of very minimal pores. Interestingly it showed hysteresis upon desorption. A reason for this may be strong interaction between the surface and the nitrogen molecules. Adsorption runs at higher temperatures did not show any adsorption of krypton on diamond. Use of a GC with HID detector to determine the adsorption capacity from the breakthrough curves was attempted, but experimental difficulties were encountered.

36 MATERIALS SCIENCE↗

Multi-component testing using HZ-PAN and AgZ-PAN Sorbents for OSPREY Model validation

In efforts to further develop the capability of the Off-gas SeParation and RecoverY (OSPREY) model, multi-component tests were completed using both HZ-PAN and AgZ-PAN sorbents. The primary purpose of this effort was to obtain multi-component xenon and krypton capacities for comparison to future OSPREY predicted multi-component capacities using previously acquired Langmuir equilibrium parameters determined from single component isotherms. Experimental capacities were determined for each sorbent using two feed gas compositions of 1000 ppmv xenon and 150 ppmv krypton in either a helium or air balance. Test temperatures were consistently held at 220 K and the gas flowrate was 50 sccm. Capacities were calculated from breakthrough curves using TableCurve® 2D software by Jandel Scientific. The HZ-PAN sorbent was tested in the custom designed cryostat while the AgZ-PAN was tested in a newly installed cooling apparatus. Previous modeling validation efforts indicated the OSPREY model can be used to effectively predict single component xenon and krypton capacities for both engineered form sorbents. Results indicated good agreement with the experimental and predicted capacity values for both krypton and xenon on the sorbents. Overall, the model predicted slightly elevated capacities for both gases which can be partially attributed to the estimation of the parameters and the uncertainty associated with the experimental measurements. Currently, OSPREY is configured such that one species adsorbs and one does not (i.e. krypton in helium). Modification of OSPREY code is currently being performed to incorporate multiple adsorbing species and non-ideal interactions of gas phase species with the sorbent and adsorbed phases. Once these modifications are complete, the sorbent capacities determined in the present work will be used to validate OSPREY multicomponent adsorption predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adsorption Properties of Lignin-derived Activated Carbon Fibers (LACF)

The object of this CRADA project between Oak Ridge National Laboratory (ORNL) and United Technologies Research Center (UTRC) is the characterization of lignin-derived activated carbon fibers (LACF) and determination of their adsorption properties for volatile organic compounds (VOC). Carbon fibers from lignin raw materials were manufactured at Oak Ridge National Laboratory (ORNL) using the technology previously developed at ORNL. These fibers were physically activated at ORNL using various activation conditions, and their surface area and pore-size distribution were characterized by gas adsorption. Based on these properties, ORNL did down-select five differently activated LACF materials that were delivered to UTRC for measurement of VOC adsorption properties. UTRC used standard techniques based on breakthrough curves to measure and determine the adsorption properties of indoor air pollutants (IAP) - namely formaldehyde and carbon dioxide - and to verify the extent of saturated fiber regenerability by thermal treatments. The results are summarized as follows: (1) ORNL demonstrated that physical activation of lignin-derived carbon fibers can be tailored to obtain LACF with surface areas and pore size distributions matching the properties of activated carbon fibers obtained from more expensive, fossil-fuel precursors; (2) UTRC investigated the LACF potential for use in air cleaning applications currently pursued by UTRC, such as building ventilation, and demonstrated their regenerability for CO2 and formaldehyde, (3) Both partners agree that LACF have potential for possible use in air cleaning applications.

36 MATERIALS SCIENCE↗

Removal of High Specific Activity Fission Products from Uranyl Sulfate Waste Solutions

The Savannah River National Laboratory (SRNL) is currently providing support to SHINE Medical Technologies (SHINE) which plans to deploy a low energy, accelerator-based neutron source to fission low enriched U in a uranyl sulfate target solution for 99 Mo production. The 99 Mo is initially separated from the fission products and target solution by an extraction column. Subsequent washing of the column will generate waste solutions containing residual U and fission product activity. A small number of high specific activity fission products (e.g., 90 Sr, 137 Cs, and 144 Ce) in these streams will likely control the classification of the low level waste (LLW). If a sufficient amount of the high specific activity isotopes are separated from the SHINE waste streams and concentrated into a waste form, it would be possible to treat a majority of the wash solutions from the column operations as a lower class of LLW (Class A versus Class B or C or Class B versus Class C). The high specific activity fission product elements could then be dispositioned as a much smaller volume of waste rather than requiring the disposal of the entire waste stream at the higher waste classification. The Savannah River Site (SRS) has experience with using monosodium titanate (MST) and crystalline silicotitanate (CST) to remove Cs and Sr from high salt content waste solutions generated during the reprocessing of nuclear fuels and targets. Both of these materials have worked very well for their intended purposes at the SRS where the fission product elements are removed from highly alkaline waste. On the other hand, SHINE waste streams from the extraction column contain H 2 SO 4 which makes the solution acidic. Additionally, the SRS waste streams do not contain other fission product elements such as transition metals or lanthanides because they precipitate upon neutralization of the SRS waste and are not present in the supernate which is dispositioned as LLW following treatment. As such, there are inherent differences between SHINE and SRS waste treatment strategies. Savannah River National Laboratory was tasked with performing scoping studies to see if MST and CST would remove Sr, Cs, and Ce from an acidic mixed metal simulant solution. Batch contact experiments were performed using MST and two CST type materials. The MST material is a 15 wt % powder in 0.15 M NaOH slurry. The MST showed low adsorption for elements of interest from acidic solution. Furthermore, the powder size makes MST non-ideal for column operations. A CST IE-911 ion exchange material had high Cs adsorption, moderate Sr, and marginal Ce adsorption. Based on adsorption of all species, the ion exchange capacity was found to be 0.032 meq/mL. A bench-top column experiment to measure elemental breakthrough curves was performed using CST IE-911 where chromatographic separations of the mixed simulant were expected to occur. While most elements behaved as expected, the lanthanide series, containing Ce, broke though the column earlier than expected. The second CST material, CST R9120, displayed high adsorption for all elements in the acidic mixed simulant solutions in a batch contact study, and had a calculated loading capacity of 0.091 meq/mL. Future studies to develop a waste treatment flowsheet should focus on CST R9120 to treat SHINE waste solutions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

Multi-Scale Microfluidics for Transport in Shale Fabric

We develop a microfluidic experimental platform to study solute transport in multi-scale fracture networks with a disparity of spatial scales ranging between two and five orders of magnitude. Using the experimental scaling relationship observed in Marcellus shales between fracture aperture and frequency, the microfluidic design of the fracture network spans all length scales from the micron (1 μ) to the dm (10 dm). This intentional `tyranny of scales’ in the design, a determining feature of shale fabric, introduces unique complexities during microchip fabrication, microfluidic flow-through experiments, imaging, data acquisition and interpretation. Here, we establish best practices to achieve a reliable experimental protocol, critical for reproducible studies involving multi-scale physical micromodels spanning from the Darcy- to the pore-scale (dm to μm). With this protocol, two fracture networks are created: a macrofracture network with fracture apertures between 5 and 500 μm and a microfracture network with fracture apertures between 1 and 500 μm. The latter includes the addition of 1 μm ‘microfractures’, at a bearing of 55°, to the backbone of the former. Comparative analysis of the breakthrough curves measured at corresponding locations along primary, secondary and tertiary fractures in both models allows one to assess the scale and the conditions at which microfractures may impact passive transport.

58 GEOSCIENCES↗

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗

Reactive transport modeling for supporting climate resilience at groundwater contamination sites

Abstract. Climate resilience is an emerging issue at contaminated sites and hazardous waste sites, since projected climate shifts (e.g., increased/decreased precipitation) and extreme events (e.g., flooding, drought) could affect ongoing remediation or closure strategies. In this study, we develop a reactive transport model (Amanzi) for radionuclides (uranium, tritium, and others) and evaluate how different scenarios under climate change will influence the contaminant plume conditions and groundwater well concentrations. We demonstrate our approach using a two-dimensional (2D) reactive transport model for the Savannah River Site F-Area, including mineral reaction and sorption processes. Different recharge scenarios are considered by perturbing the infiltration rate from the base case as well as considering cap-failure and climate projection scenarios. We also evaluate the uranium and nitrate concentration ratios between scenarios and the base case to isolate the sorption effects with changing recharge rates. The modeling results indicate that the competing effects of dilution and remobilization significantly influence pH, thus changing the sorption of uranium. At the maximum concentration on the breakthrough curve, higher aqueous uranium concentration implies that sorption is reduced with lower pH due to remobilization. To better evaluate the climate change impacts in the future, we develop the workflow to include the downscaled CMIP5 (Coupled Model Intercomparison Project) climate projection data in the reactive transport model and evaluate how residual contamination evolves through 2100 under four climate Representative Concentration Pathway (RCP) scenarios. The integration of climate modeling data and hydrogeochemistry models enables us to quantify the climate change impacts, assess which impacts need to be planned for, and therefore assist climate resiliency efforts and help guide site management.

54 ENVIRONMENTAL SCIENCES↗

Biochar alters hydraulic conductivity and impacts nutrient leaching in two agricultural soils

Biochar is purported to provide agricultural benefits when added to the soil, through changes in saturated hydraulic conductivity (K sat ) and increased nutrient retention through chemical or physical means. Despite increased interest and investigation, there remains uncertainty regarding the ability of biochar to deliver these agronomic benefits due to differences in biochar feedstock, production method, production temperature, and soil texture. In this project, a suite of experiments was carried out using biochars of diverse feedstocks and production temperatures, in order to determine the biochar parameters which may optimize agricultural benefits. Sorption experiments were performed with seven distinct biochars to determine sorption efficiencies for ammonium and nitrate. Only one biochar effectively retained nitrate, while all biochars bound ammonium. The three biochars with the highest binding capacities (produced from almond shell at 500 and 800 °C (AS500 and AS800) and softwood at 500 °C (SW500)) were chosen for column experiments. Biochars were amended to a sandy loam and a silt loam at 0 % and 2 % (w/w), and K sat was measured. Biochars reduced K sat in both soils by 64 %–80 %, with the exception of AS800, which increased K sat by 98 % in the silt loam. Breakthrough curves for nitrate and ammonium, as well as leachate nutrient concentration, were also measured in the sandy loam columns. All biochars significantly decreased the quantity of ammonium in the leachate, by 22 % to 78 %, and slowed its movement through the soil profile. Ammonium retention was linked to high cation exchange capacity and a high oxygen-to-carbon ratio, indicating that the primary control of ammonium retention in biochar-amended soils is the chemical affinity between biochar surfaces and ammonium. Biochars had little to no effect on the timing of nitrate release, and only SW500 decreased total quantity, by 27 % to 36 %. The ability of biochar to retain nitrate may be linked to high micropore specific surface area, suggesting a physical entrapment rather than a chemical binding. Together, this work sheds new light on the combined chemical and physical means by which biochar may alter soils to impact nutrient leaching and hydraulic conductivity for agricultural production.

54 ENVIRONMENTAL SCIENCES↗