Dear All,
I am trying to implement multidoamin/poromechanics/el2p for heterogeneous
problem.
Considering the value of porosity and permeability are stored in vector
initPoro_ and initK_ respectively. Porosity and permeability is returned
using element Idx in their respective methods.
//! Returns the porosity of the porous medium.
template<class ElementSolution>
Scalar porosity(const Element& element,
const SubControlVolume& scv,
const ElementSolution& elemSol) const
{
const auto eIdx =
this->fvGridGeometry().elementMapper().index(element);
return
PorosityDeformation<Scalar>::evaluatePorosity(this->fvGridGeometry(),
element, scv, elemSol, initPoro_[eIdx]);
}
The porosity output for 2p and poroelastic are different (can be viewed
using below link):
https://home.iitb.ac.in/index.php/s/M7b0ZpTNFt5Wesr (2p)
https://home.iitb.ac.in/index.php/s/uJFpUoq6rtxQo25 (el)
Could you please give some idea the fix the above issue?
Thank You,
Regards,
Ranjeet
_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
