A hybrid boundary element‐finite element analysis procedure for fluid flow simulation in fractured rock masses

Research output: Contribution to journalArticle

37 Citations (Scopus)

Abstract

Extension of numerical techniques to the analysis of fissure flow in three‐dimensional rock masses of realistic complexity and extent constitutes an important facility in civil engineering and resource exploitation practice. Migration velocities of contaminants and fluid discharge will always be higher for a true three‐dimensional network over a two‐dimensional representation in that the effect of critical, out‐of‐plane intersections may be correctly accommodated. Revised direct boundary element formulations are developed that are capable of discretizing individual fissure discs and their intersections with adjacent fissures. Discretization coverage, by definition, is limited to the fissure edges and internal intersections, with this factor alone being a major advantage of the technique over the more conventionally utilized domain methods. Appropriate manipulation of the resulting set of equations is shown to yield a fully populated, positive definite, symmetric tensor representing the geometric conductivity of a single fissure disc. The retained degrees‐of‐freedom for the ‘super element’ are purely in terms of the fissure intersections with a minimum of 1 degree‐of‐freedom required per intersection. Global matrix assembly and solution is accomplished by standard finite element techniques, the global matrix being, in general, sparsely populated. The procedure is ideally suited to micro‐computer solution in that a reduced degree‐of‐freedom system is obtained from a much larger and computationally, intractable system. The advantages of boundary solution procedures are realized with minimal data input requirements and effective representation of high potential gradients. The sparsely populated and symmetric from of the global matrix retains one of the more favourable assets of the secondary finite element formulation.

Original languageEnglish (US)
Pages (from-to)569-584
Number of pages16
JournalInternational Journal for Numerical and Analytical Methods in Geomechanics
Volume10
Issue number6
DOIs
StatePublished - Jan 1 1986

Fingerprint

Insulator Elements
Flow simulation
fissure
fluid flow
Flow of fluids
Rocks
rock
simulation
matrix
Civil engineering
Discharge (fluid mechanics)
Microcomputers
civil engineering
Tensors
analysis
Impurities
conductivity
Fluids
fluid
pollutant

All Science Journal Classification (ASJC) codes

  • Computational Mechanics
  • Materials Science(all)
  • Geotechnical Engineering and Engineering Geology
  • Mechanics of Materials

Cite this

@article{317289ce29544922bed3fc24c9ddb7be,
title = "A hybrid boundary element‐finite element analysis procedure for fluid flow simulation in fractured rock masses",
abstract = "Extension of numerical techniques to the analysis of fissure flow in three‐dimensional rock masses of realistic complexity and extent constitutes an important facility in civil engineering and resource exploitation practice. Migration velocities of contaminants and fluid discharge will always be higher for a true three‐dimensional network over a two‐dimensional representation in that the effect of critical, out‐of‐plane intersections may be correctly accommodated. Revised direct boundary element formulations are developed that are capable of discretizing individual fissure discs and their intersections with adjacent fissures. Discretization coverage, by definition, is limited to the fissure edges and internal intersections, with this factor alone being a major advantage of the technique over the more conventionally utilized domain methods. Appropriate manipulation of the resulting set of equations is shown to yield a fully populated, positive definite, symmetric tensor representing the geometric conductivity of a single fissure disc. The retained degrees‐of‐freedom for the ‘super element’ are purely in terms of the fissure intersections with a minimum of 1 degree‐of‐freedom required per intersection. Global matrix assembly and solution is accomplished by standard finite element techniques, the global matrix being, in general, sparsely populated. The procedure is ideally suited to micro‐computer solution in that a reduced degree‐of‐freedom system is obtained from a much larger and computationally, intractable system. The advantages of boundary solution procedures are realized with minimal data input requirements and effective representation of high potential gradients. The sparsely populated and symmetric from of the global matrix retains one of the more favourable assets of the secondary finite element formulation.",
author = "Derek Elsworth",
year = "1986",
month = "1",
day = "1",
doi = "10.1002/nag.1610100603",
language = "English (US)",
volume = "10",
pages = "569--584",
journal = "International Journal for Numerical and Analytical Methods in Geomechanics",
issn = "0363-9061",
publisher = "John Wiley and Sons Ltd",
number = "6",

}

TY - JOUR

T1 - A hybrid boundary element‐finite element analysis procedure for fluid flow simulation in fractured rock masses

AU - Elsworth, Derek

PY - 1986/1/1

Y1 - 1986/1/1

N2 - Extension of numerical techniques to the analysis of fissure flow in three‐dimensional rock masses of realistic complexity and extent constitutes an important facility in civil engineering and resource exploitation practice. Migration velocities of contaminants and fluid discharge will always be higher for a true three‐dimensional network over a two‐dimensional representation in that the effect of critical, out‐of‐plane intersections may be correctly accommodated. Revised direct boundary element formulations are developed that are capable of discretizing individual fissure discs and their intersections with adjacent fissures. Discretization coverage, by definition, is limited to the fissure edges and internal intersections, with this factor alone being a major advantage of the technique over the more conventionally utilized domain methods. Appropriate manipulation of the resulting set of equations is shown to yield a fully populated, positive definite, symmetric tensor representing the geometric conductivity of a single fissure disc. The retained degrees‐of‐freedom for the ‘super element’ are purely in terms of the fissure intersections with a minimum of 1 degree‐of‐freedom required per intersection. Global matrix assembly and solution is accomplished by standard finite element techniques, the global matrix being, in general, sparsely populated. The procedure is ideally suited to micro‐computer solution in that a reduced degree‐of‐freedom system is obtained from a much larger and computationally, intractable system. The advantages of boundary solution procedures are realized with minimal data input requirements and effective representation of high potential gradients. The sparsely populated and symmetric from of the global matrix retains one of the more favourable assets of the secondary finite element formulation.

AB - Extension of numerical techniques to the analysis of fissure flow in three‐dimensional rock masses of realistic complexity and extent constitutes an important facility in civil engineering and resource exploitation practice. Migration velocities of contaminants and fluid discharge will always be higher for a true three‐dimensional network over a two‐dimensional representation in that the effect of critical, out‐of‐plane intersections may be correctly accommodated. Revised direct boundary element formulations are developed that are capable of discretizing individual fissure discs and their intersections with adjacent fissures. Discretization coverage, by definition, is limited to the fissure edges and internal intersections, with this factor alone being a major advantage of the technique over the more conventionally utilized domain methods. Appropriate manipulation of the resulting set of equations is shown to yield a fully populated, positive definite, symmetric tensor representing the geometric conductivity of a single fissure disc. The retained degrees‐of‐freedom for the ‘super element’ are purely in terms of the fissure intersections with a minimum of 1 degree‐of‐freedom required per intersection. Global matrix assembly and solution is accomplished by standard finite element techniques, the global matrix being, in general, sparsely populated. The procedure is ideally suited to micro‐computer solution in that a reduced degree‐of‐freedom system is obtained from a much larger and computationally, intractable system. The advantages of boundary solution procedures are realized with minimal data input requirements and effective representation of high potential gradients. The sparsely populated and symmetric from of the global matrix retains one of the more favourable assets of the secondary finite element formulation.

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

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

U2 - 10.1002/nag.1610100603

DO - 10.1002/nag.1610100603

M3 - Article

VL - 10

SP - 569

EP - 584

JO - International Journal for Numerical and Analytical Methods in Geomechanics

JF - International Journal for Numerical and Analytical Methods in Geomechanics

SN - 0363-9061

IS - 6

ER -