Hi Nikolai,
you can have a look at the fracture exercise in the DuMuX course: https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-course/tree/master/exercises/exercise-fractures Therein, we use the Gmsh domain markers to distribute fracture apertures. The explanation would be a bit lengthy and I guess you can find all the information you need within the exercise description. The main concept is to use the GridData object from the GridManager. In the main file of the exercise you can see how this data is retrieved and then passed to the sub-problems and spatial parameters classes. In short words, in Gmsh you can assign an index (domain marker) to each fracture line (2d domain) or fracture surface (3d domain). All elements that lie on the line/surface will then get this domain marker in the mesh file and it is possible to retrieve this again in DuMuX. Rergarding you question of the interface. As for the other parameter interface functions, for the extrusion factor there is also the version with the suffix AtPos(), but there is also the version taking element/scv/elemsol (see fvproblem.hh): template<class ElementSolution> Scalar extrusionFactor(const Element& element, const SubControlVolume& scv, const ElementSolution& elemSol) const {} Note that the element index alone is not enough to distribute your parameters because depending on the grid manager (ALUGrid, UG, ...), the indices might not be equal to the order of the .msh file. You would need the insertion index of the element. The way to go are the domain markers I would say. I hope this helps! Best wishes, Dennis ________________________________ Von: Dumux <[email protected]> im Auftrag von Nikolai Andrianov <[email protected]> Gesendet: Mittwoch, 12. Dezember 2018 18:40:21 An: DuMuX User Forum Betreff: [DuMuX] Variable fracture aperture Dear DuMuX experts, Please advise what is the best way to assign variable fracture aperture for a multidomain/fracture flow application. I would like to have apertures be pre-computed and read by DuMuX from an external file, so that the fracture apertures correspond to the numbering of fracture elements in the gmsh mesh file. Currently the apertures are read via extrusionAtPos(const GlobalPosition& globalPos) in the fractureproblem.hh. It would be great if I could get not the position of the element’s center as parameter, but say the corresponding SubControlVolume which apparently does have the elementIndex. But how can I do this without modifying dumux/common/fvproblem.hh ? Thanks, Nikolai
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
