Search NASA⌕ Search

SEARCH · Search NASA

Results for “Process Control”

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 739 records · Page 41

The evolution of the Human Systems and Simulation Laboratory in nuclear power research

The events at Three Mile Island in the United States brought about fundamental changes in the ways that simulation would be used in nuclear operations. The need for research simulators was identified to scientifically study human-centered risk and make recommendations for process control system designs. This paper documents the human factors research conducted at the Human Systems and Simulation Laboratory (HSSL) since its inception in 2010 at Idaho National Laboratory. The facility’s primary purposes are to provide support to utilities for system upgrades and to validate modernized control room concepts. In the last decade, however, as nuclear industry needs have evolved, so too have the purposes of the HSSL. Thus, beyond control room modernization, human factors researchers have evaluated the security of nuclear infrastructure from cyber adversaries and evaluated human-in-the-loop simulations for joint operations with an integrated hydrogen generation plant. Lastly, our review presents research using human reliability analysis techniques with data collected from HSSL-based studies and concludes with potential future directions for the HSSL, including severe accident management and advanced control room technologies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Bulk Attenuation in the Earth’s Mantle Due to Phase Transitions

The calculated bulk attenuation due to phase transition from mineral physics data is reported here. With relaxation time less than 1 s, the calculated value for a pyrolite mantle is consistent with the inverted bulk attenuation of the upper mantle from seismic observations. The two important mechanisms of phase transitions, diffusion-controlled and kinetics-controlled, have different relaxation time as indicated by the models here. The diffusion controlled is more likely to contribute to the observed seismic bulk attenuation than the kinetic-controlled process based on the available diffusivity and kinetics data. The correlation between the bulk attenuation and relaxation time emphasizes the importance of a number of parameters in the mineral physics database such as Fe-Mg diffusivity and kinetics in olivine-wadsleyite-ringwoodite-perovskite, Mg-Ca-Al-Si diffusivity and kinetics in pyroxene-garnet-Ca perovskite, some of which are still unknown.

58 GEOSCIENCES↗

Nanocomposite thermite powders with improved flowability prepared by mechanical milling

Nanocomposite thermite powders can be used in reactive parts and components. Manufacturing these components requires tailoring powder particle size distributions, particle shapes, and powder flowability. Specifically, an improved flowability is desired to use such powders as feedstock in additive manufacturing. Arrested reactive milling (ARM) offers a versatile and practical approach for preparing nanocomposite thermites with fully dense particles, which will retain their structures and mixedness between reactive components while being stored, handled, and processed. However, ARM products usually have broad particle size distributions, rock-like particle shapes, and poor flowability. Here in this study, ARM is modified to include a low-energy milling step to tune the shapes and flowability of the prepared powders. Experiments are performed with aluminum-rich Al·Fe 2 O 3 thermites. After the initial nanocomposite thermite is prepared in a planetary mill, it is additionally milled at a reduced rotation rate, replacing milling balls with smaller glass beads, and adding different liquid process control agents. Powders with modified particle shapes and size distributions are obtained, which have substantially improved flowability compared to the initial material. The reactivity of the initial and modified powders is evaluated using their ignition on a heated wire, by electro-static discharge, by thermal analysis (DSC) and in constant volume explosion tests. Finally, the reactivity of the modified powders is not diminished; instead, an improved reactivity is observed for selected samples.

36 MATERIALS SCIENCE↗

Mean annual precipitation predicts primary production resistance and resilience to extreme drought

Extreme drought is increasing in frequency and intensity in many regions globally, with uncertain consequences for the resistance and resilience of ecosystem functions, including primary production. Primary production resistance, the capacity to withstand change during extreme drought, and resilience, the degree to which production recovers, vary among and within ecosystem types, obscuring generalized patterns of ecological stability. Theory and many observations suggest forest production is more resistant but less resilient than grassland production to extreme drought; however, studies of production sensitivity to precipitation variability indicate that the processes controlling resistance and resilience may be influenced more by mean annual precipitation (MAP) than ecosystem type. Here, we conducted a global meta-analysis to investigate primary production resistance and resilience to extreme drought in 64 forests and grasslands across a broad MAP gradient. We found resistance to extreme drought was predicted by MAP; however, grasslands (positive) and forests (negative) exhibited opposing resilience relationships with MAP. Our findings indicate that common plant physiological mechanisms may determine grassland and forest resistance to extreme drought, whereas differences among plant residents in turnover time, plant architecture, and drought adaptive strategies likely underlie divergent resilience patterns. Finally, the low resistance and resilience of dry grasslands suggests that these ecosystems are the most vulnerable to extreme drought – a vulnerability that is expected to compound as extreme drought frequency increases in the future.

54 ENVIRONMENTAL SCIENCES↗

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

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

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Metal-embedded fiber optic sensor packaging and signal demodulation scheme towards high-frequency dynamic measurements in harsh environments

Optical fiber-based sensors are commonly used for industrial monitoring and control processes due to their small size, resistance to electromagnetic interference, and ability to perform a wide variety of high-precision measurements. However, implementing optical fiber sensors in harsh environments is challenging because they are small and fragile. Proper packaging of fiber-optic sensors could extend their use to harsh environments, including at high temperature and under high radiation. Furthermore, conventional fiber optic-based measurement systems often use computationally expensive signal processing algorithms which hinder their use in high-frequency dynamic applications. This work reports on the design of an optical fiber-pressure sensor system based on low-coherence interferometry that uses a metal-embedded optical fiber to provide a robust sensor package. Additionally, a novel phase demodulation scheme is proposed to extract length changes from the deformation of a thin diaphragm within a Fabry–Pérot cavity to measure pressure. The sensor has been tested to ±100 kPa, has a theoretical linear response over a 270 kPa dynamic range (24 dB maximum signal-to-noise ratio), and can resolve pressure transients up to 3910 kPa/s. Sampling at 100 kHz, the present sensor can resolve 2 kPa dynamic pressures at frequencies up to 2 kHz. Faster transients on the order of tens to hundreds of kHz can theoretically be resolved at the expense of decreasing the maximum resolvable amplitude. A methodology for designing a LCI pressure sensor for a given application is outlined based on the limitations imposed by the Nyquist criterion, the diaphragm's resonant frequency, the LCI optoelectronics, and the phase demodulation scheme. The sensor is the first to implement a low-coherence light source and a Fabry–Pérot interferometer designed to provide real-time high-fidelity pressure measurements using a metal-embedded optical fiber. The demonstrated sensor provides a platform for sensing in harsh conditions such as in nuclear and energy applications.

47 OTHER INSTRUMENTATION↗

Developmental exposure to 50 parts-per-billion arsenic influences histone modifications and associated epigenetic machinery in a region- and sex-specific manner in the adult mouse brain

Epidemiological studies report that arsenic exposure via drinking water adversely impacts cognitive development in children and, in adults, can lead to greater psychiatric disease susceptibility, among other conditions. While it is known that arsenic toxicity has a profound effect on the epigenetic landscape, very few studies have investigated its effects on chromatin architecture in the brain. We have previously demonstrated that exposure to a low level of arsenic (50 ppb) during all three trimesters of fetal/neonatal development induces deficits in adult hippocampal neurogenesis in the dentate gyrus (DG), depressive-like symptoms, and alterations in gene expression in the adult mouse brain. As epigenetic processes control these outcomes, here we assess the impact of our developmental arsenic exposure (DAE) paradigm on global histone posttranslational modifications and associated chromatin-modifying proteins in the dentate gyrus and frontal cortex (FC) of adult male and female mice. DAE influenced histone 3 K4 trimethylation with increased levels in the male DG and FC and decreased levels in the female DG (no change in female FC). The histone methyltransferase MLL exhibited a similar sex- and region-specific expression profile as H3K4me3 levels, while histone demethylase KDM5B expression trended in the opposite direction. DAE increased histone 3 K9 acetylation levels in the male DG along with histone acetyltransferase (HAT) expression of GCN5 and decreased H3K9ac levels in the male FC along with decreased HAT expression of GCN5 and PCAF. DAE decreased expression of histone deacetylase enzymes HDAC1 and HDAC2, which were concurrent with increased H3K9ac levels but only in the female DG. Levels of H3 and H3K9me3 were not influenced by DAE in either brain region of either sex. These findings suggest that exposure to a low, environmentally relevant level of arsenic during development leads to long-lasting changes in histone methylation and acetylation in the adult brain due to aberrant expression of epigenetic machinery based on region and sex. - Highlights: • Brain tissue from adult mice with developmental arsenic exposure (DAE) was used. • DAE impacted histone methylation and associated methyltransferases based on sex. • DAE differentially altered histone acetylation based on brain region. • DAE altered HATs in males and HDACs in females. • Epigenetic modifier expression correlated with the associated histone modification.

60 APPLIED LIFE SCIENCES↗

Kinetics of steel slag leaching: Batch tests and modeling

Reusing steel slag as an aggregate for road construction requires to characterize the leaching kinetics and metal releases. In this study, basic oxygen furnace (BOF) steel slag were subjected to batch leaching tests at liquid to solid ratios (L/S) of 10 and 100 over 30 days; the leachate chemistry being regularly sampled in time. A geochemical model of the steel slag is developed and validated from experimental data, particularly the evolution with leaching of mineralogical composition of the slag and trace element speciation. Kinetics is necessary for modeling the primary phase leaching, whereas a simple thermodynamic equilibrium approach can be used for secondary phase precipitation. The proposed model simulates the kinetically-controlled dissolution (hydrolysis) of primary phases, the precipitation of secondary phases (C-S-H, hydroxide and spinel), the pH and redox conditions, and the progressive release of major elements as well as the metals Cr and V. Modeling indicates that the dilution effect of the L/S ratio is often coupled to solubility-controlled processes, which are sensitive to both the pH and the redox potential. A sensitivity analysis of kinetic uncertainties on the modeling of element releases is performed.

36 MATERIALS SCIENCE↗

Community structure and function during periods of high performance and system upset in a full-scale mixed microalgal wastewater resource recovery facility

Microalgae have the potential to exceed current nutrient recovery limits from wastewater, enabling water resource recovery facilities (WRRFs) to achieve increasingly stringent effluent permits. The use of photobioreactors (PBRs) and the separation of hydraulic retention and solids residence time (HRT/SRT) further enables increased biomass in a reduced physical footprint while allowing operational parameters (e.g., SRT) to select for desired functional communities. However, as algal technology transitions to full-scale, there is a need to understand the effect of operational and environmental parameters on complex microbial dynamics among mixotrophic microalgae, bacterial groups, and pests (i.e., grazers and pathogens) and to implement robust process controls for stable long-term performance. Here, we examine a full-scale, intensive WRRF utilizing mixed microalgae for tertiary treatment in the US (EcoRecover, Clearas Water Recovery Inc.) during a nine-month monitoring campaign. We investigated the temporal variations in microbial community structure (18S and 16S rRNA genes), which revealed that stable system performance of the EcoRecover system was marked by a low-diversity microalgal community (D INVSIMPSON = 2.01) dominated by Scenedesmus sp. (MRA = 55 %-80 %) that achieved strict nutrient removal (effluent TP < 0.04 mg·L -1 ) and steady biomass concentration (TSS monthly avg . = 400–700 mg·L −1 ). Operational variables including pH, alkalinity, and influent ammonium (NH 4 + ), correlated positively (p < 0.05, method = Spearman) with algal community during stable performance. Further, the use of these parameters as operational controls along with N/P loading and SRT allowed for system recovery following upset events. Importantly, the presence or absence of bacterial nitrification did not directly impact algal system performance and overall nutrient recovery, but partial nitrification (potentially resulting from NO 2 − accumulation) inhibited algal growth and should be considered during long-term operation. The microalgal communities were also adversely affected by zooplankton grazers (ciliates, rotifers) and fungal parasites (Aphelidium), particularly during periods of upset when algal cultures were experiencing culture turnover or stress conditions (e.g., nitrogen limitation, elevated temperature). Altogether, the active management of system operation in order to maintain healthy algal cultures and high biomass productivity can result in significant periods (>4 months) of stable system performance that achieve robust nutrient recovery, even in winter months in northern latitudes (WI, USA).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reimagining pH Measurement: Utilizing Raman Spectroscopy for Enhanced Accuracy in Phosphoric Acid Systems

Measurement of pH is an integral component of chemical studies and process control; however traditional pH probes are difficult to utilize in harsh or complex chemical systems. Optical spectroscopy-based on-line monitoring offers a powerful and novel route for characterizing system parameters, such as pH, and is well adapted to deployment in harsh environments or chemically complex systems. Specifically, Raman spectroscopy combined with chemometric analysis can provide an improved method of on-line pH measurement. Therefore, multivariate curve resolution analysis of Raman spectra can be utilized to determine speciation as a function of pH and subsequent chemometric modeling can be used to correlate spectral response to pH. This was demonstrated with phosphoric acid, a chemical system known to challenge traditional pH probes. Raman spectra exhibit clear changes with pH due to changing speciation, and chemometric modeling can be successfully utilized to correlate those fingerprints to pH. Using this approach, true pH of the phosphoric acid system can be measured with greater accuracy than with traditional pH probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulated Aquifer Heterogeneity Leads to Enhanced Attenuation and Multiple Retention Processes of Zinc

Alluvial aquifers serve as one of the main water sources for domestic, agricultural, and industrial purposes globally. Groundwater quality, however, can be threatened by naturally occurring and anthropogenic metal contaminants. Differing hydrologic and biogeochemical conditions between predominantly coarse-grained aquifer sediments and embedded layers or lenses of fine-grained materials lead to variation in metal behavior. Here, we examine processes controlling Zn partitioning within a dual-pore domain-reconstructed alluvial aquifer. Natural coarse aquifer sediments from the Wind River-Little Wind River floodplain near Riverton, WY, were used in columns with or without fine-grained lenses to examine biogeochemical controls on Zn concentrations, retention mechanisms, and transport. Furthermore, following the introduction of Zn to the groundwater source, Zn preferentially accumulated in the fine-grained lenses, despite their small volumetric contributions. While the clay fraction dominated Zn retention in the sandy aquifer, the lenses supported additional reaction pathways of retention—the reducing conditions within the lenses resulted in ZnS precipitation, overriding the contribution of organic matter. Zinc concentration in the groundwater controlled the formation of Zn-clays and Zn-layered double hydroxides, whereas the extent of sulfide production controlled precipitation of ZnS. Our findings illustrate how both spatial and compositional heterogeneities govern the extent and mechanisms of Zn retention in intricate groundwater systems, with implications for plume behavior and groundwater quality.

54 ENVIRONMENTAL SCIENCES↗

Effect of Organic Coatings, Humidity and Aerosol Acidity on Multiphase Chemistry of Isoprene Epoxydiols

Multiphase chemistry of isomeric isoprene epoxydiols (IEPOX) has been shown to be the dominant source of isoprene-derived secondary organic aerosol (SOA). Recent studies have reported particles composed of ammonium bisulfate (ABS) mixed with model organics exhibit slower rates of IEPOX uptake. In the present study, we investigate the effect of atmospherically-relevant organic coatings of α-pinene (AP) SOA on the reactive uptake of trans-β-IEPOX onto ABS particles under different conditions and coating thicknesses. Single particle mass spectrometry was used to characterize in real-time particle size, shape, density, and quantitative composition before and after reaction with IEPOX. We find that IEPOX uptake by pure sulfate particles is a volume-controlled process, which results in particles with uniform concentration of IEPOX-derived SOA across a wide range of sizes. Aerosol acidity was shown to enhance IEPOX-derived SOA formation, consistent with recent studies. The presence of water has a weaker impact on IEPOX-derived SOA yield, but significantly enhanced formation of 2-methyltetrols, consistent with offline filter analysis. In contrast, IEPOX uptake by ABS particles coated by AP-derived SOA is strongly dependent on particle size and composition. IEPOX uptake occurred only when weight fraction of AP-derived SOA dropped below 50 %, effectively limiting IEPOX uptake to larger particles.

Environmental Molecular Sciences Laboratory↗

Hexavalent Chromium Generation within Naturally Structured Soils and Sediments

Chromium(VI) produced from the oxidation of indigenous Cr(III) minerals is increasingly being recognized as a threat to groundwater quality. A critical determinant of Cr(VI) generation within soils and sediments is the necessary interaction of two low-solubility phases$-$Cr(III) silicates or (hydr)oxides and Mn(III/IV) oxides—that lead to its production. Here in this paper, we investigate the potential for Cr(III) oxidation by Mn oxides within fixed solid matrices common to soils and sediments. Artificial aggregates were constructed from Cr(OH) 3 - and Cr 0.25 Fe 0.75 (OH) 3 -coated quartz grains and either mixed with synthetic birnessite or inoculated with the Mn(II)-oxidizing bacterium Leptothrix cholodnii. In aggregates simulating low organic carbon environments, we observe Cr(VI) concentrations within advecting solutes at levels more than twenty-times the California drinking water standard. Chromium(VI) production is highly dependent on Cr-mineral solubility; increasing Fe-substitution (x = 0 to x = 0.75) decreases the solubility of the solid and concomitantly decreases total Cr(VI) generation by 37%. In environments with high organic carbon, reducing conditions within aggregate cores (microbially) generate sufficient Fe(II) to suppress Cr(VI) efflux. Our results illustrate Cr(VI) generation from reaction with Mn oxides within structured media simulating soils and sediments and provide insight into how fluctuating hydrologic and redox conditions impact coupled processes controlling Cr and Mn cycling.

54 ENVIRONMENTAL SCIENCES↗

Kinetics of Cation and Oxyanion Adsorption and Desorption on Ferrihydrite: Roles of Ferrihydrite Binding Sites and a Unified Model

Understanding the kinetics of toxic ion reactions with ferrihydrite is crucial for predicting the dynamic behavior of contaminants in soil environments. In this study, the kinetics of As(V), Cr(VI), Cu, and Pb adsorption and desorption on ferrihydrite were investigated with a combination of laboratory macroscopic experiments, microscopic investigation and mechanistic modeling. The rates of As(V), Cr(VI), Cu, and Pb adsorption and desorption on ferrihydrite, as systematically studied using a stirred-flow method, was highly dependent on the reaction pH and metal concentrations and varied significantly among four metals. Spherical aberration-corrected scanning transmission electron microscopy (Cs-STEM) showed, at sub-nano scales, all four metals were distributed within the ferrihydrite particle aggregates homogeneously after adsorption reactions, with no evidence of surface diffusion-controlled processes. Based on experimental results, we developed a unifying kinetics model for both cation and oxyanion adsorption/desorption on ferrihydrite based on the mechanistic-based equilibrium model CD-MUSIC. Overall, the model described the kinetic results well, and we quantitatively demonstrated how the equilibrium properties of the cation and oxyanion binding to various ferrihydrite sites affected the adsorption and desorption rates. Our results provided a unifying quantitative modeling method for the kinetics of both cation and oxyanion adsorption/desorption on iron minerals.

Environmental Molecular Sciences Laboratory↗

Calcium-Facilitated Aggregation and Precipitation of the Uranyl Peroxide Nanocluster U 60 in the Presence of Na-Montmorillonite

The unique and diverse features of uranyl peroxide nanoclusters may contribute to the enhanced mobility of uranium in the environment. This study examines the sorption of the uranyl peroxide nanocluster [UO 2 (O 2 )(OH)] 60 60- (U 60 ) to Na-montmorillonite (SWy-2), plagioclase (anorthite), and quartz (SiO 2 ) as a function of time, U 60 concentration, and mineral concentration. SWy-2 was studied in both its untreated form as well as after two different pretreatments, denoted as partially treated SWy-2 and fully treated SWy-2. U 60 was removed (~99%) from solution in the presence of untreated and partially treated SWy-2. However, U 60 was not removed from suspensions containing anorthite, quartz, or fully treated SWy-2, even after several months. The removal of U 60 from suspensions containing untreated SWy-2 is promoted in part by the exchange of Li + counter-ions, normally weakly associated with U 60 in solution, for Ca 2+ ions naturally present in the clay. In solution, Ca 2+ ions induce the aggregation of nanoclusters, which precipitate on the surface of SWy-2. Ca-rich U 60 aggregates associated with SWy-2 were identified and characterized by scanning electron microscopy with energy dispersive spectroscopy, Raman spectroscopy, and X-ray photoelectron spectroscopy. This research enhances our understanding of the molecular-scale processes controlling U 60 behavior at the mineral–water interface.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Trace Element Mass Flow Rates from U.S. Coal Fired Power Plants

Trace elements (TEs) exit coal-fired power plants (CFPPs) via solid, liquid, and gaseous waste streams. Estimating the TE concentrations of these waste streams is essential to selecting pollution controls and estimating emission reduction benefits. Our report presents a generalizable mass balance model for estimating TE mass flow rates in CFPP waste streams and evaluates model accuracy for the U.S. coal fleet given current data constraints. We stochastically estimate, using a bootstrapping approach, the 2015 plant-level mass flow rates of Hg, Se, As, and Cl to solid, liquid, and gas phase waste streams by combining publicly available data for combusted coal TE concentrations with estimates of TE partitioning within installed air pollution control processes. When compared with measured and reported data on TE mass flow rates, this model generally overestimates masses by 30–50%, with larger errors for Hg. The partitioning estimates are consistent for Se, As, and Cl removal from flue gas, but tend to underestimate Hg removal. While our model is suitable for first-order estimates of TE mass flows, future work to enhance model performance should focus on collecting and using new data on TE concentrations in the coal blend, where data quality is the weakest.

01 COAL, LIGNITE, AND PEAT↗

Development and Validation of a Dynamic Response Model for a Cold Flow Circulating Fluidized Bed

A dynamic response model has been created of the Cold Flow Circulating Fluidized Bed at NETL in Morgantown, WV. The model is considered dynamic, as opposed to steady state, because it is used to predict the transient responses of gas and solids flows throughout the unit as process conditions vary. The model was built in ProTrax with Advanced Continuous Simulation Language (ACSL) custom code. Dynamic experiments, varying both gas and solids flows, in the CFB unit have been used to validate the accuracy of the model. The advantage of this model is its utility in predicting gas hold up and solids flow throughout the unit during transitions of process conditions. Accurate predictions of these phenomena are important in process control and optimization. The dynamic experiments used in validating the model include modulations of the system pressure, move air (solids flow), and riser gas velocity with step and sinusoidal changes.

01 COAL, LIGNITE, AND PEAT↗

Molecular Dynamics Simulations of the Adsorption of Phthalate Esters on Smectite Clay Surfaces

The partitioning of organic contaminants between water and solid surfaces is a key process controlling their fate and transport in natural environments. A novel methodology was developed to predict the adsorption of organic contaminants by smectite clay minerals (high specific surface area adsorbents abundant in natural soils) using molecular dynamics (MD) simulations. The methodology models a stack of flexible Ca- montmorillonite lamellae in direct contact with a bulk aqueous reservoir and uses the metadynamics technique to facilitate the exploration of the free energy landscape. The methodology was tested and validated in the case of six phthalate esters, widely used chemical plasticizers with endocrine disrupting properties. Simulation predictions reveal strong phthalate adsorption, especially for the larger and more hydrophobic phthalates. Predicted partition coefficients (K d values) are consistent with collected batch experimental data. Adsorption was observed on both the exterior basal surfaces and within the interlayer nanopore, with phthalate molecules predominately adopting a flat orientation on the clay surface. Intercalation was also detected in complementary X-ray diffraction (XRD) experiments. A strong inverse relationship between extent of adsorption and clay surface charge density was observed, as phthalate molecules preferentially occupied the more hydrophobic uncharged patches on each surface. Detailed analysis of the free energy of adsorption revealed that phthalate affinity for the clay surface results from a small favorable van der Waals contribution and a large favorable entropic contribution. Overall, this research demonstrates the substantial affinity of smectite clays for phthalate esters and establishes a computational methodology capable of predicting the water-clay partition coefficients of organic contaminants, a key parameter in environmental fate and transport models.

36 MATERIALS SCIENCE↗