Dear Dmitry,
unfortunately, my knowledge about the sequential MPFA-L implementation is very limited but I have two remarks: 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. 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). Best regards Kilian ________________________________ Von: Dumux <[email protected]> im Auftrag von Dmitry Pavlov <[email protected]> Gesendet: Montag, 23. März 2020 15:30 An: [email protected] Betreff: [DuMuX] "boundary shape not available as interaction volume shape" with gmsh-made grid Hello, I am a new DuMux user and I am trying to simulate a 2p2c flow on a 2D rectangular area with one border being the influx and the opposite one being the outflow (an the two others being walls, of course). I used the example dumux/test/porousmediumflow/2p/sequential/test_mpfa2p.cc as the starting point of my program, using MPFA-L method. All went fine on a regular grid. Now I am willing to try an unstructured grid, too. From the docs, I figured that one recommended way of creating unstructured grids is gmsh, am I correct? .msh files obtained from gmsh 4 did not work because Dune::IOError [read:.../dune-grid/dune/grid/io/file/gmshreader.hh:297]: can only read Gmsh version 2 files Fair enough. I installed gmsh 3 and got past this error. Now I have another one that I can not understand. The error message, with some log before it, follows. Reading parameters from file mpfa2p.input. Reading 2d Gmsh grid... version 2.2 Gmsh file detected file contains 21574 nodes file contains 43150 elements number of real vertices = 21574 number of boundary elements = 500 number of elements = 42646 Initializing problem 'mpfal2p' Dune reported error: Dune::NotImplemented [storeInteractionVolumeInfo:.../dumux/dumux/porousmediumflow/2p/sequential/diffusion/mpfa/lmethod/2dpressure.hh:1156]: fvmpfao2pfaboundpressure2p.hh, l. 1164: boundary shape not available as interaction volume shape So I thought that I am misusing gmsh and tried one of the pre-installed .msh grids, namely dune-grid/doc/grids/gmsh/curved2d.msh. I got the same error. So probably I am misusing MPFA-L implementation and not gmsh. What should I do to have an unstructured grid in this kind of simulation? Thank you. 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
