porepy.numerics.fracture_deformation.conforming_propagation module
Module indended for combining fracture propagation with complex multiphysics, as represented in the model classes.
- WARNING: This should be considered experimental code and should be used with
extreme caution. In particular, the code likely contains bugs, possibly of a severe character. Moreover, simulation of fracture propagation may cause numerical stability issues that it will likely take case-specific adaptations to resolve.
- WARNING 2: At the moment there are assumptions of purely tensional propagation.
Violation of this assumption is not recommended with the current implementation.
- Contains:
ConformingFracturePropagation - class to be used together with a pp Model for propagation simulation.
- Required additional mechanical parameters:
poisson_ratio (data_primary) shear_modulus (data_primary) SIFs_critical (data_secondary)
- Literature:
Thomas et al. 2020: Growth of three-dimensional fractures, arrays, and networks in brittle rocks under tension and compression Nejati et al, 2015: On the use of quarter-point tetrahedral finite elements in linear elastic fracture mechanics Richard et al. 2005: Theoretical crack path prediction
- class ConformingFracturePropagation(params)[source]
Bases:
FracturePropagationClass for fracture propagation along existing faces in the higher-dimensional grid.
The propagation criteria is based on stress intensity factors, computed from a displacement correlation approach.
Should be used in combination with a Model class, i.e. assumptions on methods and fields are made.
- WARNING: This should be considered experimental code and should be used with
extreme caution. In particular, the code likely contains bugs, possibly of a severe character. Moreover, simulation of fracture propagation may cause numerical stability issues that it will likely take case-specific adaptations to resolve.
The code structure for fracture propagation cannot be considered fixed, and it may be fundamentally restructured at unknown points in the future. If you use this functionality, please notify the maintainers (Eirik.Keilegavlen@uib.no), so that we may keep your usecases in mind if a major overhaul of the code is undertaken.
- evaluate_propagation()[source]
Evaluate propagation for all fractures based on the current solution.
Computes SIFs using the Displacement Correlation method described in Nejati et al. based on the displacement jumps of the previous iterate. Then, propagation onset and angles are evaluated as described in Thomas et al.
Returns
None.
- Return type
None