Dear Gergely,
sorry for the late response. Without having seen your code it is quite
hard to tell why your program stops.
Have you tried running valgrind?
Anyway, there might be a different solution for getting the fluxes
across the boundary (unfortunately, no standard method is implemented yet):
You can call
localResidual.evalFluxes(element, elemVolVars);
localResidual.evalStorage(element);
and then sum up the entries of the element residual
localResidual.residual() (in case you are using box).
If you have no object of localResidual at hand, you can simply
instantiate one and initialize it with the problem.
LocalResidual localResidual;
localResidual.init(problem);
This has to be repeated for all elements on the boundary of interest.
Best regards
Kilian
On 03/01/2017 03:26 PM, Gergely Schmidt wrote:
Dear Dumux community,
I’d like to calculate the mass flow through a boundary face for an
implicit 2p2c case for multiple time steps. Is there an already
implemented way? If not, is it possible to duplicate the
localresidual.hh/computeFlux function and extend its parameter list
with specific entities of the types Problem, Element,
FVElementGeometry and ElementVolumeVariables which can then be used to
update fluxVars? The new computeFlux function could then be called
with onBoundary=true, couldn’t it? But how can I proceed from there?
My compiled program stops without errors after the first call of
asImp_()->computeAdvectiveFlux(flux, fluxVars) inside the new
computeFlux function.
Many thanks and best regards,
Gergő Schmidt
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Kilian Weishaupt M.Sc.
Institut für Wasser- und Umweltsystemmodellierung (IWS)
Lehrstuhl für Hydromechanik und Hydrosystemmodellierung
Universität Stuttgart, Pfaffenwaldring 61, 70569 Stuttgart
Email: [email protected]
Telefon: 0049 711 685-60461 ** fax: 0049-711-685-60430
http://www.hydrosys.uni-stuttgart.de
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
