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

Reply via email to