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

Reordering 2p transport solver now takes full saturation (all phases) as input.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-30 14:39:09 +02:00
parent ad2dff642b
commit 299ea28afe
2 changed files with 10 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
opm/core/transport/reorder/TransportModelTwophase.cpp
View File

@ -51,7 +51,7 @@ namespace Opm
darcyflux_(0), darcyflux_(0),
source_(0), source_(0),
dt_(0.0), dt_(0.0),
saturation_(0), saturation_(grid.number_of_cells, -1.0),
fractionalflow_(grid.number_of_cells, -1.0), fractionalflow_(grid.number_of_cells, -1.0),
mob_(2*grid.number_of_cells, -1.0) mob_(2*grid.number_of_cells, -1.0)
#ifdef EXPERIMENT_GAUSS_SEIDEL #ifdef EXPERIMENT_GAUSS_SEIDEL
@ -79,13 +79,13 @@ namespace Opm
const double* porevolume, const double* porevolume,
const double* source, const double* source,
const double dt, const double dt,
double* saturation) std::vector<double>& saturation)
{ {
darcyflux_ = darcyflux; darcyflux_ = darcyflux;
porevolume_ = porevolume; porevolume_ = porevolume;
source_ = source; source_ = source;
dt_ = dt; dt_ = dt;
saturation_ = saturation; toWaterSat(saturation, saturation_);
#ifdef EXPERIMENT_GAUSS_SEIDEL #ifdef EXPERIMENT_GAUSS_SEIDEL
std::vector<int> seq(grid_.number_of_cells); std::vector<int> seq(grid_.number_of_cells);
@ -103,6 +103,7 @@ namespace Opm
#endif #endif
reorderAndTransport(grid_, darcyflux); reorderAndTransport(grid_, darcyflux);
toBothSat(saturation_, saturation);
} }
// Residual function r(s) for a single-cell implicit Euler transport // Residual function r(s) for a single-cell implicit Euler transport
@ -617,9 +618,7 @@ namespace Opm
cells[c] = c; cells[c] = c;
} }
mob_.resize(2*nc); mob_.resize(2*nc);
std::vector<double> boths; props_.relperm(cells.size(), &saturation[0], &cells[0], &mob_[0], 0);
Opm::toBothSat(saturation, boths);
props_.relperm(cells.size(), &boths[0], &cells[0], &mob_[0], 0);
const double* mu = props_.viscosity(); const double* mu = props_.viscosity();
for (int c = 0; c < nc; ++c) { for (int c = 0; c < nc; ++c) {
mob_[2*c] /= mu[0]; mob_[2*c] /= mu[0];
@ -629,7 +628,7 @@ namespace Opm
// Set up other variables. // Set up other variables.
porevolume_ = porevolume; porevolume_ = porevolume;
dt_ = dt; dt_ = dt;
saturation_ = &saturation[0]; toWaterSat(saturation, saturation_);
// Solve on all columns. // Solve on all columns.
int num_iters = 0; int num_iters = 0;
@ -639,6 +638,8 @@ namespace Opm
} }
std::cout << "Gauss-Seidel column solver average iterations: " std::cout << "Gauss-Seidel column solver average iterations: "
<< double(num_iters)/double(columns.size()) << std::endl; << double(num_iters)/double(columns.size()) << std::endl;
toBothSat(saturation_, saturation);
} }
} // namespace Opm } // namespace Opm

4
opm/core/transport/reorder/TransportModelTwophase.hpp
View File

@ -43,7 +43,7 @@ namespace Opm
const double* porevolume, const double* porevolume,
const double* source, const double* source,
const double dt, const double dt,
double* saturation); std::vector<double>& saturation);
virtual void solveSingleCell(const int cell); virtual void solveSingleCell(const int cell);
virtual void solveMultiCell(const int num_cells, const int* cells); virtual void solveMultiCell(const int num_cells, const int* cells);
@ -71,7 +71,7 @@ namespace Opm
const double* porevolume_; // one volume per cell const double* porevolume_; // one volume per cell
const double* source_; // one source per cell const double* source_; // one source per cell
double dt_; double dt_;
double* saturation_; // one per cell std::vector<double> saturation_; // one per cell, only water saturation!
std::vector<double> fractionalflow_; // = m[0]/(m[0] + m[1]) per cell std::vector<double> fractionalflow_; // = m[0]/(m[0] + m[1]) per cell
// For gravity segregation. // For gravity segregation.
std::vector<double> gravflux_; std::vector<double> gravflux_;