Saturation values above 1

Open paveltomin opened this issue 2 years ago • 1 comments

GEOS formulation allows saturation values to go above 1 (not sure about below zero). The computed saturation values are used as is to compute rel perms that might be not really physical thing to do and also may mislead the solver blocking some flow or allowing too much flow. More reasonable way would be to normalize saturation values used for rel perm computes, see, e.g., https://drive.google.com/open?id=10dHD58pO8dy5hSz4RjVo8m0Mjs6C0a69

Sep 21 '23 15:09 paveltomin

Attaching small test case. The first nonlinear steps for this case show saturation changes of more than 1 indicating that (at least the intermediate) saturations are greater than 1. caseFiles.tar.gz

Jun 06 '24 01:06 dkachuma

Should be addressed by new formulation in https://github.com/GEOS-DEV/GEOS/pull/3443

Dec 05 '24 17:12 paveltomin