Hi Anwar,
you read in a scalar permeabilty value into the variable
`permeability_`, but your function `permeabilityAtPos` returns a
`DimWorldMatrix`. That means that this scalar value will be in each
entry of the full permeability tensor. You should be able to fix this by
changing line 66 of your spatialparams file to
using PermeabilityType = Scalar;
Kind regards
Bernd
On 08/07/2018 10:08 AM, Anwar Al Assadi wrote:
Hello everybody,
I have a very strange unusual flow distribution after injecting a gas
inside my domain (see attachment).
First of all I have a rectangular domain. My boundary conditions are
on the left and right side: Dirichlet ( with temperature gradient and
hydrostatic pressure). The bottom and top side are considered as
Neumann, no flow. My injection is in the middle of horizontal side.
During the injection time everything seems okay, the distribution is
symmetric. After the injection my gas starts to migrate to left side.
It should actually migrate to the top side. I tried that as non
isothermal as well, but I'm still getting the same results.
I will be very thankful, if you offer my any suggestions to stop the
migration to left side.
facts:
model: 2pnc /2pncni
fluidsystem: brineco2 (from stabel with hardcoded salinity)
dumux tag: dumux-course-2018
greetings
Anwar
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
_______________________________________________________________
Bernd Flemisch phone: +49 711 685 69162
IWS, Universität Stuttgart fax: +49 711 685 60430
Pfaffenwaldring 61 email: [email protected]
D-70569 Stuttgart url: www.hydrosys.uni-stuttgart.de
_______________________________________________________________
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
