OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-15 23:01:54 -06:00
Code Issues Packages Projects Releases Wiki Activity

Simplified computation of adsorbed mass of polymer.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-09-03 12:49:04 +02:00
parent a335e1b149
commit 708c56bb94
1 changed files with 1 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/polymer/polymerUtilities.cpp
View File

@ -219,16 +219,14 @@ namespace Opm
{
const int num_cells = props.numCells();
const double rhor = polyprops.rockDensity();
std::vector<double> porevolume;
std::vector<double> porosity;
computePorevolume(grid, props.porosity(), rock_comp, state.pressure(), porevolume);
computePorosity(grid, props.porosity(), rock_comp, state.pressure(), porosity);
double abs_mass = 0.0;
const std::vector<double>& cmax = state.maxconcentration();
for (int cell = 0; cell < num_cells; ++cell) {
double c_ads;
polyprops.simpleAdsorption(cmax[cell], c_ads);
abs_mass += c_ads*porevolume[cell]*((1.0 - porosity[cell])/porosity[cell])*rhor;
abs_mass += c_ads*grid.cell_volumes[cell]*(1.0 - porosity[cell])*rhor;
}
return abs_mass;
}