Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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 1,189 records · Page 66

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗

Effect of decades of corrosion on the microstructure of altered glasses and their radiation stability

Understanding the microstructural evolution of glasses during their interaction with water and radiation is of fundamental importance in addressing the corrosion of nuclear waste forms under geological disposal conditions. Here we report the results of more than 21 years of corrosion of two borosilicate glasses showing the formation of mesoporous C–S–H gels in Ca-bearing glasses and a mainly microporous microstructure in Al-bearing glasses. These porous corroded glasses were then irradiated with heavy ions to simulate the effects of recoil nucleus damage and monitored in real time using transmission electron microscopy with in situ ion irradiation. The ballistic collisions remarkably healed the porous corroded glasses to a pore-free homogeneous microstructure. Besides providing new insights and predictions about how doped glasses and actual waste forms may evolve under corrosion and irradiation, the results highlight the non-universal nature of the existing corrosion models and the important role that the glass composition and radiation damage play in the evolution of the microstructure during corrosion.

36 MATERIALS SCIENCE↗

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Leveraging the signature of heterotrophic respiration on atmospheric CO 2 for model benchmarking

Spatial and temporal variations in atmospheric carbon dioxide (CO 2 ) reflect large-scale net carbon exchange between the atmosphere and terrestrial ecosystems. Soil heterotrophic respiration (HR) is one of the component fluxes that drive this net exchange, but, given observational limitations, it is difficult to quantify this flux or to evaluate global-scale model simulations thereof. Here, we show that atmospheric CO 2 can provide a useful constraint on large-scale patterns of soil heterotrophic respiration. We analyze three soil model configurations (CASA-CNP, MIMICS, and CORPSE) that simulate HR fluxes within a biogeochemical test bed that provides each model with identical net primary productivity (NPP) and climate forcings. We subsequently quantify the effects of variation in simulated terrestrial carbon fluxes (NPP and HR from the three soil test-bed models) on atmospheric CO 2 distributions using a three-dimensional atmospheric tracer transport model. Our results show that atmospheric CO 2 observations can be used to identify deficiencies in model simulations of the seasonal cycle and interannual variability in HR relative to NPP. In particular, the two models that explicitly simulated microbial processes (MIMICS and CORPSE) were more variable than observations at interannual timescales and showed a stronger-than-observed temperature sensitivity. Our results prompt future research directions to use atmospheric CO 2 , in combination with additional constraints on terrestrial productivity or soil carbon stocks, for evaluating HR fluxes.

54 ENVIRONMENTAL SCIENCES↗

Modeling biological nitrogen fixation in global natural terrestrial ecosystems

Biological nitrogen fixation plays an important role in the global nitrogen cycle. However, the fixation rate has been usually measured or estimated at a particular observational site. To quantify the fixation amount at the global scale, process-based models are needed. This study develops a biological nitrogen fixation model to quantitatively estimate the nitrogen fixation rate by plants in a natural environment. The revised nitrogen module better simulates the nitrogen cycle in comparison with our previous model that has not considered the fixation effects. The new model estimates that tropical forests have the highest fixation rate among all ecosystem types, which decreases from the Equator to the polar region. The estimated nitrogen fixation in global terrestrial ecosystems is 61.5 Tg N yr-1 with a range of 19.8–107.9 Tg N yr-1 in the 1990s. Our estimates are relatively low compared to some early estimates using empirical approaches but comparable to more recent estimates that involve more detailed processes in their modeling. Furthermore, the contribution of nitrogen made by biological nitrogen fixation depends on ecosystem type and climatic conditions. This study highlights that there are relatively large effects of biological nitrogen fixation on ecosystem nitrogen cycling. and the large uncertainty of the estimation calls for more comprehensive understanding of biological nitrogen fixation. More direct observational data for different ecosystems are in need to improve future quantification of fixation and its impacts.

54 ENVIRONMENTAL SCIENCES↗

Strategies for conservative and non-conservative monotone remapping on the sphere

Abstract. Monotonicity is an important property of remapping operators for coupled weather and climate models. However, it is often challenging to design highly accurate operators that avoid the generation of new extrema or keep a remapped field between physically prescribed bounds. To that end, this paper explores several traditional and novel approaches for both conservative and non-conservative monotone remapping on the sphere. The accuracy and effectiveness of these algorithms are evaluated in the context of several different real and idealized fields and meshes.

54 ENVIRONMENTAL SCIENCES↗

Hydromechanical impact of basement rock on injection-induced seismicity in Illinois Basin

Abstract The common explanation of observed injection-induced microseismicity is based on the measured stress state at the injection interval and the assumption that it remains the same in the vicinity. We argue here that representing the stress state in different geologic formations over the injection site with the single Mohr’s circle is insufficient due to local stratigraphic features and contrast in compressibilities of the involved formations. The role of hydromechanical coupling in the microseismic response is also crucial for the proper assessment of the problem. Thoroughly monitored Illinois Basin Decatur Project revealed the majority of CO 2 injection-associated microseismic events being originated in the crystalline basement. Even though basement faults can serve as the conduits for fluid flow—the predicted pressure increase seems to be insufficient to trigger seismicity. To address this issue, accurate laboratory measurements of rock properties from the involved formations are conducted. The pre-injection stress state and its evolution are evaluated with the hydromechanically coupled numerical model. It appears that the presence of an offset in a stiff competent layer affects the stress state in its vicinity. Therefore, both the pre-injection stress state and its evolution during the fluid injection should be addressed during the induced seismicity assessment.

54 ENVIRONMENTAL SCIENCES↗

Basin & Range Investigation for Developing Geothermal Energy

Hidden geothermal systems represent a potentially prolific energy resource that could support critical U.S. public and government energy priorities. Basin and Range Investigations for Developing Geothermal Energy (BRIDGE) addressed some the challenges associated with hidden system exploration by prioritizing cost-effective exploration early on through strategic workflow and informed decision-making that mitigates early risk and shifts resources to later exploration stages (e.g., drilling). Sandia National Laboratories partnered with U.S. Navy Geothermal Office, Geologic Geothermal Group, and independent consultants, with additional collaboration with U.S. Geological Survey and private industry. The primary tool of the BRIDGE project was to deploy a regional-scale airborne electromagnetic method to investigate the shallow resistivity structure in areas with high prospectivity. This was followed up at several prospects by a multidisciplinary exploration approach, including additional geologic, geophysical and geochemical studies. A central tenet to the BRIDGE methodology is that zones of low resistivity frequently occur over geothermal systems in the Basin and Range, and when paired with other data constraints, imaging these zones can enable discovery of these systems. In addition to exploring greenfield areas (i.e., Grover Point), the BRIDGE project also flew HTEM resistivity surveys over known geothermal systems including those with established power plants (Don A. Campbell and Salt Wells) and prospects that are known to the literature but remain undeveloped, at least in part, due to a lack of understanding on the location of their producible reservoirs. BRIDGE produced a comprehensive set of data from prospects identified in the Nevada Play Fairway Analysis along with conceptual models for top ranking prospects, wherein all of the observations are used to inform an interpreted model of the system. These models present a range of possible system parameters such as temperature and size, and they are further informed by system analogues in the Basin and Range province and elsewhere. The results of this work leave space for further exploration that may now occur at prospects ‘down the list’ rather than distribution exploration resources evenly across all prospects.

15 GEOTHERMAL ENERGY↗

State Charlton & MRCSP 1 Characterization Well Report

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the efforts to drill and complete a well in the St. Peter sandstone for characterization and evaluation of storage potential in norther Michigan. However, as presented in this report, the well was not drilled and completed up to the intended target intervals due to drilling related difficulties. While the drilling issues were being resolved, the new Class VI injection well regulations were finalized by the US EPA. MRCSP efforts focused on evaluating CO 2 storage and containment using site characterization, modeling, and monitoring in conjunction with the ongoing CO 2 -EOR operations in the same area. This report describes work associated with the initial test well, the State Chester and MRCSP Well No. 1, drilled to a depth of about 7,968 feet TVD, and subsequently plugged and designated as temporarily abandoned. A high-level overview of data collected and analyzed is provided for shallower units.

01 COAL, LIGNITE, AND PEAT↗

Techno-economic feasibility of geothermal energy production using inactive oil and gas wells for district heating and cooling systems in Tuttle, Oklahoma

Oil and gas wells have been drilled in the United States and often encounter hot co-produced water possibly suitable for geothermal direct-use applications. Here this study evaluated geothermal resources, heating and cooling demand, and techno-economic potential of four oil and gas wells repurposed for geothermal energy production to serve two public schools and 250 nearby houses in Tuttle, Oklahoma, USA. Subsurface geology in the study area primarily consists of Permian to Mississippian sandstone and limestone formations approximately from 2 km to 3.5 km depth. With a geothermal gradient of 25.3 °C/km, the temperature of geothermal energy production was estimated as 68.2 °C at 2.1 km. The potential of calcite scaling within boreholes and pipes was minimal. Using the characterized reservoir properties, techno-economic analyses were conducted for evaluating levelized costs of geothermal energy production for direct-use heating and cooling and district heating systems with five different production scenarios. Heating and cooling demand in two schools and 250 houses and additional 0.3 MW and 0.6 MW base heating loads for potential geothermal direct-use applications in the study area were also incorporated in the modeling. Results indicated that the levelized cost of heating of geothermal district heating systems utilizing existing wells were significantly lower than those of conventional geothermal energy systems that include well drillings. The geothermal energy production with doublet and quartet configurations was generally sufficient to supply total annual heating demand in the district, while a peaking boiler was used to supply peak loads. Doublet configurations showed higher production temperature with a lower thermal drawdown than the quartet configurations. The doublet system at 2.1 km production depth with 9.3 kg/s flow rate represented the optimal levelized cost of $71/MWh and 91.2% utilization factor. These results imply that the Tuttle geothermal district heating system is techno-economically feasible and comparable to geothermal district heating systems in the United States and the regional natural gas rate for residentials by recycling inactive oil and gas wells with no drillings.

15 GEOTHERMAL ENERGY↗

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗

Thermoelastic Properties of Fe 3+ ‐Rich Jeffbenite and Application to Superdeep Diamond Barometry

Abstract Jeffbenite (Mg 3 Al 2 Si 3 O 12 ) is a tetragonal phase found in so far only in superdeep diamonds, and its thermoelastic parameters are a prerequisite for determining entrapment pressures as it is regarded as a potential indicator for superdeep diamonds. In this study, the thermoelastic properties of synthetic Fe 3+ ‐jeffbenite were measured up to 33.7 GPa and 750 K. High‐temperature static compression data were fitted, giving ( ∂K T 0 / ∂ T ) P = −0.0107 (4) GPa/K and α T = 3.50 (3) × 10 −5 K −1 . The thermoelastic properties and phase stability are applied to modeling isomekes, or P‐T paths intersecting possible conditions of entrapment in diamond. We calculate that under ideal exhumation, jeffbenite entrapped at mantle transition zone conditions will exhibit a high remnant pressure at 300 K ( P inc ) of ∼5.0 GPa. Elastic geobarometry on future finds of jeffbenite inclusions can use the new equation of state to estimate entrapment pressures for this phase with still highly uncertain stability field in the mantle.

Geology↗

NETL’s Perspective on Storage Efficiency and CO2-SCREEN

Carbon capture and storage (CCS) is a process that captures carbon dioxide (CO2) by separating it from anthropogenic emissions sources before atmospheric release and storing that CO2 in deep geologic reservoirs. CCS is a powerful method for reducing anthropogenic CO2 which can ultimately diminish the effects of climate change. Prospective CO2 storage resource is the amount of carbon dioxide that can be stored in a given geologic formation typically given as a mass (e.g., metric tons). Obtaining accurate estimates of CO2 storage resources is necessary for governments and industries to make energy-related policy decisions. Researchers at the National Energy Technology Laboratory (NETL) under the Department of Energy (DOE) developed a methods and a tool [CO2-SCREEN (Storage prospeCtive Resource Estimation Excel aNalysis) to estimate prospective carbon storage resources for saline formations, unconventional shale formations, and residual oil zones. The methods and tool provide CO2 storage and efficiency outputs in the form of probability estimates (i.e. P10 and P90) as well as partitioning storage and efficiency estimates based on storage mechanism (total, free phase, sorbed phase, and dissolution phase). This presentation will focus on how storage efficiency is calculated based on numerical modeling efforts, how it’s applied in the storage methods and tool, and then highlighting needs for future development.

Hanson, Angela Goodman↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

Conceptual Model Testing Related to SDU 6 Drainwell Observations

From its inception in the early 1950s through the end of the Cold War in the early 1990s, the Savannah River Site (SRS) produced nuclear materials for national defense in five reactors. Additionally, irradiated reactor fuel and target tubes were dissolved in nitric acid to recover plutonium and uranium using the PUREX (Plutonium Uranium Reduction EXtraction) process. Liquid waste from these chemical separations processes was then stored onsite in 51 underground tanks. Eight waste storage tanks have been operationally closed (i.e. cleaned and grouted) and the remaining tanks hold a mixture of liquids, insoluble solids, and precipitated salts (SRMC-LWP-2022-00001), the latter generated by evaporating water from the liquid waste. Waste is currently being retrieved from tanks and separated into 1) high-radioactivity, low-volume, and 2) low-radioactivity, high-volume components, principally through the Salt Waste Processing Facility (SWPF) (SRMC-LWP-2023-00001). The former waste stream is vitrified in the Defense Waste Processing Facility (DWPF), stored onsite, and destined for offsite disposal in a deep geologic repository. The latter stream is mixed with dry cementitious materials in the Saltstone Production Facility (SPF) and the wet slurry placed in onsite Saltstone Disposal Units (SDUs) within the Saltstone Disposal Facility (SDF), where it hardens into a cement waste form termed saltstone. A low-infiltration surface cover system will be placed over the SDF at closure, where SDUs will then be in the subsurface post-closure (SRR-CWDA-2019-00001).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Physical Modeling of Coupled Thermohydraulic Behavior of Compacted MX80 Bentonite during Heating

Abstract This paper describes a tank-scale test setup and experimental methodology used to investigate coupled heat transfer and water flow processes during heating of compacted MX80 bentonite to high temperatures. Specifically, a temperature of 200°C was maintained by a cylindrical heating element at the center of a compacted bentonite layer containing an array of temperature, dielectric, and relative humidity sensors. In addition to providing an evaluation of the spatiotemporal variations in temperature, relative humidity, degree of saturation, and global volume, the coupled thermohydraulic properties of the bentonite were assessed. A wetting front was initially observed to move away from the central heater, followed by a drying process until reaching thermohydraulic equilibrium. The soil–water retention curve (SWRC) of the bentonite followed a wetting scanning path before following the primary drying path exhibiting a shift in water retention with elevated temperature. Results from the tank-scale test can be used for validation of numerical simulations of drying processes in the engineered barrier system of a high-level radioactive waste geological disposal repository and confirm that a temperature-dependent hysteretic SWRC with scanning paths is required to accurately capture the bentonite response.

Engineering↗

toughio v1.2.2

TOUGH is a widely used general purpose numerical simulator designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media, which has been applied to many real-world problems such as underground geological storage, geothermal reservoir engineering or nuclear waste disposal, to name a few. `toughio` is a lightweight, object-oriented and vectorized Python library that aims to provide user-friendly routines to facilitate pre- and post-processing of a TOUGH simulation. Although its primary target is the latest version TOUGH3, `toughio` can also be used to read/write TOUGH2 input/output files. Currently, to the best of our knowledge, only PyTOUGH offers an exhaustive number of features to carry out a complete TOUGH simulation using a scripting language. `toughio` and PyTOUGH share the same objectives, yet with different approaches and philosophies. On the one hand, a PyTOUGH mesh is represented as a MULGRAPH geometry where elements can be unstructured horizontally but only layered vertically. It also checks that the model parameters are valid to ensure that a TOUGH simulation runs smoothly. On the other hand, although it provides basic meshing features, `toughio` mostly relies on common third-party softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT) to generate the mesh by importing and converting it to a TOUGH mesh, which also conveniently facilitates the coupling of TOUGH with any other simulator that also supports the same mesh formats. Besides, `toughio` is designed to be just a TOUGH input and output file reader/writer. It assumes that the user knows what he is doing and does not verify the validity of the model parameters (`toughio` does verify the functional arguments).

Luu, Keurfon↗

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗