Pavel Tomin

Unify convergence check for lagrange-contact solvers

3

Main thing: unify convergence check for lagrange-contact solvers - do not assign max force to zero in `SinglePhasePoromechanicsConformingFractures::assembleElementBasedContributions` to ensure the same normalization across the solvers and FIM/SEQ - do...

Separate stick, slip, open parts of contact residual and try to normalize them properly

6

Separately compute contact residuals for elements in stick, slip, and open states and normalize: - Stick: normalize locally rhs by element area and compute L2-norm of that - Slip and...

Simplify poromechanics initialization setup

1

Simplify definition of poromechanics initialization in the input file. Instead of: ``` ... ``` one now needs to specify only: ``` ... ``` No results changes are expected.

Preparatory work for fractures + wells

2

For now, assuming wells are not directly connected to the fractures. In this version, only `SinglePhasePoromechanicsConformingFracturesReservoir` solver can run, although I am having hard time to build a well-behaving example....

Freeze porosity when doing poromechanics initialization

6

This PR attempts to fix the following issue: - during poromechanics initialization step the porosity is being updated which affects the mechanical solution trough changes in body force - once...

Make sequential poromechanics to conserve mass

1

Has to go after https://github.com/GEOS-DEV/GEOS/pull/2906 and https://github.com/GEOS-DEV/GEOS/pull/2897 With changes in this PR explicit coupling seems to work fine and I am able to match FIM reference for SEAM case: ...

pknViscosityDominated_benchmark does not run

12

File: inputFiles\hydraulicFracturing\pknViscosityDominated_benchmark.xml Output: [run.604958.txt](https://github.com/GEOS-DEV/GEOS/files/13395574/run.604958.txt)

HDF5 errors when running with MPI

5

Running this case https://geosx-geosx.readthedocs-hosted.com/en/latest/docs/sphinx/basicExamples/multiphaseFlowWithWells/Example.html Serial runs just fine, when run with mpirun: `OMP_NUM_THREADS=1 mpirun -n 2 /data/rpo_ptls/GEOSX/residual_flash/GEOS/GEOS/build-CPU-OPTO2-Hypre-GCC_10.2.0-ompi_hpcx-OMP-relwithdebinfo/bin/geosx -i deadOilEgg_benchmark.xml` HDF5 complains at the end of simulation when timeHistoryOutput is triggered:...

Make sure all documentation examples are included in integrated tests

1

Documentation examples should run as expected, this is first thing that users start to look at when getting familiar with GEOS. Therefore all of them must be in integrated tests...

Single Fracture Under Shear Compression from Advanced Examples does not work with MPI

3

``` [ptls@lxsphou0099 lagrangianContactMechanics]$ OMP_NUM_THREADS=1 mpirun -n 2 /data/rpo_ptls/GEOSX/residual_flash/GEOS/GEOS/build-CPU-OPTO2-Hypre-GCC_10.2.0-ompi_hpcx-OMP-relwithdebinfo/bin/geosx -i ContactMechanics_SingleFracCompression_benchmark.xml Num ranks: 2 Max threads: 1 GEOSX version: 0.2.0 (develop, sha1: 0eb9b90c2) - c++ compiler: gcc 10.4.0 - openmp version:...