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

Use ADB for solvent fraction

Browse Source
This commit is contained in:
Tor Harald Sandve
2015-09-11 14:17:30 +02:00
parent 64c5b29b87
commit 853e1f0dc1
1 changed files with 9 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
opm/autodiff/BlackoilSolventModel_impl.hpp
View File

@@ -250,17 +250,19 @@ namespace Opm {
? state.saturation[ pu.phase_pos[ Gas ] ]
: zero);
Selector<double> zero_selector(ss.value(), Selector<double>::Zero);
V F_solvent = zero_selector.select(ss, ss / (ss + sg)).value();
const std::vector<int> well_cells(wells().well_cells, wells().well_cells + nperf);
const int nw = wells().number_of_wells;
V wellSolventFraction = Eigen::Map<const V>(&xw.solventFraction()[0], nperf);
Selector<double> zero_selector(ss.value() + sg.value(), Selector<double>::Zero);
ADB F_solvent = subset(zero_selector.select(ss, ss / (ss + sg)),well_cells);
const int nw = wells().number_of_wells;
V injectedSolventFraction = Eigen::Map<const V>(&xw.solventFraction()[0], nperf);
V isProducer = V::Zero(nperf);
V ones = V::Constant(nperf,1.0);
for (int w = 0; w < nw; ++w) {
if(wells().type[w] == PRODUCER) {
for (int perf = wells().well_connpos[w]; perf < wells().well_connpos[w+1]; ++perf) {
wellSolventFraction[perf] = F_solvent[well_cells[perf]];
isProducer[perf] = 1;
}
}
}
@@ -271,7 +273,7 @@ namespace Opm {
// remove contribution from the dissolved gas.
// TODO compensate for gas in the oil phase
assert(!has_vapoil_);
const ADB cq_s_solvent = wellSolventFraction * (cq_s[gas_pos] - rs_perfcells * cq_s[oil_pos]);
const ADB cq_s_solvent = (isProducer * F_solvent + (ones - isProducer) * injectedSolventFraction) * (cq_s[gas_pos] - rs_perfcells * cq_s[oil_pos]);
// Solvent contribution to the mass balance equation is given as a fraction
// of the gas contribution.