OilfieldToolsNet/opm-upscaling
4
0
Fork 0
mirror of https://github.com/OPM/opm-upscaling.git synced 2025-02-25 18:45:23 -06:00
Code Issues Packages Projects Releases Wiki Activity

Initialize system with s_min < s < s_max.

Browse Source
This commit is contained in:
Ingeborg Ligaarden
2011-11-18 17:51:25 +01:00
parent 865a342306
commit cf96539d7b
1 changed files with 9 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
dune/upscaling/SteadyStateUpscalerImplicit_impl.hpp
View File

@@ -154,6 +154,13 @@ namespace Dune
// Set up initial saturation profile.
std::vector<double> saturation = initial_saturation;
for (int c = 0; c < saturation.size(); c++ ) {
double s_min = this->res_prop_.s_min(c);
double s_max = this->res_prop_.s_max(c);
saturation[c] = std::max(s_min, saturation[c]);
saturation[c] = std::min(s_max, saturation[c]);
}
// Set up boundary conditions.
setupUpscalingConditions(this->ginterf_, this->bctype_, flow_direction,
pressure_drop, boundary_saturation, this->twodim_hack_, this->bcond_);
@@ -175,6 +182,8 @@ namespace Dune
bool stationary = false;
int it_count=0;
double stepsize=init_stepsize_;
std::vector<double> init_saturation(saturation);
while((!stationary) & (it_count < max_it_)){// & transport_cost < max_transport_cost_)
//for (int iter = 0; iter < simulation_steps_; ++iter) {