A Semi-Continuum Model for Numerical Simulations of Mass Transport in 3-D Fractured Rock Masses

Guowei Ma, Tuo Li, Yang Wang, Yun Chen

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Natural discontinuities are the major concern when considering mass transport in engineered barrier systems and their host rocks. Numerical simulations of highly fractured geological formations are limited because of the contradiction between results accuracy and computational costs. To alleviate such a contradiction, this study proposes an improved fracture continuum method to simulate the radioactive spreading in a complex 3-D fracture system. With Boolean operations and the unified pipe-network method, discontinuities are mapped on structured subdomains and standardized to equivalent paths. Moreover, adaptive mesh refinement is utilized to ease the complexity further. We verify the accuracy of this method in two metric cases, and results show that perfect agreement is achieved with analytical solutions. This method demonstrates its applicability to the simulation of a nuclear leak in a repository for high-level radioactive waste. Effects of the maximum refinement level are discussed. For a room-scale problem, flow rates and mass fluxes on boundary surface converge to stable values when the maximum refinement level is larger than 4. An extensive case with complex fracture networks is modeled and compared with the conventional finite-difference method. The proposed method is capable of conducting robust results with significantly lower computational complexity and negligible errors by avoiding ill-conditioned mesh elements.

Original languageEnglish
JournalRock Mechanics and Rock Engineering
DOIs
Publication statusE-pub ahead of print - 26 Aug 2019

Fingerprint

mass transport
Mass transfer
Rocks
Computer simulation
rock
simulation
Radioactive wastes
discontinuity
Finite difference method
Computational complexity
Pipe
Flow rate
fracture network
finite difference method
repository
radioactive waste
host rock
pipe
method
Costs

Cite this

@article{84567b9932cf47eb9c1900db6593b6f3,
title = "A Semi-Continuum Model for Numerical Simulations of Mass Transport in 3-D Fractured Rock Masses",
abstract = "Natural discontinuities are the major concern when considering mass transport in engineered barrier systems and their host rocks. Numerical simulations of highly fractured geological formations are limited because of the contradiction between results accuracy and computational costs. To alleviate such a contradiction, this study proposes an improved fracture continuum method to simulate the radioactive spreading in a complex 3-D fracture system. With Boolean operations and the unified pipe-network method, discontinuities are mapped on structured subdomains and standardized to equivalent paths. Moreover, adaptive mesh refinement is utilized to ease the complexity further. We verify the accuracy of this method in two metric cases, and results show that perfect agreement is achieved with analytical solutions. This method demonstrates its applicability to the simulation of a nuclear leak in a repository for high-level radioactive waste. Effects of the maximum refinement level are discussed. For a room-scale problem, flow rates and mass fluxes on boundary surface converge to stable values when the maximum refinement level is larger than 4. An extensive case with complex fracture networks is modeled and compared with the conventional finite-difference method. The proposed method is capable of conducting robust results with significantly lower computational complexity and negligible errors by avoiding ill-conditioned mesh elements.",
keywords = "Adaptive mesh refinement, Boolean operations, Fractured rock masses, Radioactive transport, Semi-continuum model, Unified pipe-network method",
author = "Guowei Ma and Tuo Li and Yang Wang and Yun Chen",
year = "2019",
month = "8",
day = "26",
doi = "10.1007/s00603-019-01950-1",
language = "English",
journal = "Journal Rock Mechanics and Rock Engineering",
issn = "0723-2632",
publisher = "Springer",

}

A Semi-Continuum Model for Numerical Simulations of Mass Transport in 3-D Fractured Rock Masses. / Ma, Guowei; Li, Tuo; Wang, Yang; Chen, Yun.

In: Rock Mechanics and Rock Engineering, 26.08.2019.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A Semi-Continuum Model for Numerical Simulations of Mass Transport in 3-D Fractured Rock Masses

AU - Ma, Guowei

AU - Li, Tuo

AU - Wang, Yang

AU - Chen, Yun

PY - 2019/8/26

Y1 - 2019/8/26

N2 - Natural discontinuities are the major concern when considering mass transport in engineered barrier systems and their host rocks. Numerical simulations of highly fractured geological formations are limited because of the contradiction between results accuracy and computational costs. To alleviate such a contradiction, this study proposes an improved fracture continuum method to simulate the radioactive spreading in a complex 3-D fracture system. With Boolean operations and the unified pipe-network method, discontinuities are mapped on structured subdomains and standardized to equivalent paths. Moreover, adaptive mesh refinement is utilized to ease the complexity further. We verify the accuracy of this method in two metric cases, and results show that perfect agreement is achieved with analytical solutions. This method demonstrates its applicability to the simulation of a nuclear leak in a repository for high-level radioactive waste. Effects of the maximum refinement level are discussed. For a room-scale problem, flow rates and mass fluxes on boundary surface converge to stable values when the maximum refinement level is larger than 4. An extensive case with complex fracture networks is modeled and compared with the conventional finite-difference method. The proposed method is capable of conducting robust results with significantly lower computational complexity and negligible errors by avoiding ill-conditioned mesh elements.

AB - Natural discontinuities are the major concern when considering mass transport in engineered barrier systems and their host rocks. Numerical simulations of highly fractured geological formations are limited because of the contradiction between results accuracy and computational costs. To alleviate such a contradiction, this study proposes an improved fracture continuum method to simulate the radioactive spreading in a complex 3-D fracture system. With Boolean operations and the unified pipe-network method, discontinuities are mapped on structured subdomains and standardized to equivalent paths. Moreover, adaptive mesh refinement is utilized to ease the complexity further. We verify the accuracy of this method in two metric cases, and results show that perfect agreement is achieved with analytical solutions. This method demonstrates its applicability to the simulation of a nuclear leak in a repository for high-level radioactive waste. Effects of the maximum refinement level are discussed. For a room-scale problem, flow rates and mass fluxes on boundary surface converge to stable values when the maximum refinement level is larger than 4. An extensive case with complex fracture networks is modeled and compared with the conventional finite-difference method. The proposed method is capable of conducting robust results with significantly lower computational complexity and negligible errors by avoiding ill-conditioned mesh elements.

KW - Adaptive mesh refinement

KW - Boolean operations

KW - Fractured rock masses

KW - Radioactive transport

KW - Semi-continuum model

KW - Unified pipe-network method

UR - http://www.scopus.com/inward/record.url?scp=85071414149&partnerID=8YFLogxK

U2 - 10.1007/s00603-019-01950-1

DO - 10.1007/s00603-019-01950-1

M3 - Article

JO - Journal Rock Mechanics and Rock Engineering

JF - Journal Rock Mechanics and Rock Engineering

SN - 0723-2632

ER -