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

renamed saltwater-->brine and removed keywords from missingFeatures

Browse Source
This commit is contained in:
Trine S. Mykkeltvedt 2019-12-03 17:48:17 +01:00 committed by Tor Harald Sandve
parent 5749615bd8
commit 99e7785fb5
20 changed files with 80 additions and 83 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
CMakeLists.txt
View File

@ -168,7 +168,7 @@ opm_add_test(flow
flow/flow_ebos_oilwater.cpp
flow/flow_ebos_polymer.cpp
flow/flow_ebos_foam.cpp
flow/flow_ebos_saltwater.cpp
flow/flow_ebos_brine.cpp
flow/flow_ebos_solvent.cpp
flow/flow_ebos_energy.cpp
flow/flow_ebos_oilwater_polymer.cpp
@ -224,7 +224,7 @@ endif()
# the research oriented general-purpose ECL simulator ("ebos" == &ecl
# &black-&oil &simulator)
set(MEBOS_TARGETS "")
foreach(OBJ blackoil solvent polymer foam gasoil oilwater oilwaterpolymer thermal)
foreach(OBJ blackoil solvent polymer foam brine gasoil oilwater oilwaterpolymer thermal)
add_library(ebos_lib${OBJ} OBJECT EXCLUDE_FROM_ALL ebos/ebos_${OBJ}.cc)
list(APPEND MEBOS_TARGETS $<TARGET_OBJECTS:ebos_lib${OBJ}>)
endforeach()
@ -248,7 +248,7 @@ else()
set(EBOS_EXTENSIONS_DEFAULT_ENABLE_IF "TRUE")
endif()
foreach(OBJ solvent polymer foam saltwater gasoil oilwater oilwaterpolymer thermal)
foreach(OBJ solvent polymer foam brine gasoil oilwater oilwaterpolymer thermal)
opm_add_test(ebos_${OBJ}
ONLY_COMPILE
DEFAULT_ENABLE_IF ${EBOS_EXTENSIONS_DEFAULT_ENABLE_IF}
@ -288,7 +288,7 @@ opm_add_test(ebos_plain
LIBRARIES opmsimulators)
if (BUILD_EBOS_EXTENSIONS)
foreach(TGT ebos_solvent ebos_polymer ebos_foam ebos_saltwater ebos_gasoil ebos_oilwater ebos_oilwaterpolymer ebos_thermal mebos)
foreach(TGT ebos_solvent ebos_polymer ebos_foam ebos_brine ebos_saltwater ebos_gasoil ebos_oilwater ebos_oilwaterpolymer ebos_thermal mebos)
install(TARGETS ${TGT} DESTINATION bin)
opm_add_bash_completion(${TGT})
endforeach()

14
ebos/ebos_saltwater.cc → ebos/ebos_brine.cc
View File

@ -32,21 +32,21 @@
BEGIN_PROPERTIES
NEW_TYPE_TAG(EbosSaltWaterTypeTag, INHERITS_FROM(EbosTypeTag));
NEW_TYPE_TAG(EbosBrineTypeTag, INHERITS_FROM(EbosTypeTag));
// enable the salt extension of the black oil model
SET_BOOL_PROP(EbosSaltWaterTypeTag, EnableSaltWater, true);
// enable the brine extension of the black oil model
SET_BOOL_PROP(EbosBrineTypeTag, EnableBrine, true);
END_PROPERTIES
namespace Opm {
void ebosSaltWaterSetDeck(Opm::Deck* deck,
void ebosBrineSetDeck(Opm::Deck* deck,
Opm::ParseContext* parseContext,
Opm::ErrorGuard* errorGuard,
double externalSetupTime)
{
typedef TTAG(EbosSaltWaterTypeTag) ProblemTypeTag;
typedef TTAG(EbosBrineTypeTag) ProblemTypeTag;
typedef GET_PROP_TYPE(ProblemTypeTag, Vanguard) Vanguard;
Vanguard::setExternalSetupTime(externalSetupTime);
@ -55,9 +55,9 @@ void ebosSaltWaterSetDeck(Opm::Deck* deck,
Vanguard::setExternalDeck(deck);
}
int ebosSaltWaterMain(int argc, char **argv)
int ebosBrineMain(int argc, char **argv)
{
typedef TTAG(EbosSaltWaterTypeTag) ProblemTypeTag;
typedef TTAG(EbosBrineTypeTag) ProblemTypeTag;
return Opm::startEbos<ProblemTypeTag>(argc, argv);
}

10
ebos/ebos_saltwater.hh → ebos/ebos_brine.hh
View File

@ -23,22 +23,22 @@
/*!
* \file
*
* \brief The function prototypes required to start the salt variant of ebos
* \brief The function prototypes required to start the brine variant of ebos
*/
#ifndef EBOS_SALTWATER_HH
#define EBOS_SALTWATER_HH
#ifndef EBOS_BRINE_HH
#define EBOS_BRINE_HH
#include <opm/parser/eclipse/Deck/Deck.hpp>
#include <opm/parser/eclipse/Parser/ParseContext.hpp>
#include <opm/parser/eclipse/Parser/ErrorGuard.hpp>
namespace Opm {
void ebosSaltWaterSetDeck(Opm::Deck* deck,
void ebosBrineSetDeck(Opm::Deck* deck,
Opm::ParseContext* parseContext,
Opm::ErrorGuard* errorGuard,
double externalSetupTime);
int ebosSaltWaterMain(int argc, char** argv);
int ebosBrineMain(int argc, char** argv);
}
#endif

8
ebos/ebos_saltwater_main.cc → ebos/ebos_brine_main.cc
View File

@ -23,15 +23,15 @@
/*!
* \file
*
* \brief The main function for the stand alone salt variant of ebos.
* \brief The main function for the stand alone brine variant of ebos.
*
* This only calls the ebosSaltWaterMain() function.
* This only calls the ebosBrineMain() function.
*/
#include "config.h"
#include "ebos_saltwater.hh"
#include "ebos_brine.hh"
int main(int argc, char** argv)
{
return Opm::ebosSaltWaterMain(argc, argv);
return Opm::ebosBrineMain(argc, argv);
}

4
ebos/ecloutputblackoilmodule.hh
View File

@ -312,7 +312,7 @@ public:
cPolymer_.resize(bufferSize, 0.0);
if (GET_PROP_VALUE(TypeTag, EnableFoam))
cFoam_.resize(bufferSize, 0.0);
if (GET_PROP_VALUE(TypeTag, EnableSaltWater))
if (GET_PROP_VALUE(TypeTag, EnableBrine))
cSalt_.resize(bufferSize, 0.0);
if (simulator_.problem().vapparsActive())
@ -570,7 +570,7 @@ public:
}
if (cSalt_.size() > 0) {
cSalt_[globalDofIdx] = fs.saltconcentration().value();
cSalt_[globalDofIdx] = fs.saltConcentration().value();
}
if (bubblePointPressure_.size() > 0) {

12
flow/flow.cpp
View File

@ -26,7 +26,7 @@
#include <flow/flow_ebos_solvent.hpp>
#include <flow/flow_ebos_polymer.hpp>
#include <flow/flow_ebos_foam.hpp>
#include <flow/flow_ebos_saltwater.hpp>
#include <flow/flow_ebos_brine.hpp>
#include <flow/flow_ebos_energy.hpp>
#include <flow/flow_ebos_oilwater_polymer.hpp>
#include <flow/flow_ebos_oilwater_polymer_injectivity.hpp>
@ -429,10 +429,10 @@ int main(int argc, char** argv)
Opm::flowEbosFoamSetDeck(externalSetupTimer.elapsed(), *deck, *eclipseState, *schedule, *summaryConfig);
return Opm::flowEbosFoamMain(argc, argv, outputCout, outputFiles);
}
// Salt case
else if ( phases.active( Opm::Phase::SALTWATER ) ) {
Opm::flowEbosSaltWaterSetDeck(externalSetupTimer.elapsed(), *deck, *eclipseState, *schedule, *summaryConfig);
return Opm::flowEbosSaltWaterMain(argc, argv, outputCout, outputFiles);
// Brine case
else if ( phases.active( Opm::Phase::BRINE ) ) {
Opm::flowEbosBrineSetDeck(externalSetupTimer.elapsed(), *deck, *eclipseState, *schedule, *summaryConfig);
return Opm::flowEbosBrineMain(argc, argv, outputCout, outputFiles);
}
// Solvent case
else if ( phases.active( Opm::Phase::SOLVENT ) ) {
@ -452,7 +452,7 @@ int main(int argc, char** argv)
else
{
if (outputCout)
std::cerr << "No suitable configuration found, valid are Twophase, polymer, foam, salt, solvent, energy, blackoil." << std::endl;
std::cerr << "No suitable configuration found, valid are Twophase, polymer, foam, brine, solvent, energy, blackoil." << std::endl;
return EXIT_FAILURE;
}
}

14
flow/flow_ebos_saltwater.cpp → flow/flow_ebos_brine.cpp
View File

@ -16,7 +16,7 @@
*/
#include "config.h"
#include <flow/flow_ebos_saltwater.hpp>
#include <flow/flow_ebos_brine.hpp>
#include <opm/material/common/ResetLocale.hpp>
#include <opm/grid/CpGrid.hpp>
@ -31,14 +31,14 @@
namespace Opm {
namespace Properties {
NEW_TYPE_TAG(EclFlowSaltWaterProblem, INHERITS_FROM(EclFlowProblem));
SET_BOOL_PROP(EclFlowSaltWaterProblem, EnableSaltWater, true);
NEW_TYPE_TAG(EclFlowBrineProblem, INHERITS_FROM(EclFlowProblem));
SET_BOOL_PROP(EclFlowBrineProblem, EnableBrine, true);
}}
namespace Opm {
void flowEbosSaltWaterSetDeck(double setupTime, Deck &deck, EclipseState& eclState, Schedule& schedule, SummaryConfig& summaryConfig)
void flowEbosBrineSetDeck(double setupTime, Deck &deck, EclipseState& eclState, Schedule& schedule, SummaryConfig& summaryConfig)
{
typedef TTAG(EclFlowSaltWaterProblem) TypeTag;
typedef TTAG(EclFlowBrineProblem) TypeTag;
typedef GET_PROP_TYPE(TypeTag, Vanguard) Vanguard;
Vanguard::setExternalSetupTime(setupTime);
@ -50,7 +50,7 @@ void flowEbosSaltWaterSetDeck(double setupTime, Deck &deck, EclipseState& eclSta
// ----------------- Main program -----------------
int flowEbosSaltWaterMain(int argc, char** argv, bool outputCout, bool outputFiles)
int flowEbosBrineMain(int argc, char** argv, bool outputCout, bool outputFiles)
{
// we always want to use the default locale, and thus spare us the trouble
// with incorrect locale settings.
@ -63,7 +63,7 @@ int flowEbosSaltWaterMain(int argc, char** argv, bool outputCout, bool outputFil
Dune::MPIHelper::instance(argc, argv).rank();
#endif
Opm::FlowMainEbos<TTAG(EclFlowSaltWaterProblem)> mainfunc;
Opm::FlowMainEbos<TTAG(EclFlowBrineProblem)> mainfunc;
return mainfunc.execute(argc, argv, outputCout, outputFiles);
}

10
flow/flow_ebos_saltwater.hpp → flow/flow_ebos_brine.hpp
View File

@ -14,8 +14,8 @@
You should have received a copy of the GNU General Public License
along with OPM. If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef FLOW_EBOS_SALTWATER_HPP
#define FLOW_EBOS_SALTWATER_HPP
#ifndef FLOW_EBOS_BRINE_HPP
#define FLOW_EBOS_BRINE_HPP
#include <opm/parser/eclipse/Deck/Deck.hpp>
#include <opm/parser/eclipse/EclipseState/EclipseState.hpp>
@ -24,8 +24,8 @@
namespace Opm {
void flowEbosSaltWaterSetDeck(double setupTime, Deck &deck, EclipseState& eclState, Schedule& schedule, SummaryConfig& summaryConfig);
int flowEbosSaltWaterMain(int argc, char** argv, bool outputCout, bool outputFiles);
void flowEbosBrineSetDeck(double setupTime, Deck &deck, EclipseState& eclState, Schedule& schedule, SummaryConfig& summaryConfig);
int flowEbosBrineMain(int argc, char** argv, bool outputCout, bool outputFiles);
}
#endif // FLOW_EBOS_SALTWATER_HPP
#endif // FLOW_EBOS_BRINE_HPP

2
flow/flow_ebos_gasoil.cpp
View File

@ -53,7 +53,7 @@ public:
GET_PROP_VALUE(TypeTag, EnablePolymer),
GET_PROP_VALUE(TypeTag, EnableEnergy),
GET_PROP_VALUE(TypeTag, EnableFoam),
GET_PROP_VALUE(TypeTag, EnableSaltWater),
GET_PROP_VALUE(TypeTag, EnableBrine),
/*PVOffset=*/0,
/*disabledCompIdx=*/FluidSystem::waterCompIdx> type;
};

2
flow/flow_ebos_oilwater.cpp
View File

@ -53,7 +53,7 @@ public:
GET_PROP_VALUE(TypeTag, EnablePolymer),
GET_PROP_VALUE(TypeTag, EnableEnergy),
GET_PROP_VALUE(TypeTag, EnableFoam),
GET_PROP_VALUE(TypeTag, EnableSaltWater),
GET_PROP_VALUE(TypeTag, EnableBrine),
/*PVOffset=*/0,
/*disabledCompIdx=*/FluidSystem::gasCompIdx> type;
};

2
flow/flow_ebos_oilwater_polymer.cpp
View File

@ -54,7 +54,7 @@ public:
GET_PROP_VALUE(TypeTag, EnablePolymer),
GET_PROP_VALUE(TypeTag, EnableEnergy),
GET_PROP_VALUE(TypeTag, EnableFoam),
GET_PROP_VALUE(TypeTag, EnableSaltWater),
GET_PROP_VALUE(TypeTag, EnableBrine),
/*PVOffset=*/0,
/*disabledCompIdx=*/FluidSystem::gasCompIdx> type;
};

2
flow/flow_ebos_oilwater_polymer_injectivity.cpp
View File

@ -56,7 +56,7 @@ public:
2,
0,
GET_PROP_VALUE(TypeTag, EnableFoam),
GET_PROP_VALUE(TypeTag, EnableSaltWater),
GET_PROP_VALUE(TypeTag, EnableBrine),
/*PVOffset=*/0,
/*disabledCompIdx=*/FluidSystem::gasCompIdx> type;
};

6
opm/core/props/BlackoilPhases.hpp
View File

@ -33,11 +33,11 @@ namespace Opm
// sense that they can be active or not and canonical indices can be translated
// to and from active ones. That said, they are not considered by num_phases or
// MaxNumPhases. The crypto phases which are currently implemented are solvent,
// polymer, energy, polymer molecular weight, foam and salt.
// polymer, energy, polymer molecular weight, foam and brine.
static const int NumCryptoPhases = 6;
// enum ComponentIndex { Water = 0, Oil = 1, Gas = 2 };
enum PhaseIndex { Aqua = 0, Liquid = 1, Vapour = 2, Solvent = 3, Polymer = 4, Energy = 5, PolymerMW = 6, Foam = 7, Salt = 8 };
enum PhaseIndex { Aqua = 0, Liquid = 1, Vapour = 2, Solvent = 3, Polymer = 4, Energy = 5, PolymerMW = 6, Foam = 7, Brine = 8 };
};
struct PhaseUsage : public BlackoilPhases
@ -51,7 +51,7 @@ namespace Opm
// polymer molecular weight
bool has_polymermw;
bool has_foam;
bool has_salt;
bool has_brine;
};
/// Check or assign presence of a formed, free phase. Limited to

10
opm/core/props/phaseUsageFromDeck.hpp
View File

@ -123,14 +123,14 @@ namespace Opm
else
pu.phase_pos[BlackoilPhases::Foam] = -1;
// Add salt info
pu.has_salt = phases.active(Phase::SALTWATER);
if (pu.has_salt) {
pu.phase_pos[BlackoilPhases::Salt] = numActivePhases;
// Add brine info
pu.has_brine = phases.active(Phase::BRINE);
if (pu.has_brine) {
pu.phase_pos[BlackoilPhases::Brine] = numActivePhases;
++ numActivePhases;
}
else
pu.phase_pos[BlackoilPhases::Salt] = -1;
pu.phase_pos[BlackoilPhases::Brine] = -1;
return pu;
}

24
opm/simulators/flow/BlackoilModelEbos.hpp
View File

@ -84,7 +84,7 @@ SET_BOOL_PROP(EclFlowProblem, EnableSolvent, false);
SET_BOOL_PROP(EclFlowProblem, EnableTemperature, true);
SET_BOOL_PROP(EclFlowProblem, EnableEnergy, false);
SET_BOOL_PROP(EclFlowProblem, EnableFoam, false);
SET_BOOL_PROP(EclFlowProblem, EnableSaltWater, false);
SET_BOOL_PROP(EclFlowProblem, EnableBrine, false);
SET_TYPE_PROP(EclFlowProblem, EclWellModel, Opm::BlackoilWellModel<TypeTag>);
SET_TAG_PROP(EclFlowProblem, LinearSolverSplice, FlowIstlSolver);
@ -125,13 +125,13 @@ namespace Opm {
static const int contiEnergyEqIdx = Indices::contiEnergyEqIdx;
static const int contiPolymerMWEqIdx = Indices::contiPolymerMWEqIdx;
static const int contiFoamEqIdx = Indices::contiFoamEqIdx;
static const int contiSaltWaterEqIdx = Indices::contiSaltWaterEqIdx;
static const int contiBrineEqIdx = Indices::contiBrineEqIdx;
static const int solventSaturationIdx = Indices::solventSaturationIdx;
static const int polymerConcentrationIdx = Indices::polymerConcentrationIdx;
static const int polymerMoleWeightIdx = Indices::polymerMoleWeightIdx;
static const int temperatureIdx = Indices::temperatureIdx;
static const int foamConcentrationIdx = Indices::foamConcentrationIdx;
static const int saltWaterConcentrationIdx = Indices::saltConcentrationIdx;
static const int saltConcentrationIdx = Indices::saltConcentrationIdx;
typedef Dune::FieldVector<Scalar, numEq > VectorBlockType;
typedef typename SparseMatrixAdapter::MatrixBlock MatrixBlockType;
@ -167,7 +167,7 @@ namespace Opm {
, has_polymermw_(GET_PROP_VALUE(TypeTag, EnablePolymerMW))
, has_energy_(GET_PROP_VALUE(TypeTag, EnableEnergy))
, has_foam_(GET_PROP_VALUE(TypeTag, EnableFoam))
, has_salt_(GET_PROP_VALUE(TypeTag, EnableSaltWater))
, has_brine_(GET_PROP_VALUE(TypeTag, EnableBrine))
, param_( param )
, well_model_ (well_model)
, terminal_output_ (terminal_output)
@ -635,11 +635,11 @@ namespace Opm {
R_sum[ contiFoamEqIdx ] += R2;
maxCoeff[ contiFoamEqIdx ] = std::max( maxCoeff[ contiFoamEqIdx ], std::abs( R2 ) / pvValue );
}
if (has_salt_ ) {
B_avg[ contiSaltWaterEqIdx ] += 1.0 / fs.invB(FluidSystem::gasPhaseIdx).value();
const auto R2 = ebosResid[cell_idx][contiSaltWaterEqIdx];
R_sum[ contiSaltWaterEqIdx ] += R2;
maxCoeff[ contiSaltWaterEqIdx ] = std::max( maxCoeff[ contiSaltWaterEqIdx ], std::abs( R2 ) / pvValue );
if (has_brine_ ) {
B_avg[ contiBrineEqIdx ] += 1.0 / fs.invB(FluidSystem::gasPhaseIdx).value();
const auto R2 = ebosResid[cell_idx][contiBrineEqIdx];
R_sum[ contiBrineEqIdx ] += R2;
maxCoeff[ contiBrineEqIdx ] = std::max( maxCoeff[ contiBrineEqIdx ], std::abs( R2 ) / pvValue );
}
if (has_polymermw_) {
@ -730,8 +730,8 @@ namespace Opm {
if (has_foam_) {
compNames[foamConcentrationIdx] = "Foam";
}
if (has_salt_) {
compNames[saltWaterConcentrationIdx] = "Salt";
if (has_brine_) {
compNames[saltConcentrationIdx] = "Brine";
}
}
@ -885,7 +885,7 @@ namespace Opm {
const bool has_polymermw_;
const bool has_energy_;
const bool has_foam_;
const bool has_salt_;
const bool has_brine_;
ModelParameters param_;
SimulatorReport failureReport_;

3
opm/simulators/flow/MissingFeatures.cpp
View File

@ -124,7 +124,6 @@ namespace MissingFeatures {
"BPARA",
"BPIDIMS",
"BRANPROP",
"BRINE",
"BTOBALFA",
"BTOBALFV",
"CALTRAC",
@ -627,7 +626,6 @@ namespace MissingFeatures {
"SALT",
"SALTNODE",
"SALTREST",
"SALTVD",
"SCALELIM",
"SCDATAB",
"SCDETAB",
@ -809,7 +807,6 @@ namespace MissingFeatures {
"WPOLYRED",
"WPOTCALC",
"WREGROUP",
"WSALT",
"WSCCLEAN",
"WSCCLENL",
"WSCTAB",

12
opm/simulators/wells/StandardWell.hpp
View File

@ -31,7 +31,7 @@
#include <opm/models/blackoil/blackoilpolymermodules.hh>
#include <opm/models/blackoil/blackoilsolventmodules.hh>
#include <opm/models/blackoil/blackoilfoammodules.hh>
#include <opm/models/blackoil/blackoilsaltwatermodules.hh>
#include <opm/models/blackoil/blackoilbrinemodules.hh>
#include <opm/material/densead/DynamicEvaluation.hpp>
@ -69,22 +69,22 @@ namespace Opm
using Base::has_solvent;
using Base::has_polymer;
using Base::has_foam;
using Base::has_salt;
using Base::has_brine;
using Base::has_energy;
using PolymerModule = Opm::BlackOilPolymerModule<TypeTag>;
using FoamModule = Opm::BlackOilFoamModule<TypeTag>;
using SaltWaterModule = Opm::BlackOilSaltWaterModule<TypeTag>;
using BrineModule = Opm::BlackOilBrineModule<TypeTag>;
// polymer concentration and temperature are already known by the well, so
// polymer and energy conservation do not need to be considered explicitly
static const int numPolymerEq = Indices::numPolymers;
static const int numEnergyEq = Indices::numEnergy;
static const int numFoamEq = Indices::numFoam;
static const int numSaltEq = Indices::numSaltWater;
static const int numBrineEq = Indices::numBrine;
// number of the conservation equations
static const int numWellConservationEq = numEq - numPolymerEq - numEnergyEq - numFoamEq - numSaltEq;
static const int numWellConservationEq = numEq - numPolymerEq - numEnergyEq - numFoamEq - numBrineEq;
// number of the well control equations
static const int numWellControlEq = 1;
// number of the well equations that will always be used
@ -142,7 +142,7 @@ namespace Opm
using Base::contiSolventEqIdx;
using Base::contiPolymerEqIdx;
using Base::contiFoamEqIdx;
using Base::contiSaltWaterEqIdx;
using Base::contiBrineEqIdx;
static const int contiEnergyEqIdx = Indices::contiEnergyEqIdx;
StandardWell(const Well& well, const int time_step,

8
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -691,16 +691,16 @@ namespace Opm
connectionRates_[perf][contiFoamEqIdx] = Base::restrictEval(cq_s_foam);
}
if (has_salt) {
if (has_brine) {
// TODO: the application of well efficiency factor has not been tested with an example yet
const unsigned waterCompIdx = Indices::canonicalToActiveComponentIndex(FluidSystem::waterCompIdx);
EvalWell cq_s_sm = cq_s[waterCompIdx] * well_efficiency_factor_;
if (this->isInjector()) {
cq_s_sm *= wsalt();
} else {
cq_s_sm *= extendEval(intQuants.fluidState().saltconcentration());
cq_s_sm *= extendEval(intQuants.fluidState().saltConcentration());
}
connectionRates_[perf][contiSaltWaterEqIdx] = Base::restrictEval(cq_s_sm);
connectionRates_[perf][contiBrineEqIdx] = Base::restrictEval(cq_s_sm);
}
// Store the perforation pressure for later usage.
@ -2528,7 +2528,7 @@ namespace Opm
{
// the following implementation assume that the polymer is always after the w-o-g phases
// For the polymer, energy and foam cases, there is one more mass balance equations of reservoir than wells
assert((int(B_avg.size()) == num_components_) || has_polymer || has_energy || has_foam || has_salt);
assert((int(B_avg.size()) == num_components_) || has_polymer || has_energy || has_foam || has_brine);
const double tol_wells = param_.tolerance_wells_;
const double maxResidualAllowed = param_.max_residual_allowed_;

6
opm/simulators/wells/WellInterface.hpp
View File

@ -98,13 +98,13 @@ namespace Opm
// flag for polymer molecular weight related
static const bool has_polymermw = GET_PROP_VALUE(TypeTag, EnablePolymerMW);
static const bool has_foam = GET_PROP_VALUE(TypeTag, EnableFoam);
static const bool has_salt = GET_PROP_VALUE(TypeTag, EnableSaltWater);
static const bool has_brine = GET_PROP_VALUE(TypeTag, EnableBrine);
static const int contiSolventEqIdx = Indices::contiSolventEqIdx;
static const int contiPolymerEqIdx = Indices::contiPolymerEqIdx;
// index for the polymer molecular weight continuity equation
static const int contiPolymerMWEqIdx = Indices::contiPolymerMWEqIdx;
static const int contiFoamEqIdx = Indices::contiFoamEqIdx;
static const int contiSaltWaterEqIdx = Indices::contiSaltWaterEqIdx;
static const int contiBrineEqIdx = Indices::contiBrineEqIdx;
// For the conversion between the surface volume rate and reservoir voidage rate
using RateConverterType = RateConverter::
@ -115,7 +115,7 @@ namespace Opm
has_temperature,
has_energy,
compositionSwitchEnabled,
has_salt,
has_brine,
Indices::numPhases >;
/// Constructor
WellInterface(const Well& well, const int time_step,

6
opm/simulators/wells/WellInterface_impl.hpp
View File

@ -378,15 +378,15 @@ namespace Opm
WellInterface<TypeTag>::
wsalt() const
{
if (!has_salt) {
if (!has_brine) {
return 0.0;
}
auto injectorType = well_ecl_.injectorType();
if (injectorType == Well::InjectorType::WATER) {
WellSaltwaterProperties fprop = well_ecl_.getSaltwaterProperties();
return fprop.m_saltwaterConcentration;
WellBrineProperties fprop = well_ecl_.getBrineProperties();
return fprop.m_saltConcentration;
} else {
// Not a water injection well => no salt (?).
return 0.0;