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

Added TransportModelPolymer::setPreferredMethod().

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2012-06-13 15:44:46 +02:00
parent b6001c8dac
commit 42e862e7ec
2 changed files with 18 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
opm/polymer/TransportModelPolymer.cpp
View File

@@ -214,6 +214,14 @@ namespace Opm
void TransportModelPolymer::setPreferredMethod(SingleCellMethod method)
{
method_ = method;
}
void TransportModelPolymer::solve(const double* darcyflux, void TransportModelPolymer::solve(const double* darcyflux,
const double* porevolume, const double* porevolume,
const double* source, const double* source,

16
opm/polymer/TransportModelPolymer.hpp
View File

@@ -50,6 +50,9 @@ namespace Opm
const double tol, const double tol,
const int maxit); const int maxit);
/// \TODO document me.
void setPreferredMethod(SingleCellMethod method);
/// Solve transport eqn with implicit Euler scheme, reordered. /// Solve transport eqn with implicit Euler scheme, reordered.
/// \TODO Now saturation is expected to be one sw value per cell, /// \TODO Now saturation is expected to be one sw value per cell,
/// change to [sw so] per cell. /// change to [sw so] per cell.
@@ -62,6 +65,13 @@ namespace Opm
std::vector<double>& concentration, std::vector<double>& concentration,
std::vector<double>& cmax); std::vector<double>& cmax);
void solveGravity(const std::vector<std::vector<int> >& columns,
const double* porevolume,
const double dt,
std::vector<double>& saturation,
std::vector<double>& concentration,
std::vector<double>& cmax);
public: // But should be made private...
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);
void solveSingleCellBracketing(int cell); void solveSingleCellBracketing(int cell);
@@ -73,12 +83,6 @@ namespace Opm
const int pos, const int pos,
const double* gravflux); const double* gravflux);
int solveGravityColumn(const std::vector<int>& cells); int solveGravityColumn(const std::vector<int>& cells);
void solveGravity(const std::vector<std::vector<int> >& columns,
const double* porevolume,
const double dt,
std::vector<double>& saturation,
std::vector<double>& concentration,
std::vector<double>& cmax);
private: private: