Dear Martin,
thank you very much, evaluating in half-cell distance worked.
Best regards,
Gergely
Am 16.02.2021 um 16:48 schrieb Martin Schneider:
Dear Gergely,
when using the Box method, this is not straightforward, because
calculating fluxes for faces at the boundary
will not result in the expected value, because Dirichlet conditions
are not weakly enforced.
Therefore, elemVolVars[scvf.insideScvIdx()] will return the volVars
which are calculated from the boundary values.
Since in your case Sn=0, the mobility will also be 0.
Instead, you could calculate fluxes over interior scvfs such that
scvf.outsideScvIdx() is related to your boundary DOF
but scvf.insideScvIdx() is not located at the boundary.
I hope this helps.
Best regards,
Martin
On 16.02.21 16:22, Gergely Schmidt wrote:
Dear Dumux community,
I'd like to calculate the the volumetric fluxes through a Dirichlet
boundary. Does somebody know, how I can access the upwind mobility of
an outflowing phase?
elemVolVars[scvf.insideScvIdx()].mobility(phaseIdx) gives me the
downwind value.
As an illustration, I can get the advective fluxes in a postTimeStep
function acc. to box/darcyslaw.hh for a porousmediumflow/2p test
problem, which looks like this:
Noflow BC
________________________
| |
Dirichlet BC for Sn = 0 | | Neumann BC <- gas
inflow
|________________________|
--> x Noflow BC
Left-side gas outflow must equal right inflow here, but the
calculated volumetric flux is always 0 m/s with
elemVolVars[scvf.insideScvIdx()].mobility(phaseIdx).
Thank you and best regards,
Gergely
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
