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;
double smax = 0.7856;
double exp_arg = (0.5 -
pow((pos[0]-x_0)/(2.+cos(3.*sqrt(std::abs(pos[1]-75.0)))), 2.) );
double exp_value = exp_arg >= -100 ? exp(exp_arg) : 0;
double small_number = 1e-6;
* values[Indices::snIdx] = std::min(smax, smin + exp_value +
small_number);*
but this last line, should be imposed on values[Indices::swIdx] but if I
write it I get an error, not declared variable.
Thank you,
Lorenzo
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
