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

Deleting some spaces.

Browse Source 
FullyImplicitBlackoilSolver_impl.hpp
This commit is contained in:
Kai Bao 2014-05-07 14:54:44 +02:00
parent 930beabd21
commit ccc1577f09
1 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@ -1744,23 +1744,23 @@ namespace {
ADB tempBW = fluidReciprocFVF(Water, press, rs, rv, cond, cells_); ADB tempBW = fluidReciprocFVF(Water, press, rs, rv, cond, cells_);
V BW = 1./tempBW.value(); V BW = 1./tempBW.value();
V RW = residual_.material_balance_eq[Water].value(); V RW = residual_.material_balance_eq[Water].value();
BW_avg = BW.sum()/nc; BW_avg = BW.sum()/nc;
const V tempV = RW.abs()/pv; const V tempV = RW.abs()/pv;
CNVW = BW_avg * dt * tempV.maxCoeff(); CNVW = BW_avg * dt * tempV.maxCoeff();
RW_sum = RW.sum(); RW_sum = RW.sum();
std::cout << " CNVW " << CNVW << " RW_sum "<< RW_sum << std::endl; std::cout << " CNVW " << CNVW << " RW_sum "<< RW_sum << std::endl;
} }
if (active_[Oil]) { if (active_[Oil]) {
// Omit the disgas here. We should add it. // Omit the disgas here. We should add it.
ADB tempBO = fluidReciprocFVF(Oil, press, rs, rv, cond, cells_); ADB tempBO = fluidReciprocFVF(Oil, press, rs, rv, cond, cells_);
V BO = 1./tempBO.value(); V BO = 1./tempBO.value();
V RO = residual_.material_balance_eq[Oil].value(); V RO = residual_.material_balance_eq[Oil].value();
BO_avg = BO.sum()/nc; BO_avg = BO.sum()/nc;
const V tempV = RO.abs()/pv; const V tempV = RO.abs()/pv;
CNVO = BO_avg * dt * tempV.maxCoeff(); CNVO = BO_avg * dt * tempV.maxCoeff();
RO_sum = RO.sum(); RO_sum = RO.sum();
std::cout << " CNVO " << CNVO << " RO_sum " << RO_sum << std::endl; std::cout << " CNVO " << CNVO << " RO_sum " << RO_sum << std::endl;
} }
@ -1770,15 +1770,15 @@ namespace {
ADB tempBG = fluidReciprocFVF(Gas, press, rs, rv, cond, cells_); ADB tempBG = fluidReciprocFVF(Gas, press, rs, rv, cond, cells_);
V BG = 1./tempBG.value(); V BG = 1./tempBG.value();
V RG = residual_.material_balance_eq[Gas].value(); V RG = residual_.material_balance_eq[Gas].value();
BG_avg = BG.sum()/nc; BG_avg = BG.sum()/nc;
const V tempV = RG.abs()/pv; const V tempV = RG.abs()/pv;
CNVG = BG_avg * dt * tempV.maxCoeff(); CNVG = BG_avg * dt * tempV.maxCoeff();
RG_sum = RG.sum(); RG_sum = RG.sum();
std::cout << " CNVG " << CNVG << " RG_sum " << RG_sum << std::endl; std::cout << " CNVG " << CNVG << " RG_sum " << RG_sum << std::endl;
} }
double tempValue = tol_mb * pvSum /dt; double tempValue = tol_mb * pvSum /dt;
bool converged_MB = (fabs(BW_avg*RW_sum) < tempValue) bool converged_MB = (fabs(BW_avg*RW_sum) < tempValue)
&& (fabs(BO_avg*RO_sum) < tempValue) && (fabs(BO_avg*RO_sum) < tempValue)
@ -1788,7 +1788,7 @@ namespace {
double residualWellFlux = residual_.well_flux_eq.value().matrix().lpNorm<Eigen::Infinity>(); double residualWellFlux = residual_.well_flux_eq.value().matrix().lpNorm<Eigen::Infinity>();
double residualWell = residual_.well_eq.value().matrix().lpNorm<Eigen::Infinity>(); double residualWell = residual_.well_eq.value().matrix().lpNorm<Eigen::Infinity>();
const double day = 24 * 60 * 60; const double day = 24 * 60 * 60;
const double barsa = 1.e5; const double barsa = 1.e5;