Hi Ranjeet,
I don't get the "issue" you mentioned. What exactly is the problem? To
me your implementation seems fine.
If you refer to the porosity output you obtained in the el-problem, it
just looks like this is because you are using vertex-wise output of an
element-wise field.
Best wishes,
Dennis
On 13.02.19 05:00, Ranjeet kumar wrote:
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
_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux