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;
