Hello,
I am a bit confused about my meshgrid choice because I get completely
different results depending on the small changes I may do.
I am injecting a gas (using the addPointSources method) into a fully
saturated 2D cell. The cell is mostly filled with coarse sand + two lenses
of a lower permeability sand:
[image: image.png]
Now with:
#The grid type
SET_TYPE_PROP(HeterogeneousProblem, Grid, Dune::ALUGrid<2, 2, Dune::*cube*,
Dune::nonconforming>);
and
const MaterialLawParams& materialLawParams(const Element &element,
const FVElementGeometry
&fvGeometry,
const int scvIdx) const
{
* const GlobalPosition &globalPos =
fvGeometry.subContVol[scvIdx].global;*
I get this:
[image: image.png]
- When I change *cube* to *simplex* I get this:
[image: image.png]
the gas shifts to the left ! and penetrate the low permeability lenses !!
- When I change *fvGeometry.subContVol[scvIdx].global;* to
*element.geometry().center();* I get:
[image: image.png]
PS: my BC and gas injection rate didn't change (same for the three cases
above).
Based on lab results, we see a vertical migration of the gas into the cell
with widening of the plume (which I get with a *simplex* mesh!)
However, when the gas reaches the lens of the finer sand, it accumulates at
the interface between the lens and the coarse sand and migrates laterally
to the left then upwards (which I get with the *cubic* mesh)!
I am confused about how the mesh influences the transport of the gas inside
the cell and why using *element.geometry().center() *forces the gas to
penetrate the low permeability lens?
Thank you,
Kenza
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
