// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
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 2 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 .
Consult the COPYING file in the top-level source directory of this
module for the precise wording of the license and the list of
copyright holders.
*/
/*!
* \file
*
* \brief This is the unit test for the black oil PVT classes
*
* This test requires the presence of opm-parser.
*/
#include "config.h"
#if !HAVE_ECL_INPUT
#error "The test for the black oil PVT classes requires eclipse input support in opm-common"
#endif
#include
#define BOOST_TEST_MODULE EclBlackOilPvt
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
// values of strings based on the first SPE1 test case of opm-data. note that in the
// real world it does not make much sense to specify a fluid phase using more than a
// single keyword, but for a unit test, this saves a lot of boiler-plate code.
static constexpr const char* deckString1 =
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
" * 2 /\n"
"\n"
"OIL\n"
"GAS\n"
"WATER\n"
"\n"
"DISGAS\n"
"\n"
"METRIC\n"
"\n"
"GRID\n"
"\n"
"DX\n"
" 300*1000 /\n"
"DY\n"
" 300*1000 /\n"
"DZ\n"
" 100*20 100*30 100*50 /\n"
"\n"
"TOPS\n"
" 100*1234 /\n"
"\n"
"PORO\n"
" 300*0.15 /\n"
"PROPS\n"
"\n"
"DENSITY\n"
" 859.5 1033.0 0.854 /\n"
" 860.04 1033.0 0.853 /\n"
"\n"
"PVTW\n"
" 1.0 1.1 1e-6 1.1 2.0e-9 /\n"
" 2.0 1.2 1e-7 1.2 3.0e-9 /\n"
"\n"
"PVCDO\n"
" 1.0 1.1 1e-6 1.1 2.0e-9 /\n"
" 2.0 1.2 1e-7 1.2 3.0e-9 /\n"
"\n"
"PVDG\n"
"1.0 1.0 10.0\n"
"2.0 * *\n"
"3.0 1e-10 30.0 /\n"
"\n"
"4.0 1.0 40.0\n"
"5.0 * *\n"
"6.0 1e-10 60.0 /\n"
"\n"
"PVTG\n"
"\n"
"-- PVT region 2 --\n"
"-- PRESSURE RV BG VISCOSITY\n"
" 1.00 1.1e-3 1.1 0.01\n"
" 1.0e-3 1.15 0.005 /\n"
"\n"
" 500.00 0.9e-3 1.2 0.02\n"
" 0.8e-3 1.25 0.015 /\n"
"/\n"
"\n"
"-- PVT region 2 --\n"
"-- PRESSURE RV BG VISCOSITY\n"
" 2.00 2.1e-3 2.1 0.02\n"
" 2.0e-3 2.15 0.015 /\n"
"\n"
" 502.00 1.2e-3 2.2 2.02\n"
" 1.1e-3 2.25 2.015 /\n"
"/\n"
"\n";
template
void ensurePvtApi(const OilPvt& oilPvt, const GasPvt& gasPvt, const WaterPvt& waterPvt)
{
// we don't want to run this, we just want to make sure that it compiles
while (0) {
Evaluation temperature = 273.15 + 20.0;
Evaluation pressure = 1e5;
Evaluation saltconcentration = 0.0;
Evaluation Rs = 0.0;
Evaluation Rsw = 0.0;
Evaluation Rv = 0.0;
Evaluation Rvw = 0.0;
Evaluation So = 0.5;
Evaluation maxSo = 1.0;
[[maybe_unused]] Evaluation tmp;
/////
// water PVT API
/////
tmp = waterPvt.viscosity(/*regionIdx=*/0,
temperature,
pressure,
Rsw,
saltconcentration);
tmp = waterPvt.inverseFormationVolumeFactor(/*regionIdx=*/0,
temperature,
pressure,
Rsw,
saltconcentration);
/////
// oil PVT API
/////
tmp = oilPvt.viscosity(/*regionIdx=*/0,
temperature,
pressure,
Rs);
tmp = oilPvt.inverseFormationVolumeFactor(/*regionIdx=*/0,
temperature,
pressure,
Rs);
tmp = oilPvt.saturatedViscosity(/*regionIdx=*/0,
temperature,
pressure);
tmp = oilPvt.saturatedInverseFormationVolumeFactor(/*regionIdx=*/0,
temperature,
pressure);
tmp = oilPvt.saturationPressure(/*regionIdx=*/0,
temperature,
Rs);
tmp = oilPvt.saturatedGasDissolutionFactor(/*regionIdx=*/0,
temperature,
pressure);
tmp = oilPvt.saturatedGasDissolutionFactor(/*regionIdx=*/0,
temperature,
pressure,
So,
maxSo);
/////
// gas PVT API
/////
tmp = gasPvt.viscosity(/*regionIdx=*/0,
temperature,
pressure,
Rv,
Rvw);
tmp = gasPvt.inverseFormationVolumeFactor(/*regionIdx=*/0,
temperature,
pressure,
Rv,
Rvw);
tmp = gasPvt.saturatedViscosity(/*regionIdx=*/0,
temperature,
pressure);
tmp = gasPvt.saturatedInverseFormationVolumeFactor(/*regionIdx=*/0,
temperature,
pressure);
tmp = gasPvt.saturationPressure(/*regionIdx=*/0,
temperature,
Rv);
tmp = gasPvt.saturatedOilVaporizationFactor(/*regionIdx=*/0,
temperature,
pressure);
tmp = gasPvt.saturatedOilVaporizationFactor(/*regionIdx=*/0,
temperature,
pressure,
So,
maxSo);
// prevent GCC from producing a "variable assigned but unused" warning
tmp = 2.0*tmp;
}
}
struct Fixture {
Fixture()
: python(std::make_shared())
, deck(Opm::Parser().parseString(deckString1))
, eclState(deck)
, schedule(deck, eclState, python)
{
}
std::shared_ptr python;
Opm::Deck deck;
Opm::EclipseState eclState;
Opm::Schedule schedule;
};
BOOST_FIXTURE_TEST_SUITE(Generic, Fixture)
using Types = boost::mpl::list;
BOOST_AUTO_TEST_CASE_TEMPLATE(ApiConformance, Scalar, Types)
{
Opm::GasPvtMultiplexer gasPvt;
Opm::OilPvtMultiplexer oilPvt;
Opm::WaterPvtMultiplexer waterPvt;
gasPvt.initFromState(eclState, schedule);
oilPvt.initFromState(eclState, schedule);
waterPvt.initFromState(eclState, schedule);
using FooEval = Opm::DenseAd::Evaluation;
ensurePvtApi(oilPvt, gasPvt, waterPvt);
ensurePvtApi(oilPvt, gasPvt, waterPvt);
}
BOOST_AUTO_TEST_CASE_TEMPLATE(ConstantCompressibilityWater, Scalar, Types)
{
constexpr Scalar tolerance = std::numeric_limits::epsilon()*1e3;
Opm::ConstantCompressibilityWaterPvt constCompWaterPvt;
constCompWaterPvt.initFromState(eclState, schedule);
// make sure that the values at the reference points are the ones specified in the
// deck.
Scalar refTmp, tmp;
refTmp = 1.1e-3; // the deck value is given in cP, while the SI units use Pa s...
tmp = constCompWaterPvt.viscosity(/*regionIdx=*/0,
/*temperature=*/273.15 + 20.0,
/*pressure=*/1e5,
/*disgas_in_water*/0.0,
/*saltconcentration=*/0.0);
BOOST_CHECK_MESSAGE(std::abs(tmp - refTmp) <= tolerance,
"The reference water viscosity at region 0 is supposed to be " <<
refTmp << ". (is " << tmp << ")");
refTmp = 1.2e-3;
tmp = constCompWaterPvt.viscosity(/*regionIdx=*/1,
/*temperature=*/273.15 + 20.0,
/*pressure=*/2e5,
/*disgas_in_water*/0.0,
/*saltconcentration=*/0.0);
BOOST_CHECK_MESSAGE(std::abs(tmp - refTmp) <= tolerance,
"The reference water viscosity at region 1 is supposed to be " <<
refTmp << ". (is " << tmp << ")");
}
BOOST_AUTO_TEST_SUITE_END()