Hello,
I am doing a porous medium flow simulation and I would like to specify
custom laws for relative permeability and capillary pressure that depend
on water saturation, surfactant concentration and pressure gradient.
I have my own spatial parameters class, of course, and I gather that the
following method should help me in what I am trying to achieve:
template<class ElementSolution>
MaterialLawParams materialLawParams(const Element& element,
const SubControlVolume& scv,
const ElementSolution& elemSol)
const
So there is no problem in getting water saturation and surfactant
concentration for the element, and also the pressure. However, I am
having trouble estimating the pressure gradient. To do that, I need to
know (in this method) the pressure in neighbor elements. But how do I
access the current solution values for those elements, while given only
an ElementSolution for the element of interest?
Best regards,
Dmitry
_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
