Dear experts,
I would like to calculate the flux across the fracture-matrix interfaces, but
get a runtime error when using the code below (a part of my matrixproblem.hh):
FluxVariables fluxVars;
fluxVars.init(*this, element, fvGeometry, elemVolVars, scvf,
elemFluxVarsCache);
if (couplingManager().isOnInteriorBoundary(element, scvf))
{
Scalar flux =
fluxVars.advectiveFlux(FluidSystem::phase0Idx, upwindTerm); // <-- Runtime
error here
std::cout << ", flux=" << flux << std::endl;
mfFlux += flux;
}
The code successfully prints the flux for a first scvf but gives the following
error message on a second scvf: "Assertion `std::find(s.begin(), s.end(),
lowDimElemIdx) != s.end()' failed.'
The error message occurs in couplingManager().getLowDimVolVars(). It seems that
the low-dimensional element can be found, but the corresponding volume
variables cannot be linked.
Appreciate your help 😊
Thanks,
Nikolai
PS I am (still) using Dumux 2.12..
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
