We reuse the grid creation facility of the facet coupling framework and
require grid files that have the information on the fracture facets in
them. This enables us to reuse the codimension one grid adapter of the
facet coupling framework to identify which matrix element facets
coincide with a fracture.

This commit also introduces a test on the same geometry as in the 2.12
release.

TODO: add test for quads and tetrahedra

Test in test/geomechanics/elastic.

Bernd Flemisch authored 6 years ago and Timo Koch committed 6 years ago

There are enough other ni convection tests in
1d and there's the 1pni conduction test in 2d.

* add source (gravity) term to residual, also when face has Neumann BCs
* correct coupling condition

* reference solution needs to be updated because test uses time-depended Dirichlet BC
* old behaviour was actually wrong

Melanie Lipp authored 6 years ago and Kilian Weishaupt committed 6 years ago

[freeflow][outflow BC] Introduce that there is no gradient normal to the outflow boundary in cases where the outflow boundary is the normal face with respect to the assembled scvf.

Bernd Flemisch authored 6 years ago and Timo Koch committed 6 years ago

Also fix the test itself in the following sense. By default, in the
box method, material law parameters are evaluated at scv centers,
while primary variables are evaluated at vertices. If no interface
condition is used, the secondary variables handed to the Vtk
writer might now differ depending on the order in which the scvs
are evaluated for each vertex. In order to saty consistent, this
test evaluates the material law parameters at the vertex position.

* test was passing but vtk #39 was compared although there are 40 steps

This makes use of the introduced fluid state adapter for mapping to the right
phase and component indices.

This fixes:
- calculation of the liquid mixture density (for the case of it being composition-depent)
- make fugacity coefficient of NaCl zero (we assume NaCl to stay in the
water phase anyway throughout this fluid system)
- adapt reference solution

previous implementation could potentially result in unwanted forward calls in the fluid system

Sina Ackermann authored 6 years ago and Kilian Weishaupt committed 6 years ago

* consistent ordering of component/phase names: x^componentName_phaseName
* consistent use of underline: S_w, p_w, etc.
* change respective data field names in reference solutions

Sina Ackermann authored 6 years ago and Timo Koch committed 6 years ago

* horizontal flow (evaporation)

Sina Ackermann authored 6 years ago and Timo Koch committed 6 years ago

* vertical flow (conservation problem)

Sina Ackermann authored 6 years ago and Timo Koch committed 6 years ago

* vertical flow
* horizontal flow

[test][multidomain][embedded] Add test and references for 1p1p 2d3d coupling (embedded fractures example)