c8ff9503a2
Foam cases should run without FOAMOPTS available. That means we need to default the mobility model to a valid value. This also adds enum value 'FUNC' and handles setting the default model if only transport phase is unspecified.
237 lines
7.3 KiB
C++
237 lines
7.3 KiB
C++
/*
|
|
Copyright 2019 SINTEF Digital, Mathematics and Cybernetics.
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#include <opm/parser/eclipse/Deck/Deck.hpp>
|
|
#include <opm/parser/eclipse/Deck/DeckItem.hpp>
|
|
#include <opm/parser/eclipse/Deck/DeckKeyword.hpp>
|
|
#include <opm/parser/eclipse/Deck/DeckRecord.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/InitConfig/FoamConfig.hpp>
|
|
#include <opm/parser/eclipse/Parser/ParserKeywords/F.hpp>
|
|
|
|
namespace Opm
|
|
{
|
|
|
|
// FoamData member functions.
|
|
|
|
FoamData::FoamData()
|
|
: FoamData(0.0,
|
|
ParserKeywords::FOAMFSC::EXPONENT::defaultValue,
|
|
ParserKeywords::FOAMFSC::MIN_SURF_CONC::defaultValue,
|
|
false, 1.0)
|
|
{
|
|
}
|
|
|
|
FoamData::FoamData(const DeckRecord& FOAMFSC_record, const DeckRecord& FOAMROCK_record)
|
|
: reference_surfactant_concentration_(FOAMFSC_record.getItem(0).getSIDouble(0))
|
|
, exponent_(FOAMFSC_record.getItem(1).getSIDouble(0))
|
|
, minimum_surfactant_concentration_(FOAMFSC_record.getItem(2).getSIDouble(0))
|
|
, allow_desorption_(true) // will be overwritten below
|
|
, rock_density_(FOAMROCK_record.getItem(1).getSIDouble(0))
|
|
{
|
|
// Check validity of adsorption index and set allow_desorption_ member.
|
|
const int ads_ind = FOAMROCK_record.getItem(0).get<int>(0);
|
|
if (ads_ind < 1 || ads_ind > 2) {
|
|
throw std::runtime_error("Illegal adsorption index in FOAMROCK, must be 1 or 2.");
|
|
}
|
|
allow_desorption_ = (ads_ind == 1);
|
|
}
|
|
|
|
FoamData::FoamData(const DeckRecord& FOAMROCK_record)
|
|
: reference_surfactant_concentration_(0.0)
|
|
, exponent_(0.0)
|
|
, minimum_surfactant_concentration_(1e-20)
|
|
, allow_desorption_(true) // will be overwritten below
|
|
, rock_density_(FOAMROCK_record.getItem(1).getSIDouble(0))
|
|
{
|
|
// Check validity of adsorption index and set allow_desorption_ member.
|
|
const int ads_ind = FOAMROCK_record.getItem(0).get<int>(0);
|
|
if (ads_ind < 1 || ads_ind > 2) {
|
|
throw std::runtime_error("Illegal adsorption index in FOAMROCK, must be 1 or 2.");
|
|
}
|
|
allow_desorption_ = (ads_ind == 1);
|
|
}
|
|
|
|
FoamData::FoamData(double reference_surfactant_concentration,
|
|
double exponent,
|
|
double minimum_surfactant_concentration,
|
|
bool allow_desorption,
|
|
double rock_density)
|
|
: reference_surfactant_concentration_(reference_surfactant_concentration)
|
|
, exponent_(exponent)
|
|
, minimum_surfactant_concentration_(minimum_surfactant_concentration)
|
|
, allow_desorption_(allow_desorption)
|
|
, rock_density_(rock_density)
|
|
{
|
|
}
|
|
|
|
double
|
|
FoamData::referenceSurfactantConcentration() const
|
|
{
|
|
return this->reference_surfactant_concentration_;
|
|
}
|
|
|
|
double
|
|
FoamData::exponent() const
|
|
{
|
|
return this->exponent_;
|
|
}
|
|
|
|
double
|
|
FoamData::minimumSurfactantConcentration() const
|
|
{
|
|
return this->minimum_surfactant_concentration_;
|
|
}
|
|
|
|
bool
|
|
FoamData::allowDesorption() const
|
|
{
|
|
return this->allow_desorption_;
|
|
}
|
|
|
|
double
|
|
FoamData::rockDensity() const
|
|
{
|
|
return this->rock_density_;
|
|
}
|
|
|
|
bool
|
|
FoamData::operator==(const FoamData& data) const
|
|
{
|
|
return reference_surfactant_concentration_ ==
|
|
data.reference_surfactant_concentration_ &&
|
|
exponent_ == data.exponent_ &&
|
|
minimum_surfactant_concentration_ ==
|
|
data.minimum_surfactant_concentration_ &&
|
|
allow_desorption_ == data.allow_desorption_ &&
|
|
rock_density_ == data.rock_density_;
|
|
}
|
|
|
|
// FoamConfig member functions.
|
|
|
|
FoamConfig::FoamConfig(const Deck& deck)
|
|
{
|
|
if (deck.hasKeyword<ParserKeywords::FOAMOPTS>()) {
|
|
// We only support the default (GAS transport phase, TAB mobility reduction model)
|
|
// setup for foam at this point, so we detect and deal with it here even though we
|
|
// do not store any data related to it.
|
|
const auto& kw_foamopts = deck.getKeyword<ParserKeywords::FOAMOPTS>();
|
|
transport_phase_ = get_phase(kw_foamopts.getRecord(0).getItem(0).get<std::string>(0));
|
|
std::string mobModel = kw_foamopts.getRecord(0).getItem(1).get<std::string>(0);
|
|
if (mobModel.empty()) {
|
|
if (transport_phase_ == Phase::GAS) {
|
|
mobility_model_ = MobilityModel::TAB;
|
|
} else if (transport_phase_ == Phase::WATER) {
|
|
mobility_model_ = MobilityModel::FUNC;
|
|
}
|
|
} else if (mobModel == "TAB") {
|
|
mobility_model_ = MobilityModel::TAB;
|
|
} else if (mobModel == "FUNC") {
|
|
mobility_model_ = MobilityModel::FUNC;
|
|
}
|
|
}
|
|
if (deck.hasKeyword<ParserKeywords::FOAMFSC>()) {
|
|
if (!deck.hasKeyword<ParserKeywords::FOAMROCK>()) {
|
|
throw std::runtime_error("FOAMFSC present but no FOAMROCK keyword found.");
|
|
}
|
|
// We have both FOAMFSC and FOAMROCK.
|
|
const auto& kw_foamfsc = deck.getKeyword<ParserKeywords::FOAMFSC>();
|
|
const auto& kw_foamrock = deck.getKeyword<ParserKeywords::FOAMROCK>();
|
|
if (kw_foamfsc.size() != kw_foamrock.size()) {
|
|
throw std::runtime_error("FOAMFSC and FOAMROCK keywords have different number of records.");
|
|
}
|
|
const int num_records = kw_foamfsc.size();
|
|
for (int record_index = 0; record_index < num_records; ++record_index) {
|
|
this->data_.emplace_back(kw_foamfsc.getRecord(record_index), kw_foamrock.getRecord(record_index));
|
|
}
|
|
} else if (deck.hasKeyword<ParserKeywords::FOAMROCK>()) {
|
|
// We have FOAMROCK, but not FOAMFSC.
|
|
const auto& kw_foamrock = deck.getKeyword<ParserKeywords::FOAMROCK>();
|
|
for (const auto& record : kw_foamrock) {
|
|
this->data_.emplace_back(record);
|
|
}
|
|
}
|
|
}
|
|
|
|
FoamConfig::FoamConfig(const std::vector<FoamData>& data,
|
|
Phase transport_phase,
|
|
MobilityModel mobility_model)
|
|
: data_(data)
|
|
, transport_phase_(transport_phase)
|
|
, mobility_model_(mobility_model)
|
|
{
|
|
}
|
|
|
|
const FoamData&
|
|
FoamConfig::getRecord(std::size_t index) const
|
|
{
|
|
return this->data_.at(index);
|
|
}
|
|
|
|
const std::vector<FoamData>&
|
|
FoamConfig::records() const
|
|
{
|
|
return this->data_;
|
|
}
|
|
|
|
std::size_t
|
|
FoamConfig::size() const
|
|
{
|
|
return this->data_.size();
|
|
}
|
|
|
|
bool
|
|
FoamConfig::empty() const
|
|
{
|
|
return this->data_.empty();
|
|
}
|
|
|
|
FoamConfig::const_iterator
|
|
FoamConfig::begin() const
|
|
{
|
|
return this->data_.begin();
|
|
}
|
|
|
|
FoamConfig::const_iterator
|
|
FoamConfig::end() const
|
|
{
|
|
return this->data_.end();
|
|
}
|
|
|
|
Phase
|
|
FoamConfig::getTransportPhase() const
|
|
{
|
|
return this->transport_phase_;
|
|
}
|
|
|
|
FoamConfig::MobilityModel
|
|
FoamConfig::getMobilityModel() const
|
|
{
|
|
return this->mobility_model_;
|
|
}
|
|
|
|
bool
|
|
FoamConfig::operator==(const FoamConfig& data) const
|
|
{
|
|
return transport_phase_ == data.transport_phase_ &&
|
|
mobility_model_ == data.mobility_model_ &&
|
|
data_ == data.data_;
|
|
}
|
|
|
|
} // namespace Opm
|