Critical review of coupled flux formulations for clay membranes based on nonequilibrium thermodynamics

Extensive research conducted over the past several decades has indicated that semipermeable membrane behavior (i.e., the ability of a porous medium to restrict the passage of solutes) may have a significant influence on solute migration through a wide variety of clay-rich soils, including both natural clay formations (aquitards, aquicludes) and engineered clay barriers (e.g., landfill liners and vertical cutoff walls). Restricted solute migration through clay membranes generally has been described using coupled flux formulations based on nonequilibrium (irreversible) thermodynamics. However, these formulations have differed depending on the assumptions inherent in the theoretical development, resulting in some confusion regarding the applicability of the formulations. Accordingly, a critical review of coupled flux formulations for liquid, current, and solutes through a semipermeable clay membrane under isothermal conditions is undertaken with the goals of explicitly resolving differences among the formulations and illustrating the significance of the differences from theoretical and practical perspectives. Formulations based on single-solute systems (i.e., uncharged solute), single-salt systems, and general systems containing multiple cations or anions are presented. Also, expressions relating the phenomenological coefficients in the coupled flux equations to relevant soil properties (e.g., hydraulic conductivity and effective diffusion coefficient) are summarized for each system. A major difference in the formulations is shown to exist depending on whether counter diffusion or salt diffusion is assumed. This difference between counter and salt diffusion is shown to affect the interpretation of values for the effective diffusion coefficient in a clay membrane based on previously published experimental data. Solute transport theories based on both counter and salt diffusion then are used to re-evaluate previously published column test data for the same clay membrane. The results indicate that, despite the theoretical inconsistency between the counter-diffusion assumption and the salt-diffusion conditions of the experiments, the predictive ability of solute transport theory based on the assumption of counter diffusion is not significantly different from that based on the assumption of salt diffusion, provided that the input parameters used in each theory are derived under the same assumption inherent in the theory. Nonetheless, salt-diffusion theory is fundamentally correct and, therefore, is more appropriate for problems involving salt diffusion in clay membranes. Finally, the fact that solute diffusion cannot occur in an ideal or perfect membrane is not explicitly captured in any of the theoretical expressions for total solute flux in clay membranes, but rather is generally accounted for via inclusion of an effective porosity, n(e), or a restrictive tortuosity factor, tau(r), in the formulation of Fick's first law for diffusion. Both n(e) and tau(r) have been correlated as a linear function of membrane efficiency. This linear correlation is supported theoretically by pore-scale modeling of solid-liquid interactions, but experimental support is limited. Additional data are needed to bolster the validity of the linear correlation for clay membranes. Copyright 2012 Elsevier B.V. All rights reserved.

Critical Review of Coupled Flux Formulations for Clay Membranes Based on Nonequilibrium Thermodynamics

Extensive research conducted over the past several decades has indicated that semipermeable membrane behavior (i.e., the ability of a porous medium to restrict the passage of solutes) may have a significant influence on solute migration through a wide variety of clay-rich soils, including both natural clay formations (aquitards, aquicludes) and engineered clay barriers (e.g., landfill liners and vertical cutoff walls). Restricted solute migration through clay membranes generally has been described using coupled flux formulations based on nonequilibrium (irreversible) thermodynamics. However, these formulations have differed depending on the assumptions inherent in the theoretical development, resulting in some confusion regarding the applicability of the formulations. Accordingly, a critical review of coupled flux formulations for liquid, current, and solutes through a semipermeable clay membrane under isothermal conditions is undertaken with the goals of explicitly resolving differences among the formulations and illustrating the significance of the differences from theoretical and practical perspectives. Formulations based on single-solute systems (i.e., uncharged solute), single-salt systems, and general systems containing multiple cations or anions are presented. Also, expressions relating the phenomenological coefficients in the coupled flux equations to relevant soil properties (e.g., hydraulic conductivity and effective diffusion coefficient) are summarized for each system. A major difference in the formulations is shown to exist depending on whether counter diffusion or salt diffusion is assumed. This difference between counter and salt diffusion is shown to affect the interpretation of values for the effective diffusion coefficient in a clay membrane based on previously published experimental data. Solute transport theories based on both counter and salt diffusion then are used to re-evaluate previously published column test data for the same clay membrane. The results indicate that, despite the theoretical inconsistency between the counter-diffusion assumption and the salt-diffusion conditions of the experiments, the predictive ability of solute transport theory based on the assumption of counter diffusion is not significantly different from that based on the assumption of salt diffusion, provided that the input parameters used in each theory are derived under the same assumption inherent in the theory. Nonetheless, salt-diffusion theory is fundamentally correct and, therefore, is more appropriate for problems involving salt diffusion in clay membranes. Finally, the fact that solute diffusion cannot occur in an ideal or perfect membrane is not explicitly captured in any of the theoretical expressions for total solute flux in clay membranes, but rather is generally accounted for via inclusion of an effective porosity, ne, or a restrictive tortuosity factor, tr, in the formulation of Fick's first law for diffusion. Both ne and tr have been correlated as a linear function of membrane efficiency. This linear correlation is supported theoretically by pore-scale modeling of solid-liquid interactions, but experimental support is limited. Additional data are needed to bolster the validity of the linear correlation for clay membranes.

Investigation of Solid-State Shear Pulverization Processing Parameters for Polyethylene-Clay Nanocomposites

Investigates multiple processing parameters, includingpolymer type, filler type, processing technique, severity of SSSP (Solid-state shear pulverization)processing, and postprocessing, of SSSP. HDPE and LLDPE polymers with pristine clay and organo-clay samples are explored. Effects on crystallization, high-temperature behavior, mechanicalproperties, and gas barrier properties are examined. Thermal, mechanical, and morphological characterization is conducted to determine polymer/filler compatibility and superior processing methods for the polymer-clay nanocomposites.

Hoisting C Structures Into Clay In Device Drivers

This project addresses the unreliability of operating system code, in particular in device drivers. Device driver software is the interface between the operating system and the device's hardware. Device drivers are written in low level code, making them difficult to understand. Almost all device drivers are written in the programming language C which allows for direct manipulation of memory. Due to the complexity of manual movement of data, most mistakes in operating systems occur in device driver code. The programming language Clay can be used to check device driver code at compile-time. Clay does most of its error checking statically to minimize the overhead of run-time checks in order to stay competitive with C's performance time. The Clay compiler can detect a lot more types of errors than the C compiler like buffer overflows, kernel stack overflows, NULL pointer uses, freed memory uses, and aliasing errors. Clay code that successfully compiles is guaranteed to run without failing on errors that Clay can detect. Even though C is unsafe, currently most device drivers are written in it. Not only are device drivers the part of the operating system most likely to fail, they also are the largest part of the operating system. As rewriting every existing device driver in Clay by hand would be impractical, this thesis is part of a project to automate translation of existing drivers from C to Clay. Although C and Clay both allow low level manipulation of data and fill the same niche for developing low level code, they have different syntax, type systems, and paradigms. This paper explores how C can be translated into Clay. It identifies what part of C device drivers cannot be translated into Clay and what information drivers in Clay will require that C cannot provide. It also explains how these translations will occur by explaining how each C structure is represented in the compiler and how these structures are changed to represent a Clay structure.

Monitoring Stormwater Redistribution Into Low Resistivity Soils Using Noninvasive Geophysical Techniques

Water held in the unsaturated zone is important for agriculture and construction and is replenished by infiltrating rainwater. Monitoring the soil water content of clay soils using ground-penetrating radar (GPR) has not been researched, as clay soils cause attenuation of GPR signal. In this study, GPR common-midpoint soundings (CMPs) are used in the clayey soils of the Miller Run floodplain to monitor changes in the soil water content (SWC) before and after rainfall events. GPR accomplishes this task because increases in water content will increase the dielectric constant of the subsurface material, and decrease the velocity of the GPR wave. Using an empirical relationship between dielectric constant and SWC, the Topp relation, we are able to calculate a SWC from these velocity measurements. Non-invasive electromagnetics, resistivity, and seismic were performed, and from these surveys, the layering at the field site was delineated. EM characterized the horizontal variation of the soil, allowing us to target the most clay rich area. At the CMP location, resistivity indicates the vertical structure of the subsurface consists of a 40 cm thick layer with a resistivity of 100 ohm*m. Between 40 cm and 1.5 m is a layer with a resistivity of 40 ohm*m. The thickness estimates were confirmed with invasive auger and trenching methods away from the CMP location. GPR CMPs were collected relative to a July 2013 and September 2013 storm. The velocity observations from the CMPs had a precision of +/- 0.001 m/ns as assessed by repeat analysis. In the case of both storms, the GPR data showed the expected relationship between the rainstorms and calculated SWC, with the SWC increasing sharply after the rainstorm and decreasing as time passed. We compared these data to auger core samples collected at the same time as the CMPs were taken, and the volumetric analysis of the cores confirmed the trend seen in the GPR, with SWC values between 3 and 5 percent lower than the GPR estimates. Our data shows that we can, with good precision, monitor changes in the SWC of conductive soils in response to rainfall events, despite the attenuation induced by the clay.

Microscopic and Spectroscopic Methods for Probing the Clay Water Interface

Clay minerals have a fundamental importance in many processes in soils and sediments such as the bioavailability of nutrients, water retention, the adsorption of common pollutants, and the formation of an impermeable barrier upon swelling. Many of the properties of clay minerals are due to the unique environment present at the clay mineral/water interface. Traditional techniques such as X-ray diffraction (XRD) and absorption isotherms have provided a wealth of information about this interface but have suffered from limitations. The methods and results presented herein are designed to yield new experimental information about the clay mineral/water interface.A new method of studying the swelling dynamics of clay minerals was developed using in situ atomic force microscopy (AFM). The preliminary results presented here demonstrate that this technique allows one to study individual clay mineral unit layers, explore the natural heterogeneities of samples, and monitor swelling dynamics of clay minerals in real time. Cation exchange experiments were conducted monitoring the swelling change of individual nontronite quasi-crystals as the chemical composition of the surrounding environment was manipulated several times. A proof of concept study has shown that the changes in swelling are from the exchange of interlayer cations and not from the mechanical force of replacing the solution in the fluid cell. A series of attenuated total internal reflection Fourier transform infrared spectroscopy (ATR-FTIR) experiments were performed to gain a better understanding of the organization of water within the interlayer region of two Fe-bearing clay minerals. These experiments made use of the Subtractive Kramers-Kronig (SKK) Transform and the calculation of difference spectra to obtain information about interfacial water hidden within the absorption bands of bulk water. The results indicate that the reduction of structural iron disrupts the organization of water around a strongly hydrated cation such as sodium as the cation transitions from an outer-sphere complex with the mineral surface to an inner-sphere complex. In the case of a less strongly hydrated cation such as potassium, reduction of structural iron actually increases the ordering of water molecules at the mineral surface. These effects were only noticed with the reduction of iron in the tetrahedral sheet close to the basal surface where the increased charge density is localized closer to the cations in the interlayer.

Hydrologic and Water Quality Assessment of Miller Run: A Study of Bucknell University’s Impact

Human development causes degradation of stream ecosystems due to impacts on channel morphology, hydrology, and water quality. Urbanization, the second leading cause of stream impairment, increases the amount of impervious surface cover, thus reducing infiltration and increasing surface runoff of precipitation, which ultimately affects stream hydrologic process and aquatic biodiversity. The main objective of this study was to assess the overall health of Miller Run, a small tributary of the Bull Run and Susquehanna River watersheds, through an integrative hydrologic and water quality approach in order to determine the degree of Bucknell University’s impact on the stream. Hydrologic conditions, including stage and discharge, and water quality conditions, including total suspended solids, ion, nutrient, and dissolved metal concentrations, specific conductivity, pH, and temperature, were measured and evaluated at two sampling sites (upstream and downstream of Bucknell’s main campus) during various rain events from September 2007 to March 2008. The primary focus of the stream analysis was based on one main rain event on 26 February 2008. The results provided evidence that Miller Run is impacted by Bucknell’s campus. From a hydrologic perspective, the stream’s hydrograph showed the exact opposite pattern of what would be expected from a ‘normal’ stream. Miller run had a flashier downstream hydrograph and a broader upstream hydrograph, which was more than likely due to the increased amount of impervious surface cover throughout the downstream half of the watershed. From a water quality perspective, sediment loads increased at a faster rate and were significantly higher downstream compared to upstream. These elevated sediment concentrations were probably the combined result of sediment runoff from upstream and downstream construction sites that were being developed over the course of the study. Sodium, chloride, and potassium concentrations, in addition to specific conductivity, also significantly increased downstream of Bucknell’s campus due to the runoff of road salts. Calcium and magnesium concentrations did not appear to be impacted by urbanization, although they did demonstrate a significant dilution effect downstream. The downstream site was not directly affected by elevated nitrate concentrations; however, soluble reactive phosphorus concentrations tended to increase downstream and ammonium concentrations significantly peaked partway through the rain event downstream. These patterns suggest that Miller Run may be impacted by nutrient runoff from the golf course, athletic fields, and/or fertilizers applications on the main campus. Dissolved manganese and iron concentrations also appeared to slightly increase downstream, demonstrating the affect of urban runoff from roads and parking lots. pH and temperature both decreased farther downstream, but neither showed a significant impact of urbanization. More studies are necessary to determine how Miller Run responds to changes in season, climate, precipitation intensity, and land-use. This study represents the base-line analysis of Miller Run’s current hydrologic and water quality conditions; based on these initial findings, Bucknell should strongly consider modifications to improve storm water management practices and to reduce the campus’s overall impact on the stream in order to enhance and preserve the integrity of its natural water resources.