Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

In a repository for radioactive waste hosted in a clay formation, hydrogen and other gases may be generated due to the corrosion of metallic materials under anoxic conditions, the radioactive decay of waste and the radiolysis of water. If the gas production rate exceeds the gas diffusion rate within the pores of the clay, a discrete gas phase will form and accumulate until its pressure becomes large enough to exceed the entry pressure of the surrounding material, at which point dilatant, advective flow of gas is expected to occur. The purpose of Task An under DECOVALEX-2019 is to better represent the processes governing the advective movement of gas in both low-permeability argillaceous repository host rocks and clay-based engineered barriers within numerical codes. In this paper special attention is given to the mechanisms controlling gas entry, flow and pathway sealing and their impact on the performance of the engineered clay barrier. Previous work suggests gas flow is accompanied by the creation of dilatant pathways whose properties change temporally and spatially within the medium. Thus, four new types of approaches have been developed: (i) standard two-phase flow models (continuous techniques) incorporating a range of different mechanical deformation behaviours, (ii) enhanced two-phase flow models in which fractures are embedded within a plastic material (continuous techniques) or incorporated into the model using a rigid-body-spring network (discrete approaches), (iii) a single-phase model incorporating a creep damage function in which only gas flow is considered, and (iv) a conceptual approach used to examine the chaotic nature of gas flow. The outputs from these different approaches are compared. This is an essential step as the choice of modelling approach strongly impacts the representation and prediction of gas flow in a future repository. In addition, experience gained through this task is of direct relevance to other clay-based engineering issues where immiscible gas flow is a consideration including hydrocarbon migration, carbon capture and storage, shale gas and landfill design. This paper summarises the outcomes of work in Task A conducted between May 2016 and May 2019 and provides a brief overview of the experimental data and a synthesis of the work of the participating modelling teams.