A constrained Delaunay discretization method for adaptively meshing highly discontinuous geological media

Yang Wang, Guowei Ma, Feng Ren, Tuo Li

    Research output: Contribution to journalArticle

    26 Citations (Scopus)

    Abstract

    A constrained Delaunay discretization method is developed to generate high-quality doubly adaptive meshes of highly discontinuous geological media. Complex features such as three-dimensional discrete fracture networks (DFNs), tunnels, shafts, slopes, boreholes, water curtains, and drainage systems are taken into account in the mesh generation. The constrained Delaunay triangulation method is used to create adaptive triangular elements on planar fractures. Persson's algorithm (Persson, 2005), based on an analogy between triangular elements and spring networks, is enriched to automatically discretize a planar fracture into mesh points with varying density and smooth-quality gradient. The triangulated planar fractures are treated as planar straight-line graphs (PSLGs) to construct piecewise-linear complex (PLC) for constrained Delaunay tetrahedralization. This guarantees the doubly adaptive characteristic of the resulted mesh: the mesh is adaptive not only along fractures but also in space. The quality of elements is compared with the results from an existing method. It is verified that the present method can generate smoother elements and a better distribution of element aspect ratios. Two numerical simulations are implemented to demonstrate that the present method can be applied to various simulations of complex geological media that contain a large number of discontinuities.

    Original languageEnglish
    Pages (from-to)134-148
    Number of pages15
    JournalComputers and Geosciences
    Volume109
    DOIs
    Publication statusPublished - 1 Dec 2017

    Fingerprint Dive into the research topics of 'A constrained Delaunay discretization method for adaptively meshing highly discontinuous geological media'. Together they form a unique fingerprint.

  • Cite this