Hi Ranjeet,
the porosity is written out for the poro-elastic sub-problem, which is
based on the box method, i.e. the degrees of freedom are distributed on
the grid vertices. That is why the output looks like that.
If you open the .vtu files generated for the 2p problem, you should be
able to see the cell-wise porosities, unless you changed the output fields.
As a reference, you can run the test in
test/multidomain/poromechanics/el2p and open the file
test_md_poromechanics_el2p_twop.pvd
with paraview. There you can see the cell-wise porosities.
I hope this helps!
Best wishes,
Dennis
On 14.02.19 07:23, Ranjeet kumar wrote:
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]
<mailto:[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 list
[email protected] <mailto:[email protected]>
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
