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

Removed unused code.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-03-05 12:50:13 +01:00
parent c3d6cc429a
commit f783f9eebe
1 changed files with 0 additions and 30 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
opm/polymer/TransportModelPolymer.cpp
View File

@ -289,31 +289,6 @@ namespace Opm
}
}
void computeExplicitStep(const double* xmin, const double* xmax, double* x) const {
double ff = tm.fracFlow(s0, c0, cell);
double mc = tm.computeMc(c0);
double dps = tm.polyprops_.deadPoreVol();
//In this explicit step, we do not compute absorption and take ads0=ads
// double rhor = tm.polyprops_.rockDensity();
// double ads0 = tm.polyprops_.adsorbtion(std::max(c0, cmax0));
// double ads = tm.polyprops_.adsorbtion(std::max(c, cmax0));
x[0] = s0 - dtpv*(outflux*ff + influx);
x[1] = 1./(x[0] - dps)*((s0 - dps)*c0 - dtpv*(outflux*ff*mc + influx_polymer)); // + rhor*((1.0 - porosity)/porosity)*(ads - ads0)
// We check that the values we obtain remains admissible (this is not guaranted for an explicit step)
if (x[0] < xmin[0]) {
x[0] = xmin[0];
} else if (x[0] > xmax[0]) {
x[0] = xmax[0];
}
if (x[1] < xmin[1]) {
x[1] = xmin[1];
} else if (x[1] > xmax[1]) {
x[1] = xmax[1];
}
}
void computeResidual(const double* x, double* res) const
{
double s = x[0];
@ -653,11 +628,6 @@ namespace Opm
double x_new[2];
double res_new[2];
// // Update x=(s, c) with an explicit solver.
// residual.computeExplicitStep(x_min, x_max, x);
// residual.computeResidual(x, res);
while ((norm(res) > tol) && (iters_used_split < max_iters_split)) {
// We first try a Newton step
if (residual.solveNewtonStep(x, x_new, gradient_method)) {