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