Pavel Tomin
Pavel Tomin
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...
1. Pull out hard-codded params to be able to modify them from input: ``` - strongestTol = 1e-8 (default value) - adaptiveGamma = 0.1 - adaptiveExponent = 1.0 ``` 2....
Case: https://github.com/GEOS-DEV/MAELSTROM/tree/master/usecases/francois/HI24L/singlePhasePoromechanics Input file: HI24L_single_phase_hybrid_fvm.xml Does not run with current version: [slurm-161984.txt](https://github.com/GEOS-DEV/GEOS/files/12563743/slurm-161984.txt) because of linear solver issue.
`postRestartInitialization` function needs to be implemented similar to https://github.com/GEOS-DEV/GEOS/blob/db7cc5d62fa962a56a0eda5d3d27f8fa5ef36db2/src/coreComponents/physicsSolvers/surfaceGeneration/SurfaceGenerator.cpp#L424C1-L425C2
one is when all component concentrations are zeros - division by zero and another is when maxTime is not specified for events - there is overflow when huge double is...
not sure how really useful it is, please let me know
also make sure non-public functions are under protected
should go after https://github.com/GEOS-DEV/GEOS/pull/3383 and https://github.com/GEOS-DEV/GEOS/pull/3386 - works equivalent to `MultiphasePoromechanicsReservoir` and `SinglePhasePoromechanicsReservoir`, just more convenient to switch between FIM and SEQ within the same solver - can use the...
should go after https://github.com/GEOS-DEV/GEOS/pull/3383