porepy.numerics.fracture_deformation.propagate_fracture module

Propagation of fractures. Much in common with (and reuse of) split_grid. For now assumes:

single fracture

When this assumption is relieved, some (re)structuring will be needed. The structure for multi-fracture propagation may possibly strongly resemble that of split_grid.

WARNING: This should be considered experimental code, which cannot be assumed to be bug free.

propagate_fractures(mdg, faces)[source]

mdg - Mixed-dimensional grid with matrix and fracture grids. faces_h - list of list of faces to be split in the highest-dimensional

grid. The length of the outer list equals the number of fractures. Each entry in the list is a list containing the higher-dimensional indices of the faces to be split for the extension of the corresponding fracture.

Changes to grids done in-place. The call changes:

Geometry and connectivity fields of the two grids involved. The face_cells mapping between them Their respective face tags.

Also adds the following to subdomain data dictionaries:

new_cells and new_faces tags, for use in e.g. local discretization updates. partial_update, a boolean flag indicating that the grids have been updated.

Parameters
Return type

None

propgation_angle(K)[source]

Compute the angle of propagation from already computed SIFs. The computation is done in the local coordinate system of the fracture tips, and positive angles correspond to. Intended for a single fracture grid.

Parameters

K – array of stress intensity factors, with mode along first axis and face along second.

Returns

array (number of tip faces, ) of propagation angles in radians.

Return type

phi