GEOS
GEOS copied to clipboard
Use aperture table in poromechanics with conforming fractures
Main thing:
replace
hydraulicAperture[k] = minimumHydraulicAperture[k] + aperture[k];
with
hydraulicAperture[k] = contactWrapper.computeHydraulicAperture( aperture[k], dHydraulicAperture_dNormalJump );
for poromechanics with conforming fractures.
At least this way there will be no negative aperture and no negative volume.
Small things:
- Separately account for new slip in
computeFractureStateStatistics
to see new segments slipped. - Output normal displacement tolerance, sliding tolerance, normal traction tolerance when
logLevel>=2
. - For
logLevel>=10
- display fracture state changes with debug information. - Small fixes for
deltaVolume
logic: assign it to 0 once applied into volume (to avoid double counting), update mass with updated volume. Those are important for fractures after poromechanics initialization that changes the volume of the fractures.