Hi Dennis, Thank you for your reply. You got it correctly, difference in plot due to vertex-wise output of an element-wise field. As I haven't implemented any code to write porosity to vtk file, I think, porosity() should returns porosity with respect to vertex instead of cell. if this is true, how one can construct the value at vertex from the cell values?.
Regards, Ranjeet On Thu, Feb 14, 2019 at 3:54 AM Dennis Gläser < [email protected]> wrote: > 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 > [email protected]https://listserv.uni-stuttgart.de/mailman/listinfo/dumux > >
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
