Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2022

Fracture toughness and crack propagation in LX-17 and PBX 9502 insensitive high explosives

High explosive materials, and in particular insensitive high explosive (IHE) materials, are of significant interest. TATB-based LX-17 and PBX 9502 are two IHE’s and mechanical properties are an important metric of performance. In this work we present, for the first time, a study of fracture as a function of temperature and displacement rate and starter crack length on billet pressed IHE materials using a standard three point bend technique. We demonstrated that the technique can yield consistent results over time by measuring samples from the same lots 4 years later and observed that the fracture toughness of LX-17 and PBX 9502 are similar. We also determined the temperature dependence for different starter crack lengths and found that the fracture toughness generally decreases with larger starter crack sizes and temperatures. For the LX-17, where a wide temperature range of -60 C to 50 C was investigated, the impact of temperature on fracture toughness between -20 C and 23 C was less than expected. For the PBX 9502 we observed the difference in fracture toughness values between -20 C and 23 C did depend on the particular lot, and also there was a tendency of the fracture toughness to scale with the tensile strength. We also investigated the effect of a lower displacement rate and observed that while the impact on the fracture toughness was relatively small, the samples surprisingly did not break (fail) after exceeding the peak load. We took advantage of this to perform x-ray computed tomography (CT) of the samples to image crack propagation in the intact samples. We observed that cracks were generally straight and unaffected by the microstructure in the LX-17, in contrast to the PBX 9502 where the direction of the cracks would frequently change and were likely influenced by the underlying microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Role of water in fracture of modified silicate glasses

Decarbonizing the glass industry requires alternative melting technology, as current industrial melting practices rely heavily on fossil fuels. Hydrogen has been proposed as an alternative to carbon-based fuels, but the ensuing consequences on the mechanical behavior of the glass remain to be clarified. A critical distinction between hydrogen and carbon-based fuels is the increased generation of water during combustion, which raises the equilibrium solubility of water in the melt and alters the behavior of the resulting glass. A series of five silicate glasses with 80% silica and variable [Na 2 O]/([H 2 O] + [Na 2 O]) ratios were simulated using molecular dynamics to elucidate the effects of water on fracture. Several fracture toughness calculation methods were used in combination with atomistic fracture simulations to examine the effects of hydroxyl content on fracture behavior. Here, this study reveals that the crack propagation pathway is a key metric to understanding fracture toughness. Notably, the fracture propagation path favors hydrogen sites over sodium sites, offering a possible explanation of the experimentally observed effects of water on fracture properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2023

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2020

Microstructure and fracture toughness characterization of three 9Cr ODS EUROFER steels with different thermo-mechanical treatments

Ferritic martensitic ODS steels are one of the candidate structural materials for future Gen-IV nuclear fission and fusion reactors. Here, the dependence of fracture toughness on microstructure was investigated by comparing three 9Cr ODS EUROFER steels manufactured through different thermo-mechanical processing routes. Quasi-static fracture toughness testing was performed with sub-sized C(T) specimens and microstructural characterization was carried out using scanning electron microscopy, electron backscatter diffraction, and transmission electron microscopy. It was found that at lower test temperatures (-100 – 22 °C), the fracture toughness was primarily controlled by crack initiation at sub-micron particles and by production of secondary cracks during fracture. At higher temperatures (above 100 °C), fracture toughness was predominantly controlled by the matrix ductility and the grain boundary strength with a relatively ductile coarse-grained alloy demonstrating higher fracture toughness compared to high-strength fine-grained alloys. These results and discussion show that variations in thermomechanical treatments can produce significant differences in microstructure and fracture toughness behavior of ferritic martensitic ODS steels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fracture behavior of irradiation induced nanocrystalline UO 2 studied by in-situ mechanical testing in transmission electron microscopy

Uranium Dioxide (UO 2 ) is widely used as a fuel in current light water reactors (LWRs). Upon accumulation of radiation damage, LWR UO 2 fuel pellets start to develop a different microstructure at the pellet periphery when fuel burnup exceeds 45–50 GWd/tHM. The resulting porous, nanocrystalline microstructure is one of the most prominent microstructural changes occurring in such fuel. Its fracture mechanisms, which causes fuel fine fragmentation, could impact safety limits when the cladding breaches. Direct measurements of these properties are challenging, therefore a surrogate obtained via ion irradiation can be used. In this study, multiple microcantilevers were fabricated by focused ion beam from both fresh UO 2 and UO 2 irradiated with 84 MeV Xe 26+ ions to a peak dose of 1357 displacements per atom (dpa). Further, the irradiation produced a pseudo high burnup structure approximately 2 µm below the surface. In-situ nano-mechanical bending tests were conducted to investigate the fracture behavior and the effect of the surrogate UO 2 high burnup structure on local fracture properties. Fresh UO 2 fuel was observed to fracture in transgranular mode without nucleation or movement of dislocations. However, the Xe-irradiated nanocrystalline microcantilevers fractured along the grain boundaries, with no influence from the pre-existing micro-cracks in the microcantilever. Fracture toughness for this type of surrogate high burnup UO 2 structure is reported for the first time in literature. Both the fracture stress and toughness show degradation for UO 2 as a result of Xe-irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Statistical fracture behavior of doped UO 2 using a ball-on-ring equibiaxial flexure test method

Metal oxide dopants, such as titanium and chromium oxides, have garnered considerable attention for their potential to increase grain size (≥ 30 µm) in UO 2 fuel, purportedly enhancing fission gas retention during reactor operation. Fuel performance is significantly impacted by fuel fracture behavior, so it is important to understand the effects of enhanced grain size and dopant content on UO 2 fuel fracture. UO 2 pellets were doped with 0.1 wt% TiO 2 and 0.3 wt% Cr 2 O 3 to alter density and grain size. Inductively coupled plasma mass spectroscopy measured dopant levels pre- and post-sintering. X-ray diffraction revealed lattice changes and microstrain via Rietveld refinement. Field emission scanning electron microscopy determined grain sizes of approximately 30 µm for TiO 2 doping and 7 µm for Cr 2 O 3 doping. Transverse rupture strength tests were performed on over 30 samples per dataset to obtain characteristic strength and Weibull modulus. Results indicate no statistical difference in fracture strength between 0.1 wt% TiO 2 doped UO 2 and undoped UO 2 , while 0.3 wt% Cr 2 O 3 doped UO 2 exhibited a 20% decrease in fracture strength. Doped UO 2 samples also showed reduced Weibull modulus compared to undoped UO 2 , suggesting increased scatter in fracture strength. This study's findings suggest that titanium and chromium oxide doping in UO 2 , regardless of grain size, induce residual stresses, decreasing fracture strength and increasing variability in fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2021

Fracture characteristics and heat treatment of laser powder bed fusion additively manufactured GRCop-84 copper

Laser Powder Bed Fusion (LPBF) of Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, produces a fully dense high conductivity alloy with a yield strength of 500 MPa and ultimate tensile strength (UTS) of 740 MPa, superior to other competing copper alloys, and 20% elongation at fracture for material stressed perpendicular to the build direction. The high thermal stability of the Cr 2 Nb precipitate in the copper matrix reduces coarsening and maintains a 300 MPa yield, 520 MPa UTS and 26% elongation after a 900 °C, 5-h heat treatment, while a 3 h 450 °C heat treatment increases yield to 810 MPa, UTS to 970 MPa with 9% elongation at fracture, for samples stressed perpendicular to the build direction. Tensile strength anisotropy based on print direction was attributed to internal stress and columnar grain formation. Void nucleation during tensile fracture was initiated by brittle fracture of precipitate particles within the copper matrix. Fracture cusps contain matching precipitate fragment geometry on opposing sides located near the cusp center in at least 80% of fracture cusps. An optimal precipitate size of 100 nm is predicted for maximum tensile strength from precipitates on fracture surfaces, while tensile testing with varying heat treatments shows maximum strength with 100 nm and smaller precipitates. Cr2Nb precipitates are shown to transition between polycrystalline and monocrystalline structures at high temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stress-sensitivity of fracture conductivity of Tuscaloosa Marine Shale cores

Production decline from wells producing from shale rocks is typically quite rapid. One possible factor responsible for such decline is loss in fracture conductivity. This is true in the case of the Tuscaloosa Marine shale reservoir. In this work, the stress-dependent fracture conductivity of TMS core samples is investigated using five propped cores under varying stress conditions. In addition, the effects of static rock mechanical properties and rock mineralogy on conductivity decline are investigated. The mineralogy of each core sample was determined by X-ray Diffraction analysis. Fracture conductivity was measured using a Hassler-type pressure core holder with a pressure rating of 15,000 psi (103.42 MPa). Estimated rock mechanical properties showed that shale anisotropy could affect fracture conductivity. The time decline in conductivity revealed two possible decline trends, a) decline from the reduction in induced fracture width, and b) decline from the healing of developed micro-cracks. An exponential decline in fracture conductivity with increasing confinement pressure was observed with an average decline rate constant of 3.15 x 10 -10 mPa -1 . While rock mineralogy did not have any clear relationship with conductivity decline, the effects of rock mechanical properties were only important at 10.34 MPa confinement pressure and above. The observation showed that a high Young's modulus correlate to high conductivity, whereas a high Poisson's ratio correlate to low conductivity. This study provides useful observations for optimizing fracture design to enhance well productivity in the TMS.

02 PETROLEUM↗
DOE OSTI2022

Monitoring Fracture Saturation With Internal Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of Pore Water on the Fracture of Silica Sand at Particle Scale

The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.

cavitation↗
DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Sealing Mechanisms of a Fracture in Opalinus Clay as Revealed by in situ Synchrotron X-Ray Micro-Tomography

The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of confining pressure. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.

58 GEOSCIENCES↗