Hello Dumux team,
Is there a way in Dumux to impose zero transmissibility along an interface between 2 layers to have zero flux across that interface like in the industrial simulator Eclipse ? I'm simulating a geothermal aquifer with an upper and lower boundary layers where into the aquifer, there's an injection well (cooled water with low temperature) and a production well. I'm using the OnePNI model and CCMpfa discretization scheme. I want to monitor the heat transfer by conduction at the boundary layers and thus I need to have them inside my domain. I tried setting zero permeability for both upper and lower boundary layers but it doesn't help the numerics and the convergence becomes extremely difficult. So I increased the permeability of those boundary layers keeping it 1e^3 lower than that of the aquifer and the code became stable regarding the easiness of convergence but I could still see an anomaly in the pressure in those upper and lower boundary layers as shown in the attached image in which 3 numerical layers were used to represent each of the upper and lower boundary layers.
So Is there a way in dumux to ensure no flux across the interface between 2 layers ? Thanks in advance.
Best Regards, Mahmoud _______________________________________________ DuMux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
