OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-27 11:06:25 -06:00
Code Issues Packages Projects Releases Wiki Activity

Initialized well connection pressures first

Browse Source 
in assemble(...). This makes VFPINJ behave as expected, and
  VFPPROD for the "trivial table". For the nontrivial table,
  VFPPROD does not match expected behaviour.
This commit is contained in:
babrodtk 2015-08-17 13:05:32 +02:00
parent 657a7c58b8
commit 3fddf86eef
1 changed files with 14 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -774,6 +774,19 @@ namespace detail {
{
using namespace Opm::AutoDiffGrid;
if (true) {
// Create the primary variables.
SolutionState state = asImpl().variableState(reservoir_state, well_state);
if (true || initial_assembly) {
// Create the (constant, derivativeless) initial state.
SolutionState state0 = state;
asImpl().makeConstantState(state0);
computeWellConnectionPressures(state0, well_state);
}
}
// Possibly switch well controls and updating well state to
// get reasonable initial conditions for the wells
updateWellControls(well_state);
@ -2307,7 +2320,7 @@ namespace detail {
std::isnan(well_flux_residual[Water]) || well_flux_residual[Water] > maxResidualAllowed() ||
std::isnan(well_flux_residual[Oil]) || well_flux_residual[Oil] > maxResidualAllowed() ||
std::isnan(well_flux_residual[Gas]) || well_flux_residual[Gas] > maxResidualAllowed() ||
std::isnan(residualWell) || residualWell > maxResidualAllowed() )
std::isnan(residualWell) || residualWell > 100*maxResidualAllowed() )
{
OPM_THROW(Opm::NumericalProblem,"One of the residuals is NaN or too large!");
}