K-regime Hydro Frac in PhaseField

See merge request ogs/ogs!4299

KeitaYoshioka authored 2 years ago and Lars Bilke committed 2 years ago

ctest

ff

code change

[cmake] Build PhaseField on PETSc configs only.

See merge request ogs/ogs!4335

Added reflection functionality for setting IP data

See merge request ogs/ogs!4316

All current PhaseField use cases need PETSc. Excluding it from serial
build will speed up CI a bit.

[PL/BC] Reimplement Robin-type BC

See merge request ogs/ogs!4302

Removed <...stream> includes from header files.

See merge request ogs/ogs!4330

Load the wedge test file in Paraview, rotate X-90,
resample with the 2D mesh, store the 2D mesh result.

Both, the 2D in cylinder coordinates test and
the corresponding 3D wedge in cartesian coordinates
are updated. Verified against analytical solution
given in both project files for initial conditions.

New references are provided for numerical comparison.

Since here two mesh resolutions were already available
the grid convergence is roughly quadratic for cylinder
coordinates and cartesian coordinates cases as expected.

Both, the 2D in cylinder coordinates test and
the corresponding 3D wedge in cartesian coordinates
are updated. Verified against analytical solution
given in both project files for initial conditions.

New references are provided for numerical comparison.

Also update extracted surface for the same test in 2D
with axial symmetry.

Both project files were unified (time stepping, output,
and layout).

A method of manufactured solutions is used to verify implementation of
the Robin-type boundary condition. Grid convergence was assessed and the
observed order of convergence is 1.99, which, as expected, is close
enough to 2.

In the same setup a Neumann and Dirichlet-type boundary conditions can
(and were) tested with same convergence order as the Robin-type BC.

The test is testing both non-linear solver schemes, Newton and Picard.
With the new implementation of the Robin-type BC the convergence rates
of the non-linear solver are satisfying, usually 2 iterations per time
step for both schemes.

The former implementation did not produce expected results for 2D case
(and probably for 3D case). Rewriting and implementing the Jacobian for
the Newton scheme. Keeping K matrix for the Picard scheme as applying
the BC to the r.h.s. only significantly increases number of non-linear
iterations because the problem becomes non-linear.

[ctest/partmesh] Use correct ip data for one test

See merge request ogs/ogs!4328

Added reflection functionality for ip data output

See merge request ogs/ogs!4317

More concise implementation.