Re: [DuMux] Question on sourceAtPos function

2021-04-09 Thread Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151
Dear Timo, Really Thank you very very much. You clarified everything for me in an amazing way and finally the results are going in the correct way now and they make sense despite some numerical error modifying the pressure on the upper and lower boundaries where the injection well is

Re: [DuMux] Question on sourceAtPos function

2021-04-09 Thread Timo Koch
Hi Mahmoud, you also need to extract the correct energy carried by the fluid in the production well. Otherwise your domain will heat up there. For solution-dependent point sources you can use the split interface. First you add your point sources (in addPointSources) only providing positions as

Re: [DuMux] Question on sourceAtPos function

2021-04-09 Thread Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151
Dear Timo, I cannot thank you enough for all your help and guidance. I wrote the pointsource function as shown below. could you please let me know if it is OK now or not ? I used 3 pointsources to represent each of the production well and the injection well because I wanted to represent

Re: [DuMux] Question on sourceAtPos function

2021-04-08 Thread Timo Koch
> On 8. Apr 2021, at 13:42, Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151 > wrote: > > Dear Timo, > > Thanks a lot for your reply. I didn't use the pointSources because I didn't > know how to assign the temperature inside such function for the injection > well because I want to impose a

Re: [DuMux] Question on sourceAtPos function

2021-04-08 Thread Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151
Dear Timo, Thanks a lot for your reply. I didn't use the pointSources because I didn't know how to assign the temperature inside such function for the injection well because I want to impose a rate and also a temperature. I'm using the OnePNI model in my problem. I will try to be more

Re: [DuMux] Question on sourceAtPos function

2021-04-07 Thread Timo Koch
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