Hi Kenza, I’m sorry that we haven’t come up with an answer to your email until now. Before taking a closer look: does the problem still persist in the way that you describe it?
Kind regards Bernd Von: DuMux <[email protected]> Im Auftrag von kenza bouznari Gesendet: Freitag, 12. Februar 2021 17:34 An: [email protected] Betreff: [DuMux] Meshgrid influence on results 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.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.png] - When I change cube to simplex I get this: [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.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
