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

Added explicitTimestepLimit() method.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2011-04-11 13:37:30 +02:00
parent 58b7db31c5
commit 389d7627f4
1 changed files with 17 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
src/TPFACompressiblePressureSolver.hpp
View File

@ -341,6 +341,23 @@ public:
/// @brief
/// Explicit IMPES time step limit.
double explicitTimestepLimit(const std::vector<double>& totcompr,
const std::vector<double>& voldiscr,
const std::vector<double>& cellA, // num phases^2 * num cells, fortran ordering!
const std::vector<double>& faceA, // num phases^2 * num faces, fortran ordering!
const std::vector<double>& phasemobf,
const std::vector<double>& phasemobf_deriv,
const double* surf_dens)
{
compr_quantities cq = { 3, &totcompr[0], &voldiscr[0], &cellA[0], &faceA[0], &phasemobf[0] };
return cfs_tpfa_impes_maxtime(grid_.c_grid(), &cq, &trans_[0], &porevol_[0], data_,
&phasemobf_deriv[0], surf_dens, gravity_);
}
/// @brief
/// Explicit IMPES transport.
void explicitTransport(const double dt,