Hi Lorenzo,
sorry, I didn't realize that you were using Dumux2.12. In this case,
forget my last mail. In the future, please try to reply to previous
mails of the conversation so that they can be read (I am sure you
mentioned somewhere that you were using 2.12).
In DuMuX2.12 there was no
Hi Lorenzo,
you have to distinguish between primary variable indices and equation
indices. Your equations (2p model) are always the mass balances of the
two phases, independent of you choice of primary variables. In
2p/implicit/incompressible/problem.hh you have an example on how to
adress
Hello,
On 02.11.2018 10:45, Dennis Gläser wrote:
per default, the formulation for the 2p model is pw-sn. That means
your primary variables are the water pressure and the non-wetting
phase saturation (in your case oil I assume). Therefore,
Indices::swIdx does not exists as it is not part of
On 02.11.18 12:07, lc wrote:
Hello,
I enclose all the files that I use.
Up to now in tutorialproblem_implicit.hh I introduce a double
small_number = 1e-6; which works ok when I initialize the oil
saturation profile (snIdx) but it is not effective on swIdx.
Hi Lorenzo,
I'm sorry, I
Hello,
I enclose all the files that I use.
Up to now in tutorialproblem_implicit.hh I introduce a double
small_number = 1e-6; which works ok when I initialize the oil saturation
profile (snIdx) but it is not effective on swIdx.
Also as I general question I'm wondering which is the correct
Hi Lorenzo,
it all depends on your setup. Probably something wrong with your custom
material law.
What have you tried so far to fix your issue?
Timo
> On 2. Nov 2018, at 11:53, lc wrote:
>
> Hello Dennis,
>
> thank you. I did as you suggested but as I do this, the Newton solver does
> not
Good morning Lorenzo,
per default, the formulation for the 2p model is pw-sn. That means your
primary variables are the water pressure and the non-wetting phase
saturation (in your case oil I assume). Therefore, Indices::swIdx does
not exists as it is not part of your primary variables.
Good morning,
I'd like to ask how to impose an initial condition on the water
saturation and not on the oil saturation.
Actually, what I implemented is the following:
const auto pos = fvGeometry.subContVol[scvIdx].global;
double x_0 = 0.;
double smin = 0.2121;