OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

added: basic tests for BCCON / BCPROP

Browse Source
This commit is contained in:
Arne Morten Kvarving
2023-08-16 12:57:12 +02:00
parent 4f3d22f655
commit bd3328b15c
2 changed files with 192 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
CMakeLists_files.cmake
View File

@@ -482,6 +482,7 @@ if(ENABLE_ECL_INPUT)
tests/parser/ACTIONX.cpp
tests/parser/ADDREGTests.cpp
tests/parser/AquiferTests.cpp
tests/parser/BCConfigTests.cpp
tests/parser/BoxTests.cpp
tests/parser/CarfinTests.cpp
tests/parser/ColumnSchemaTests.cpp

191
tests/parser/BCConfigTests.cpp Normal file
View File

@@ -0,0 +1,191 @@
/*
Copyright (C) 2023 Equinor
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/>.
*/
#define BOOST_TEST_MODULE BCConfigTests
#include <boost/test/unit_test.hpp>
#include <opm/input/eclipse/Deck/Deck.hpp>
#include <opm/input/eclipse/EclipseState/SimulationConfig/BCConfig.hpp>
#include <opm/input/eclipse/Parser/Parser.hpp>
#include <opm/input/eclipse/Parser/ParserKeywords/B.hpp>
#include <opm/input/eclipse/Units/UnitSystem.hpp>
#include <opm/input/eclipse/Schedule/BCProp.hpp>
#include <string>
namespace {
Opm::Deck createDeck(const std::string& input)
{
Opm::Parser parser;
auto deck = parser.parseString(input);
return deck;
}
}
BOOST_AUTO_TEST_CASE(GasRateFree)
{
const std::string input = R"(
RUNSPEC
DIMENS
10 10 3 /
OIL
GAS
WATER
LAB
START
1 'JAN' 2015 /
GRID
DX
300*1000 /
DY
300*1000 /
DZ
300*1000 /
TOPS
100*8325 /
BCCON
1 1 1 2* 5 10 X- /
2 20 20 4* X /
/
SCHEDULE
BCPROP
1 RATE GAS -0.01 /
2 FREE /
/
)";
auto deck = createDeck(input);
Opm::BCConfig config(deck);
const auto& kw = deck.get<Opm::ParserKeywords::BCPROP>();
Opm::BCProp prop;
for (const auto& record : kw.back()) {
prop.updateBCProp(record);
}
BOOST_CHECK_EQUAL(config.size(), 2U);
const auto& c1 = *config.begin();
BOOST_CHECK_EQUAL(c1.index, 1);
BOOST_CHECK_EQUAL(c1.i1, 0);
BOOST_CHECK_EQUAL(c1.i2, 0);
BOOST_CHECK_EQUAL(c1.j1, 0);
BOOST_CHECK_EQUAL(c1.j2, 9);
BOOST_CHECK_EQUAL(c1.k1, 4);
BOOST_CHECK_EQUAL(c1.k2, 9);
BOOST_CHECK_EQUAL(c1.dir, Opm::FaceDir::XMinus);
const auto& c2 = *(config.begin()+1);
BOOST_CHECK_EQUAL(c2.index, 2);
BOOST_CHECK_EQUAL(c2.i1, 19);
BOOST_CHECK_EQUAL(c2.i2, 19);
BOOST_CHECK_EQUAL(c2.j1, 0);
BOOST_CHECK_EQUAL(c2.j2, 9);
BOOST_CHECK_EQUAL(c2.k1, 0);
BOOST_CHECK_EQUAL(c2.k2, 2);
BOOST_CHECK_EQUAL(c2.dir, Opm::FaceDir::XPlus);
BOOST_CHECK_EQUAL(prop.size(), 2U);
BOOST_CHECK(prop[1].bctype == Opm::BCType::RATE);
BOOST_CHECK(prop[1].component == Opm::BCComponent::GAS);
BOOST_CHECK_EQUAL(prop[1].rate,
deck.getActiveUnitSystem().to_si("Mass/Time*Length*Length", -0.01));
BOOST_CHECK(prop[2].bctype == Opm::BCType::FREE);
BOOST_CHECK(prop[2].component == Opm::BCComponent::NONE);
BOOST_CHECK_EQUAL(prop[2].rate, 0.0);
}
BOOST_AUTO_TEST_CASE(DirichletThermal)
{
const std::string input = R"(
RUNSPEC
DIMENS
10 10 3 /
OIL
GAS
WATER
LAB
START
1 'JAN' 2015 /
GRID
DX
300*1000 /
DY
300*1000 /
DZ
300*1000 /
TOPS
100*8325 /
BCCON
1 1* * 1 5 1 10 Y- /
2 20 20 4* Y /
/
SCHEDULE
BCPROP
1 THERMAL WATER -0.05 * 50.0 /
2 DIRICHLET OIL -0.01 279.0 /
/
)";
auto deck = createDeck(input);
Opm::BCConfig config(deck);
const auto& kw = deck.get<Opm::ParserKeywords::BCPROP>();
Opm::BCProp prop;
for (const auto& record : kw.back()) {
prop.updateBCProp(record);
}
BOOST_CHECK_EQUAL(config.size(), 2U);
const auto& c1 = *config.begin();
BOOST_CHECK_EQUAL(c1.index, 1);
BOOST_CHECK_EQUAL(c1.i1, 0);
BOOST_CHECK_EQUAL(c1.i2, 9);
BOOST_CHECK_EQUAL(c1.j1, 0);
BOOST_CHECK_EQUAL(c1.j2, 4);
BOOST_CHECK_EQUAL(c1.k1, 0);
BOOST_CHECK_EQUAL(c1.k2, 9);
BOOST_CHECK_EQUAL(c1.dir, Opm::FaceDir::YMinus);
const auto& c2 = *(config.begin()+1);
BOOST_CHECK_EQUAL(c2.index, 2);
BOOST_CHECK_EQUAL(c2.i1, 19);
BOOST_CHECK_EQUAL(c2.i2, 19);
BOOST_CHECK_EQUAL(c2.j1, 0);
BOOST_CHECK_EQUAL(c2.j2, 9);
BOOST_CHECK_EQUAL(c2.k1, 0);
BOOST_CHECK_EQUAL(c2.k2, 2);
BOOST_CHECK_EQUAL(c2.dir, Opm::FaceDir::YPlus);
BOOST_CHECK_EQUAL(prop.size(), 2U);
BOOST_CHECK(prop[1].bctype == Opm::BCType::THERMAL);
BOOST_CHECK(prop[1].component == Opm::BCComponent::WATER);
using measure = Opm::UnitSystem::measure;
BOOST_CHECK_EQUAL(deck.getActiveUnitSystem().to_si("Mass/Time*Length*Length", -0.05),
prop[1].rate);
BOOST_CHECK(!prop[1].pressure.has_value());
BOOST_CHECK_EQUAL(deck.getActiveUnitSystem().to_si(measure::temperature, 50.0),
*prop[1].temperature);
BOOST_CHECK(prop[2].bctype == Opm::BCType::DIRICHLET);
BOOST_CHECK(prop[2].component == Opm::BCComponent::OIL);
BOOST_CHECK_EQUAL(deck.getActiveUnitSystem().to_si("Mass/Time*Length*Length", -0.01),
prop[2].rate);
BOOST_CHECK_EQUAL(deck.getActiveUnitSystem().to_si(measure::pressure, 279.0),
*prop[2].pressure);
}