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
mdg (MixedDimensionalGrid) –
- 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