OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bounding Lmin and Lmax between 0. and 1.

Browse Source 
when solving Rachford Rice equations with solveRachfordRice_g_()
This commit is contained in:
Kai Bao 2023-11-01 22:04:53 +01:00
parent 3ff7fffb04
commit aacf4a4d81
1 changed files with 11 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/material/constraintsolvers/PTFlash.hpp
View File

@ -297,10 +297,8 @@ protected:
typename Vector::field_type Kmin = K[0];
typename Vector::field_type Kmax = K[0];
for (int compIdx=1; compIdx<numComponents; ++compIdx){
if (K[compIdx] < Kmin)
Kmin = K[compIdx];
else if (K[compIdx] >= Kmax)
Kmax = K[compIdx];
Kmin = std::min(Kmin, K[compIdx]);
Kmax = std::max(Kmin, K[compIdx]);
}
// Lower and upper bound for solution
@ -308,13 +306,20 @@ protected:
auto Lmax = Kmax / (Kmax - 1);
// Check if Lmin < Lmax, and switch if not
if (Lmin > Lmax)
{
if (Lmin > Lmax) {
auto Ltmp = Lmin;
Lmin = Lmax;
Lmax = Ltmp;
}
Lmin = std::clamp(Lmin, 0., 1.);
Lmax = std::clamp(Lmax, 0., 1.);
if (Lmin == Lmax) {
Lmin = 0.;
Lmax = 1.0;
}
// Initial guess
auto L = (Lmin + Lmax)/2;