GEOS icon indicating copy to clipboard operation
GEOS copied to clipboard

Use aperture table in poromechanics with conforming fractures

Open paveltomin opened this issue 7 months ago • 4 comments

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.

paveltomin avatar Jun 28 '24 16:06 paveltomin