Dear Mahmoud,
as far as I know, this is not possible out-of-the-box. If you generate
the mesh externally from a geometry, you may try duplicating the
interface line (I am assuming 2d here) between the two layers so that
they get meshed "twice". Then, these two lines end up as boundaries in
the computational domain, at which you can prescribe no-flow boundary
conditions. Gmsh, for instance, also has a feature called "crack" (or
similar), where you can split the mesh along defined lines after meshing.
If you do such thing, it only works if there are no floating
(sub-)domains around that are not connected somehow to Dirichlet
boundary conditions. I am assuming your interface does not go through
the entire domain, as otherwise you'd have two decoupled systems. In any
case, I do not see your attached image with your setup. My guess is that
the mailing list does not allow for attachments?
Cheers,
Dennis
On 28.04.21 05:52, Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151 wrote:
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
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
