OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

update PAV formulation.

Browse Source
This commit is contained in:
Liu Ming 2016-09-05 10:38:57 +08:00
parent c5a490c75a
commit 1943978dbd
1 changed files with 16 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
opm/polymer/fullyimplicit/FullyImplicitCompressiblePolymerSolver.cpp
View File

@ -564,16 +564,26 @@ namespace {
const int dims = *std::max_element(fipnum.begin(), fipnum.end()); const int dims = *std::max_element(fipnum.begin(), fipnum.end());
std::vector<V> values(dims, V::Zero(5)); std::vector<V> values(dims, V::Zero(7));
for (int d = 0; d < dims; ++d) {
for (int c = 0; c < nc; ++c) {
for (int i = 0; i < 5; ++i) { for (int i = 0; i < 5; ++i) {
if (fipnum[c] == d) { for (int c = 0; c < nc; ++c) {
values[d][i] += fip[c][i]; if (fipnum[c] != 0) {
values[fipnum[c]-1][i] += fip[i][c];
} }
} }
} }
// compute PAV and PORV or every regions.
for (int c = 0; c < nc; ++c) {
if (fipnum[c] != 0) {
values[fipnum[c]-1][5] += pv[c] * s.col(Oil)[c];
values[fipnum[c]-1][6] += pv[c] * state.pressure.value()[c] * s.col(Oil)[c];
}
}
for (auto& val : values) {
val[6] = val[6] / val[5];
} }
return values; return values;