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

Reply via email to