I see that dirichletAtPos() is available for Box method, too, but in
fact, I must use neumann() because by boundary conditions are more
complex than just the constant pressure. I have a two-phase flow and
need to implement a certain equation on the outflux of different phases,
depending on the pressure gradient.
Best regards,
Dmitry
On 26.05.2020 14:08, Dmitry Pavlov wrote:
Hello,
After transition from the CC method to the Box method in my 1d porous
medium flow problem, I am wondering how I should specify a constant
pressure boundary condition.
With the CC method, it was relatively easy: I could either implement
dirichletAtPos(), or implement neumann() instead and calculate the
flux across the boundary from the pressure gradient as
gradp = ([boundary pressure] - [pressure in the center of the cell]) /
(half-width of the cell)
With the Box method, though, it is no longer obvious to me. I will
appreciate any pointers. I see addCCDirichletFluxDerivatives in the
DuMux code, but no similar method for Box method.
Best regards,
Dmitry
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
