JIANG/EPANET
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix for PRV directly downstream of a PSV

Browse Source
This commit is contained in:
Lew Rossman
2022-08-03 19:52:01 -04:00
parent 44ad2488e4
commit 93389ad1e4
1 changed files with 22 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
src/hydcoeffs.c
View File

@@ -62,6 +62,7 @@ static void tcvcoeff(Project *pr, int k);
static void prvcoeff(Project *pr, int k, int n1, int n2); static void prvcoeff(Project *pr, int k, int n1, int n2);
static void psvcoeff(Project *pr, int k, int n1, int n2); static void psvcoeff(Project *pr, int k, int n1, int n2);
static void fcvcoeff(Project *pr, int k, int n1, int n2); static void fcvcoeff(Project *pr, int k, int n1, int n2);
static void valveconnectcoeffs(Project *pr, int k, int n1, int n2);
void resistcoeff(Project *pr, int k) void resistcoeff(Project *pr, int k)
@@ -317,9 +318,11 @@ void valvecoeffs(Project *pr)
{ {
case PRV: case PRV:
prvcoeff(pr, k, n1, n2); prvcoeff(pr, k, n1, n2);
valveconnectcoeffs(pr, k, n1, n2);
break; break;
case PSV: case PSV:
psvcoeff(pr, k, n1, n2); psvcoeff(pr, k, n1, n2);
valveconnectcoeffs(pr, k, n1, n2);
break; break;
case FCV: case FCV:
fcvcoeff(pr, k, n1, n2); fcvcoeff(pr, k, n1, n2);
@@ -329,6 +332,25 @@ void valvecoeffs(Project *pr)
} }
} }
void valveconnectcoeffs(Project *pr, int k, int n1, int n2)
/*
**--------------------------------------------------------------
** Input: k = valve's link index
** n1 = upstream node of valve
** n2 = downstream node of valve
** Output: none
** Purpose: insures that the off-diagonals in the rows for a
** PRV/PSV are non-zero to preserve connectivity.
**--------------------------------------------------------------
*/
{
Smatrix *sm = &pr->hydraul.smatrix;
double p = 1.0 / CBIG;
sm->Aij[sm->Ndx[k]] -= p;
sm->Aii[sm->Row[n1]] += p;
sm->Aii[sm->Row[n2]] += p;
}
void emittercoeffs(Project *pr) void emittercoeffs(Project *pr)
/* /*