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
--
***********************************************************
Gergely Schmidt M.Sc.
Wissenschaftl. Mitarbeiter
Institut für Strömungsmechanik und Umweltphysik im Bauwesen
Leibniz Universität Hannover
Appelstr. 9 a
30167 Hannover
tel +49 (0) 511 / 762 4205
***********************************************************
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux