Grid Project:   Discrete Fracture Networks

Overview

Intersection

Voronoi Cells

Gallery

For more pictures, take a look at our gallery.

Design Criteria

• The mesh of the fractures are conforming at the intersection.
• The grid resolution is higher near the intersection (length scale 0.1) then gradually coarsens to the surrounding area (length scale 5).

Methods

The Grid Method successfully follows the design criteria.

Given a network of intersecting fractures, we mesh each fracture separately and merge all the fractures into a network at the end. The method uses a newly developed capability that allows adaptive refinement according to a gradient field. The field chosen for this case is a simple linear function of a distance field, with a minimum floor value. The field is computed by LaGriT commands inside a separate .mlgi, .lgi, or any text file, which is needed for 'massage2' command.

A python script is also written to control the workflow and write lagrit command files for meshing and merging.

Input files

• input_n100.m
• final_dfn.py
• intersections.inp
• params.txt
• polys.inp
• user_function.lgi

Run Sequence

python dfn_final.py params.txt 8

Project Archive Location

/scratch/sft/quanb/grid_gen/dfn_examples/example_n100/

FEHMN output files

-rw-r----- 1 quanb altern  9852836 Apr  5 14:29 tri_fracture.fehmn
-rw-r----- 1 quanb altern   783146 Apr  5 14:29 tri_fracture_material.zone
-rw-r----- 1 quanb altern   834805 Apr  5 14:29 tri_fracture_interface.zone
-rw-r----- 1 quanb altern   691112 Apr  5 14:29 tri_fracture_multi_mat.zone
-rw-r----- 1 quanb altern    83835 Apr  5 14:29 tri_fracture_outside.zone
-rw-r----- 1 quanb altern   483127 Apr  5 14:29 tri_fracture_outside_vor.area
-rw-r----- 1 quanb altern 24102832 Apr  5 14:29 tri_fracture.stor

Work Directory:

/home/quanb/src/dfn_src/final/example_n100/

Grid Team Contacts