Dear Dmitry,
if I remember correctly, the sequential MPFA-L scheme requires that each
vertex is surrounded by 4 quadrilaterals (in 2D).
If this is not the case, the construction of interaction volumes fails.
Is there any reason why you are using the sequential method?
Best regards,
Martin
On 23.03.20 18:42, Dmitry Pavlov wrote:
Dear Kilian,
1.) You should be able to use more recent versions of gmsh if you go to File ->
Export. If you name your grid *.msh
there will be an option to chose the older Version II ASCII format, readable by
Dumux.
Thank you!
2.) Have you tried an unstructured grid with quadrilaterals? According to the
header where the error comes from,
only quadrilaterals are supported (at least, this is what some comments in the
code suggest). You can force gmsh to
use only quads (also unstructured ones).
I just tried an unstructured quad grid and got the same error.
However, a gmsh-generated structured quad grid went fine, which means
that there is no trouble with gmsh itself.
Good news: using FVPressureTwoPAdaptive instead of
FvMpfaL2dPressureTwoP allowed me to use an gmsh-generated unstructured
grid (both triangular and quadrilateral). I am wondering whether
MPFA-L in DuMux "officially" requires a structured quad grid.
Best regards,
Dmitry
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
