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

Use Value<double> for Connection Factor

Browse Source
This commit is contained in:
Joakim Hove 2015-02-03 13:10:38 +01:00
parent 034961c8ba
commit e435545eee
1 changed files with 21 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
opm/core/wells/WellsManager_impl.hpp
View File

@ -166,24 +166,28 @@ void WellsManager::createWellsFromSpecs(std::vector<WellConstPtr>& wells, size_t
int cell = cgit->second; int cell = cgit->second;
PerfData pd; PerfData pd;
pd.cell = cell; pd.cell = cell;
if (completion->getConnectionTransmissibilityFactor() > 0.0) { {
pd.well_index = completion->getConnectionTransmissibilityFactor(); const Value<double>& transmissibilityFactor = completion->getConnectionTransmissibilityFactorAsValueObject();
} else {
double radius = 0.5*completion->getDiameter(); if (transmissibilityFactor.hasValue()) {
if (radius <= 0.0) { pd.well_index = transmissibilityFactor.getValue();
radius = 0.5*unit::feet; } else {
OPM_MESSAGE("**** Warning: Well bore internal radius set to " << radius); double radius = 0.5*completion->getDiameter();
if (radius <= 0.0) {
radius = 0.5*unit::feet;
OPM_MESSAGE("**** Warning: Well bore internal radius set to " << radius);
}
const std::array<double, 3> cubical =
WellsManagerDetail::getCubeDim<3>(c2f, begin_face_centroids, cell);
const double* cell_perm = &permeability[dimensions*dimensions*cell];
pd.well_index =
WellsManagerDetail::computeWellIndex(radius, cubical, cell_perm,
completion->getSkinFactor(),
completion->getDirection(),
ntg[cell]);
} }
const std::array<double, 3> cubical =
WellsManagerDetail::getCubeDim<3>(c2f, begin_face_centroids, cell);
const double* cell_perm = &permeability[dimensions*dimensions*cell];
pd.well_index =
WellsManagerDetail::computeWellIndex(radius, cubical, cell_perm,
completion->getSkinFactor(),
completion->getDirection(),
ntg[cell]);
} }
wellperf_data[well_index].push_back(pd); wellperf_data[well_index].push_back(pd);
} else { } else {