Hi, I am starting to get familiarized with the solution of one-phase flow in fractured media in Dumux 3.2. I started by trying to simulate the first benchmark proposed in Flemisch et al (2018) (reference below) witht he Box and the Ctpfa methods. The paper used a preivous version of Dumux, when the multidomain coupling manager did not exist.
It is a 2d domain containning a regular fracture network (gmsh file is attached). I used an example in Dumux-pub and basically only changed the dimension of the domain from 2d to 3d and altered the boundary conditions (project is attached). First of all, I have a question about the multidomain assembler. It deals with the two domains by solving the fracture and the matrix domains separetely, treating the jacobian and solution vectors as block structures, which makes this stationary problem non-linear. This led me to convergence issues, since the Box method does not converge for any of the meshes I tested. *Is there the possibility of handling this problem in a fully coupled way ?* that is, making the multi-domain assemblage so I have only one matrix where the nodes which contain both fracture and matrix have the contributions of both domains added up. This way I would not need a non-linear solver. Any function to do that while maintaining the multi-domain context ? If not, I would like to hear some insights, if you have any, on what may be causing the non-convergence of the bok method. I use the UMFPackBackend linear solver and I have tried many other finer and coarser meshes. Reference: Flemisch, B., Berre, I., Boon, W., Fumagalli, A., Schwenck, N., Scotti, A., Stefansson, I., Tatomir, A., 2018. Benchmarks for single-phase flow in fractured porous media. Advances in Water Resources. 111, 239–258. Thank you. Best regards, Ana
benchmark1_project.rar
Description: Binary data
_______________________________________________ DuMux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
