Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CLAYS”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Encapsulation method for preparation of pellets with high attrition resistance

This disclosure provides a method to produce highly attrition resistant pellets by encapsulating reactive components in a vitrified clay outer layer. The reactive component mixture is present relative to the clay substrate in a weight ratio of part per 60-100 part to about 60 parts of the clay substrate. The reactive components are agglomerated first, and clay substrate is added to form the outer layer of the pellet. The pellets are calcined at temperatures above 1200 C to form a vitrified clay semi porous outer layer providing high strength to the pellet while facilitating the gas transfer for the reaction with the encapsulated reactive components. Pellets containing CuO—Fe 2 O 3 -alumina oxygen carrier for chemical looping combustion of fuel demonstrated high attrition resistance and high reactivity with methane.

Siriwardane, Ranjani↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
DOE OSTI2024

Induced Polarization of Clayey Rocks and Soils: Non‐Linear Complex Conductivity Models

Abstract The past decades have witnessed the increased applications of induced polarization (IP) method in the critical zone studies with ubiquitous clay minerals. Although IP outperforms traditional electrical and electromagnetic methods through its unique ability to measure quadrature conductivity, the nonlinearity that quadrature conductivity behaves with salinities and frequencies greatly tortures IP practitioners, as (a) salinity‐dependency makes the quadrature conductivity a varyingly unstable parameter to quantitatively estimate hydraulic properties and clay content; (b) frequency‐dependent Cole‐Cole and Debye/Warburg decomposition models, although mathematically sound, physically mingle the properties of pore water and clay minerals and are empirical in nature. From basic principles, we demonstrate that quadrature conductivity remains a hybrid property involving both clay and water, and develop relevant models to distinguish them. Our models are validated by theories, experiments, simulations, and comparisons, all of which proclaim considerable advantages over previous models and offer the prospect of quantitative applications.

58 GEOSCIENCES↗
DOE OSTI2022

Synthetic antibacterial minerals: harnessing a natural geochemical reaction to combat antibiotic resistance

The overuse of antibiotics in clinical and livestock settings is accelerating the selection of multidrug resistant bacterial pathogens. Antibiotic resistant bacteria result in increased mortality and financial strain on the health care and livestock industry. The development of new antibiotics has stalled, and novel strategies are needed as we enter the age of antibiotic resistance. Certain naturally occurring clays have been shown to have antimicrobial properties and kill antibiotic resistant bacteria. Harnessing the activity of compounds within these clays that harbor antibiotic properties offers new therapeutic opportunities for fighting the potentially devastating effects of the post antibiotic era. However, natural samples are highly heterogenous and exhibit variable antibacterial effectiveness, therefore synthesizing minerals of high purity with reproducible antibacterial activity is needed. Here we describe for the first time synthetic smectite clay minerals and Fe-sulfide microspheres that reproduce the geochemical antibacterial properties observed in natural occurring clays. We show that these mineral formulations are effective at killing the ESKAPE pathogens ( Enterococcus sp., Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter sp., Pseudomonas aeruginosa and Enterobacter sp.) by maintaining Fe 2+ solubility and reactive oxygen species (ROS) production while buffering solution pH, unlike the application of metals alone. Our results represent the first step in utilizing a geochemical process to treat antibiotic resistant topical or gastrointestinal infections in the age of antibiotic resistance.

58 GEOSCIENCES↗
DOE OSTI2022

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

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

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Soil properties explain tree growth and mortality, but not biomass, across phosphorus-depleted tropical forests

We observed strong positive relationships between soil properties and forest dynamics of growth and mortality across twelve primary lowland tropical forests in a phosphorus-poor region of the Guiana Shield. Average tree growth (diameter at breast height) increased from 0.81 to 2.1 mm yr -1 along a soil texture gradient from 0 to 67% clay, and increasing metal-oxide content. Soil organic carbon stocks in the top 30cm ranged from 30 to 118 tons C ha -1 , phosphorus content ranged from 7 to 600 mg kg -1 soil, and the relative abundance of arbuscular mycorrhizal fungi ranged from 0 to 50%, all positively correlating with soil clay, and iron and aluminum oxide and hydroxide content. In contrast, already low extractable phosphorus (Bray P) content decreased from 4.4 to <0.0 mg kg -1 in soil with increasing clay content. A greater prevalence of arbuscular mycorrhizal fungi in more clayey forests that had higher tree growth and mortality, but not biomass, indicates that despite the greater investment in nutrient uptake required, soils with higher clay content may actually serve to sustain high tree growth in tropical forests by avoiding phosphorus losses from the ecosystem. Our study demonstrates how variation in soil properties that retain carbon and nutrients can help to explain variation in tropical forest growth and mortality, but not biomass, by requiring niche specialization and contributing to biogeochemical diversification across this region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Soil resistance and recovery during neotropical forest succession

The recovery of soil conditions is crucial for successful ecosystem restoration and, hence, for achieving the goals of the UN Decade on Ecosystem Restoration. Here, we assess how soils resist forest conversion and agricultural land use, and how soils recover during subsequent tropical forest succession on abandoned agricultural fields. Our overarching question is how soil resistance and recovery depend on local conditions such as climate, soil type and land-use history. For 300 plots in 21 sites across the Neotropics, we used a chronosequence approach in which we sampled soils from two depths in old-growth forests, agricultural fields (i.e. crop fields and pastures), and secondary forests that differ in age (1–95 years) since abandonment. We measured six soil properties using a standardized sampling design and laboratory analyses. Soil resistance strongly depended on local conditions. Croplands and sites on high-activity clay (i.e. high fertility) show strong increases in bulk density and decreases in pH, carbon (C) and nitrogen (N) during deforestation and subsequent agricultural use. Resistance is lower in such sites probably because of a sharp decline in fine root biomass in croplands in the upper soil layers, and a decline in litter input from formerly productive old-growth forest (on high-activity clays). Soil recovery also strongly depended on local conditions. During forest succession, high-activity clays and croplands decreased most strongly in bulk density and increased in C and N, possibly because of strongly compacted soils with low C and N after cropland abandonment, and because of rapid vegetation recovery in high-activity clays leading to greater fine root growth and litter input. Furthermore, sites at low precipitation decreased in pH, whereas sites at high precipitation increased in N and decreased in C : N ratio. Extractable phosphorus (P) did not recover during succession, suggesting increased P limitation as forests age. These results indicate that no single solution exists for effective soil restoration and that local site conditions should determine the restoration strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Controlling Sustainability of Hydraulic Fracture Permeability in Ductile Shales

Hydraulic fracturing technology, along with horizontal drilling and associated completion technologies, dramatically increased US oil and gas production, by increasing the permeability of tight shales and allowing for cost-effective extraction of trapped hydrocarbons. For efficient and economical production, the targeted shale has to be fractured easily (i.e., good ‘frackability’), but also the permeability of induced fractures have to be sustained during production (good ‘sustainability’). For these reasons, currently, ductile shale with high clay content (>~40%) is difficult to exploit as a hydrocarbon resource, although hydrocarbons can still be found in it. Good frackability, high-TOC, low-clay-content resource shales—the ‘low-hanging fruit’—will be depleted, and we need to develop tools that allow tapping into currently underutilized, clay-rich ductile shales. This research project aims (1) to understand the behavior of fractures in clay-rich, ductile (and sometimes swelling) shales and (2) to begin to develop technologies for efficient and economical production from such shales. With these objectives, we examine time-dependent, coupled mechanical-hydrological behavior of open and proppant-filled fractures within different types of shales. In the preceding project, the focus of the research was to understand the behavior of the fractures and to predict their permeability changes due to fracture closure from shale matrix deformation and proppant embedment. In the current project, we investigate the possibility of chemically manipulating the sustainability of hydraulic fractures in ductile shales—specifically by altering the proppant-embedding behavior—via chemical means, at different stages of oil and gas production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Coal-based Bricks & Blocks (CBBs): Process Development to Prototype Fabrication Coupled with Techno-Economic Analysis and Market Survey

The decline of coal use for energy production provides an abundance of local feedstock for new innovative uses and value-added products. Expanding the U.S. coal-value chain to manufacture high-value carbon products can strengthen the nation’s energy and mineral security, enhance the U.S. national defense security, increase the United States’ economic prosperity, while achieving U.S. environmental objectives. The primary project goal is a relatively light-weight composite product with superior (or comparable) compressive strength. Both virgin and, where available, post-consumer recycled thermoplastic versions are tested for each thermoplastic species. A useful attribute of thermoplastics as binders is that they can be heated to their melting point, cooled, and reheated again without significant degradation. A key advantage of thermoset CBBs is that they require only mixing and molding. CBB advantages include low cost, availability, binding ability and processability. Coal-based bricks and blocks (CBBs) weigh about 50% less than clay bricks and can be manufactured with an interlocking design to promote ease of use for the novice builder. CBB formulation is evaluated according to a design-of-experiments (DoE) approach. DoE variables are a) relative weight fractions of binder, b) relative proportions of large versus small (milled) anthracite size fractions, and c) additive percentage. Fabrication methods include hot-press molding and extrusion, the later being the most commercially viable. CBBs are tested for compressive strength, modulus of rupture (by flexure test) and water absorption per ASTM C67, with density determined by the Archimedes drainage method. Fractured interfaces are examined by SEM (Scanning Electron Microscopy) to resolve fracture dynamics and interior microstructure uniformity. Differential scanning calorimetry (DSC) is used to compare plastic transition temperatures i.e., glass and melting temperatures to contrast virgin with post-consumer recycled thermoplastics and optimize their usage. These results are used in the DoE analysis to identify the binder and relative weight percentages for optimum strength, density, and porosity. Overall, CBBs possess strength comparable to clay-based bricks but are non-permeable and hydrophobic, and hence resistant to degradation by freeze-fracturing, corrosion, and efflorescence. The strongest composites have been made with the following thermoplastic binders (in order of strength): thermoset, high-density polyethylene crosslink resin, high-density polyethylene, nylon 6/6, and polypropylene. Results from a techno-economic analysis TEA show economy of scale for CBBs by modularization and reveal the binder as the cost driver for material costs. Ideally, the incorporation of post-consumer recycled thermoplastic will decrease material acquisition costs and increase product sustainability. Notably, CBBs do not require the high temperature calcination needed to produce cement, nor do they require firing in the 1600-2400 °F range for three days using natural gas, as do clay brick equivalents. Instead, CBBs are heated to a modest <600 °F according to the melt flow index of the thermoplastic binder. Existing anthracite mines can be expanded to produce CBBs to reduce aggregate transportation costs and emissions that exist for clay bricks. TEA reflects this reduced energy cost while a comparative CO2 emission analysis quantifies the reduced environmental footprint. The market survey identifies several commercialization opportunities, dependent upon the brick classification.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

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

Impacts of Crystalline Host Rock on Repository Barrier Materials at 250 °C: Hydrothermal Co-Alteration of Wyoming Bentonite and Steel in the Presence of Grimsel Granodiorite

Direct disposal of dual-purpose canisters (DPC) has been proposed to streamline the disposal of spent nuclear fuel. However, there are scenarios where direct disposal of DPCs may result in temperatures in excess of the specified upper temperature limits for some engineered barrier system (EBS) materials, which may cause alteration within EBS materials dependent on local conditions such as host rock composition, chemistry of the saturating groundwaters, and interactions between barrier materials themselves. Here we report the results of hydrothermal experiments reacting EBS materials—bentonite buffer and steel—with an analogue crystalline host rock and groundwater at 250 °C. Experiment series explored the effect of reaction time on the final products and the effects of the mineral and fluid reactants on different steel types. Post-mortem X-ray diffraction, electron microprobe, and scanning electron microscopy analyses showed characteristic alteration of both bentonite and steel, including the formation of secondary zeolite and calcium silicate hydrate minerals within the bentonite matrix and the formation of iron-bearing clays and metal oxides at the steel surfaces. Swelling clays in the bentonite matrix were not quantitatively altered to non-swelling clay species by the hydrothermal conditions. The combined results of the solution chemistry over time and post-mortem mineralogy suggest that EBS alteration is more sensitive to initial groundwater chemistry than the presence of host rock, where limited potassium concentration in the solution prohibits conversion of the smectite minerals in the bentonite matrix to non-swelling clay species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Layer size polydispersity in hydrated montmorillonite creates multiscale porosity networks

The aluminosilicate layers of the swelling clay mineral montmorillonite, and the saturated pores they delineate, control the mechanical properties and the transport of solutes in many natural and engineered environments. However, the structural basis of montmorillonite porosity remains poorly characterized due to the difficulty in visualizing hydrated samples in their native state. Here in this paper, we used cryogenic transmission electron microscopy (cryo-TEM) and cryo electron tomography (cryo-ET) to show that stacking defects within minimally altered, fully hydrated montmorillonite particles define multiscale porosity networks. Variations in layer lateral dimensions over tens to thousands of nanometers cause a range of topological and dynamic defects that generate pervasive curvature and introduce previously uncharacterized solute transport pathways. Observations of long-range rotational order between neighboring layers indicate that the layer-layer interactions that govern clay swelling involve three dimensional orienting forces that operate across nanoscale pores. These direct observations of the hierarchical structure of hydrated montmorillonite pore networks with nanoscale resolution reveal potentially general aspects of colloidal interactions in fluid-saturated clay minerals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Illitization in the Mt. Simon Sandstone, Illinois Basin, USA: Implications for carbon dioxide storage

Clay mineral cementation is one of the most important controls on sandstone reservoir properties. Here, in this paper, the diagenetic history of the Mt. Simon reservoir complex is studied to reveal the origin, timing, and controls of illitization. Samples of sandstone and shale from Mt. Simon reservoir complex were acquired from the Illinois Basin–Decatur Project (IBDP), a CO 2 storage demonstration project in the central Illinois Basin. Petrographic, SEM, TEM, XRD analyses and K/Ar age dating were completed to identify the major detrital and diagenetic components of the samples and reveal illite to be the major clay component in all samples. Illitic clay coatings in the lower Mt. Simon reservoir are identified as a major control on reservoir properties by inhibiting major precipitation of authigenic quartz during illitization, resulting in highly permeable sandstone, essential for CO 2 storage. The coating box-work morphology and mineralogy are indicative of a detrital smectite origin with subsequent illite growth associated with feldspar dissolution and kaolinite alteration. The mineralogy of bulk material and separate grain size fractions (2–0.6 μm; 0.6–0.2 μm; < 0.2 μm) with illite polytypes 2M 1 , 1M, and 1M d were quantified and age dated via 40 K- 40 Ar methods. The shales or reservoir seals contain the highest proportions of detrital illite with illite in the lower Mt. Simon Sandstone reservoir identified as solely diagenetic. Two major events of illitization are identified throughout the Mt. Simon with more porous reservoir rock exhibiting the older event from approximately 360 to 315 Ma and tighter sandstone units with reservoir properties too low to be considered reservoir exhibiting illite dates from 250 to 220 Ma. This partitioning of illitization is attributed to varied mineralogy controlled by depositional changes and evolution of the greater basin.

58 GEOSCIENCES↗
DOE OSTI2021

Molecular Determination of Organic Adsorption Sites on Smectite during Fe Redox Processes Using ToF-SIMS Analysis

Turnover of soil organic carbon (SOC) is strongly affected by a balance between mineral protection and microbial degradation. However, the mechanisms controlling the heterogeneous and preferential adsorption of different types of SOC remain elusive, especially for humic substrates (HSs) and microbial carbon (MC), two hotly-debated forms of SOC prevalent in soils. In this work, the heterogeneous adsorption of HSs and MC on a clay mineral (nontronite NAu-2) during microbial-mediated Fe redox cycling was determined with the use of time-of-flight secondary ion mass spectrometry (ToF-SIMS). The results revealed that HSs pre-adsorbed on NAu-2 would partially inhibit the microbial Fe(III) reduction and structural modification of NAu-2, thus retarding the subsequent adsorption of MC. In contrast, NAu-2 without pre-coated HSs was adsorbed by a significant amount of MC from microbial polysaccharides as a result of Fe(III) reduction. This was attributed to the deposition of a thin Al-rich layer on clay surface, which provided active sites for MC adsorption. This study provides direct and detailed molecular evidence for the first time to explain the preferential adsorption of MC over HSs on the surface of clay minerals in the iron redox processes, which could be critical for the preservation of MC in soil. The results also indicate that ToF-SIMS is a unique tool for understanding complex organic-mineral-microbe interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Dynamic Evaluation of the Upper Tyler Formation and Well Stimulation Fluid Interactions Using Micro-CT Imaging

Abundant concentrations of swelling clays in the oil-bearing upper Tyler Formation inhibit unconventional well stimulation techniques and associated long-term oil and gas production success. Laboratory evaluation of the geochemical interactions between the formation material and various stimulation fluids may help identify innovative approaches that provide a solution to successful well stimulation and subsequent oil production. The objective of this research was to understand the complexities of well stimulation fluid and clay mineral interactions within the Tyler Formation and identify potential fracturing fluid formulations that mitigate swelling properties of the clays in the reservoir to enhance stimulation success and promote long-term oil and gas production. Collaboration with the National Energy Technology Laboratory (NETL), utilizing their Tescan DynaTOM micro-CT analysis instrumentation, provided an innovative approach to understand real-time, dynamic interactions of the formation material and various potential stimulation fluids. Results are anticipated to identify key mechanisms occurring at the micro-scale level and provide insight into modified stimulation techniques uniquely suited for successful production applications.

enhanced oil recovery↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗