Dear Lorenzo,
in my opinion it makes no sense to set a Dirichlet and a Neumann
condition at the same time. Either you know the flux over your boundary
or you know the values of your primary variables on the boundary. If you
change your code to setAllNeumann on the left boundary of your
Hi Lorenzo,
you still have large capillary pressure gradients at the left boundary
at t=600 thus the saturation changes.
You also have them in the initial conditions, so they probably never go
away. Why this is so, depends on your setup. Does it also occur when
your grid is fine from the
Dear Lorenzo,
generally for cell-centered discretizations in Dumux, Dirichlet boundary
conditions are weakly enforced on the boundary. Your first degree of
freedom is half a cell away from the boundary, so it doesn't necessarily
have to be the same as the Dirichlet value.
In order to see if
Hello Dumux team,
We are, finally, in the post-processing phase, aiming at publishing soon
... but we observed a scaring issue concerning the verification of the
BC at the inlet.
I'm using Dumux 2.12 with the 2p sequential model for a rectangular
domain discretized with a cartesian