Dear Mahmoud, what you are doing seems correct to me. You can get the relevant cell volume with “scv.volume()” so you don’t have to compute that yourself. You can also get the density with "elemVolVars[scv].density(phaseIdx)". But that shouldn’t change anything if you already used the correct volume in your computation.
You can also use pointSources (see e.g https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/test/porousmediumflow/1p/pointsources/timeindependent/problem.hh <https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/test/porousmediumflow/1p/pointsources/timeindependent/problem.hh>) which is simpler for the case you describe. There you just specify the position of the well(s) and the rate(s) (this time directly in kg/s without the need for considering cell volumes (that will be done automatically)). I guess you have to recheck your scenario to see if you would expect an influence on pressure for the specified production rate. Best wishes Timo > On 6. Apr 2021, at 18:05, Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151 > <[email protected]> wrote: > > Hello Dumux team, > > > I'm defining two wells inside a structured 3D Yasp grid, one is injection and > the other is production and I'm using the MPFA discretization. The system is > a geothermal system for which I'm defining a certain low temperature for > injection. I'm defining both wells inside the sourceAtPos method and they are > both defined inside 3 cells in the z direction. each cell is 10x10x10 m^3 Now > I understood from the doxygen documentation that the unit inside that source > function is Kg/(m^3.s) and I want to impose a rate of 100 m^3/day == 1.157e-3 > m^3/s in each well so I'm dividing by (3 x 1000) and then multiplying by 1000 > kg/m^3 for water to get an estimation of the mass rate to use so I get a > value around 4e-4 Kg/(m^3.s) to use. However whenever I use that value, the > simulation proceeds really quickly like I'm not injecting or producing > anything at all (static condition) and whenever I use a random higher mass > rate like 0.1 for example, I get very high pressures and warnings during the > simulation and the simulator aborts. Could somebody please let me know how to > correctly account for the 100 m^3/day that I want to impose in each of the > production and injection wells ? tried using some middle values and the > simulation proceeded successfully but I don't know the basis. Your guidance > is much appreciated. > > Best Regards, > > Mahmoud Atef > _______________________________________________ > DuMux mailing list > [email protected] > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________ DuMux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
