OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity
master
opm-common/tests/material/test_eclmateriallawmanager.cpp

942 lines
30 KiB
C++
Raw Permalink Normal View History

various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
// -*- 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 <http://www.gnu.org/licenses/>.
clean up the licensing preable of source files this patch removes the in-file lists in favor of a global list of in the COPYING file. this is done because (a) maintaining a list of authors at the beginning of each source file is a major pain in the a**, (b) for this reason, the list of authors was not accurate in about 85% of all cases where more than one person was involved and (c) this list is not legally binding in any way (the copyright is at the person who authored a given change; if these lists had any legal relevance, one could "aquire" the copyright of the module by forking it and replacing the lists...)
2016-03-14 08:11:22 -05:00
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.
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
*/
/*!
* \file
*
* \brief This is the unit test for the class which manages the parameters for the ECL
* saturation functions.
*
* This test requires the presence of opm-parser.
*/
#include "config.h"
adjustments for moved eclipse support (now in opm-common)
2018-03-05 05:48:09 -06:00
#if !HAVE_ECL_INPUT
#error "The test for EclMaterialLawManager requires eclipse input support in opm-common"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
#endif
fixed: earlier boosts only accepts boost::mpl::list for BOOST_AUTO_TEST_TEMPLATE
2023-11-10 03:12:13 -06:00
#include <boost/mpl/list.hpp>
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
#define BOOST_TEST_MODULE EclMaterialLawManager
#include <boost/test/unit_test.hpp>
EclEpsScalingPoints(Info): put init and print in separate compile unit this allows to encapsulate EclipseState usage and iostream usage
2022-12-09 08:28:25 -06:00
#include <opm/material/fluidmatrixinteractions/EclEpsGridProperties.hpp>
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
#include <opm/material/fluidmatrixinteractions/EclMaterialLawManager.hpp>
#include <opm/material/fluidstates/SimpleModularFluidState.hpp>
Downstream of opm-common filesystem reorganisation
2021-12-14 01:24:17 -06:00
#include <opm/input/eclipse/Parser/Parser.hpp>
#include <opm/input/eclipse/Deck/Deck.hpp>
#include <opm/input/eclipse/EclipseState/EclipseState.hpp>
#include <opm/input/eclipse/EclipseState/Grid/EclipseGrid.hpp>
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
// values of strings taken from the SPE1 test case1 of opm-data
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam1DeckString =
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"WATER\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\n"
"\n"
"GRID\n"
"\n"
"DX\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 300*1000 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"DY\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 300*1000 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"DZ\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 100*20 100*30 100*50 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"TOPS\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 100*8325 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"\n"
Add PORO to testdata
2019-11-29 04:17:40 -06:00
"PORO\n"
" 300*0.15 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"PROPS\n"
"\n"
"SWOF\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"0.12 0 1 0\n"
"0.18 4.64876033057851E-008 1 0\n"
"0.24 0.000000186 0.997 0\n"
"0.3 4.18388429752066E-007 0.98 0\n"
"0.36 7.43801652892562E-007 0.7 0\n"
"0.42 1.16219008264463E-006 0.35 0\n"
"0.48 1.67355371900826E-006 0.2 0\n"
"0.54 2.27789256198347E-006 0.09 0\n"
"0.6 2.97520661157025E-006 0.021 0\n"
"0.66 3.7654958677686E-006 0.01 0\n"
"0.72 4.64876033057851E-006 0.001 0\n"
"0.78 0.000005625 0.0001 0\n"
"0.84 6.69421487603306E-006 0 0\n"
"0.91 8.05914256198347E-006 0 0\n"
"1 0.984 0 0 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"\n"
"SGOF\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"0 0 1 0\n"
"0.001 0 1 0\n"
"0.02 0 0.997 0\n"
"0.05 0.005 0.980 0\n"
"0.12 0.025 0.700 0\n"
"0.2 0.075 0.350 0\n"
"0.25 0.125 0.200 0\n"
"0.3 0.190 0.090 0\n"
"0.4 0.410 0.021 0\n"
"0.45 0.60 0.010 0\n"
"0.5 0.72 0.001 0\n"
"0.6 0.87 0.0001 0\n"
"0.7 0.94 0.000 0\n"
"0.85 0.98 0.000 0\n"
"0.88 0.984 0.000 0 /\n";
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam2DeckString =
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"WATER\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\n"
"\n"
"GRID\n"
"\n"
"DX\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 300*1000 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"DY\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 300*1000 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"DZ\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 100*20 100*30 100*50 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"TOPS\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 100*8325 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"\n"
Add PORO to testdata
2019-11-29 04:17:40 -06:00
"PORO\n"
" 300*0.15 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"PROPS\n"
"\n"
"PVTW\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 4017.55 1.038 3.22E-6 0.318 0.0 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"\n"
"SWFN\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"0.12 0 0\n"
"0.18 4.64876033057851E-008 0\n"
"0.24 0.000000186 0\n"
"0.3 4.18388429752066E-007 0\n"
"0.36 7.43801652892562E-007 0\n"
"0.42 1.16219008264463E-006 0\n"
"0.48 1.67355371900826E-006 0\n"
"0.54 2.27789256198347E-006 0\n"
"0.6 2.97520661157025E-006 0\n"
"0.66 3.7654958677686E-006 0\n"
"0.72 4.64876033057851E-006 0\n"
"0.78 0.000005625 0\n"
"0.84 6.69421487603306E-006 0\n"
"0.91 8.05914256198347E-006 0\n"
"1 0.984 0 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"\n"
"SGFN\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"0 0 0\n"
"0.001 0 0\n"
"0.02 0 0\n"
"0.05 0.005 0\n"
"0.12 0.025 0\n"
"0.2 0.075 0\n"
"0.25 0.125 0\n"
"0.3 0.190 0\n"
"0.4 0.410 0\n"
"0.45 0.60 0\n"
"0.5 0.72 0\n"
"0.6 0.87 0\n"
"0.7 0.94 0\n"
"0.85 0.98 0\n"
"0.88 0.984 0 /\n"
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
"\n"
"SOF3\n"
Fix the SOF3 input to make the KROW input similar as in SWOF.
2015-10-12 07:26:19 -05:00
" 0 0 0 \n"
" 0.03 0 0 \n"
" 0.09 0 0 \n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
" 0.16 0 0 \n"
Fix the SOF3 input to make the KROW input similar as in SWOF.
2015-10-12 07:26:19 -05:00
" 0.18 1* 0 \n"
" 0.22 0.0001 1* \n"
" 0.28 0.001 0.0001 \n"
" 0.34 0.01 1* \n"
" 0.38 1* 0.001 \n"
" 0.40 0.021 1* \n"
" 0.43 1* 0.01 \n"
" 0.46 0.09 1* \n"
" 0.48 1* 0.021 \n"
" 0.52 0.2 1* \n"
" 0.58 0.35 0.09 \n"
" 0.63 1* 0.2 \n"
" 0.64 0.7 1* \n"
" 0.68 1* 0.35 \n"
" 0.70 0.98 1* \n"
" 0.76 0.997 0.7 \n"
" 0.83 1 0.98 \n"
" 0.86 1 0.997 \n"
" 0.879 1 1 \n"
" 0.88 1 1 / \n"
"\n";
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
//Taken as a mix of the SPE1 cases above, and Norne to enable hysteresis
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* hysterDeckString =
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"WATER\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
Add PORO to testdata
2019-11-29 04:17:40 -06:00
"PORO\n"
" 300*0.15 /\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"\n"
"\n"
"EHYSTR\n"
Activate cap pressure hysteresis.
2022-10-10 11:16:59 -05:00
"0.1 0 0.1 1* BOTH /\n"
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
"\n"
"SATOPTS\n"
"HYSTER /\n"
"\n"
"PROPS\n"
"\n"
"SWOF\n"
"0.12 0 1 0\n"
"0.18 4.64876033057851E-008 1 0\n"
"0.24 0.000000186 0.997 0\n"
"0.3 4.18388429752066E-007 0.98 0\n"
"0.36 7.43801652892562E-007 0.7 0\n"
"0.42 1.16219008264463E-006 0.35 0\n"
"0.48 1.67355371900826E-006 0.2 0\n"
"0.54 2.27789256198347E-006 0.09 0\n"
"0.6 2.97520661157025E-006 0.021 0\n"
"0.66 3.7654958677686E-006 0.01 0\n"
"0.72 4.64876033057851E-006 0.001 0\n"
"0.78 0.000005625 0.0001 0\n"
"0.84 6.69421487603306E-006 0 0\n"
"0.91 8.05914256198347E-006 0 0\n"
"1 0.984 0 0 /\n"
"\n"
"\n"
"SGOF\n"
"0 0 1 0\n"
"0.001 0 1 0\n"
"0.02 0 0.997 0\n"
"0.05 0.005 0.980 0\n"
"0.12 0.025 0.700 0\n"
"0.2 0.075 0.350 0\n"
"0.25 0.125 0.200 0\n"
"0.3 0.190 0.090 0\n"
"0.4 0.410 0.021 0\n"
"0.45 0.60 0.010 0\n"
"0.5 0.72 0.001 0\n"
"0.6 0.87 0.0001 0\n"
"0.7 0.94 0.000 0\n"
"0.85 0.98 0.000 0\n"
"0.88 0.984 0.000 0 /\n";
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam1DeckStringGasOil =
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
"\n"
"\n"
Add PORO to testdata
2019-11-29 04:17:40 -06:00
"PORO\n"
" 300*0.15 /\n"
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
"PROPS\n"
"\n"
"\n"
"SGOF\n"
"0 0 1 0\n"
"0.001 0 1 0\n"
"0.02 0 0.997 0\n"
"0.05 0.005 0.980 0\n"
"0.12 0.025 0.700 0\n"
"0.2 0.075 0.350 0\n"
"0.25 0.125 0.200 0\n"
"0.3 0.190 0.090 0\n"
"0.4 0.410 0.021 0\n"
"0.45 0.60 0.010 0\n"
"0.5 0.72 0.001 0\n"
"0.6 0.87 0.0001 0\n"
"0.7 0.94 0.000 0\n"
"0.85 0.98 0.000 0\n"
"0.88 0.984 0.000 0 /\n";
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam2DeckStringGasOil =
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
Add PORO to testdata
2019-11-29 04:17:40 -06:00
"PORO\n"
" 300*0.15 /\n"
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
"\n"
"\n"
"PROPS\n"
"\n"
"PVTW\n"
" 4017.55 1.038 3.22E-6 0.318 0.0 /\n"
"\n"
"\n"
"SGFN\n"
"0 0 0\n"
"0.001 0 0\n"
"0.02 0 0\n"
"0.05 0.005 0\n"
"0.12 0.025 0\n"
"0.2 0.075 0\n"
"0.25 0.125 0\n"
"0.3 0.190 0\n"
"0.4 0.410 0\n"
"0.45 0.60 0\n"
"0.5 0.72 0\n"
"0.6 0.87 0\n"
"0.7 0.94 0\n"
"0.85 0.98 0\n"
"0.88 0.984 0 /\n"
"\n"
"SOF2\n"
"0.12 0.000 \n"
"0.15 0.000 \n"
"0.3 0.000 \n"
"0.4 0.0001 \n"
"0.5 0.001 \n"
"0.55 0.010 \n"
"0.6 0.021 \n"
"0.7 0.090 \n"
"0.8 0.350 \n"
"0.88 0.700 \n"
"0.95 0.980 \n"
"0.98 0.997 \n"
"0.999 1 \n"
"1.0 1 \n /\n";
Add option to store/restore hysteresis parameters
2017-02-24 10:09:40 -06:00
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* letDeckString =
Some LET related tests.
2022-01-28 11:02:31 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"OIL\n"
"GAS\n"
"WATER\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
"\n"
"\n"
"PORO\n"
" 300*0.15 /\n"
"PROPS\n"
"\n"
"SWOFLET\n"
" 0.1 0.2 1.5 7.0 1.5 0.5 0.05 0.15 3.5 3.0 1.3 1.0 0.7 17.0 0.95 3.8 0.04 /\n"
"\n"
"SGOFLET\n"
" 0.0 0.03 1.8 1.9 1.0 0.95 0.0 0.01 3.5 4.0 1.1 1.0 1.0 1.0 1.0 0.2 0.01 /\n"
"\n";
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam3DeckStringGasWater =
Implement WSF/GSF satfunc familiy for CO2STORE
2023-02-21 03:09:36 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"WATER\n"
"GAS\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
"\n"
"\n"
"PORO\n"
" 300*0.15 /\n"
"PROPS\n"
"\n"
"\n"
"GSF\n"
"0 0 0\n"
"0.001 0 0\n"
"0.02 0 0\n"
"0.05 0.005 0\n"
"0.12 0.025 0\n"
"0.2 0.075 0\n"
"0.25 0.125 0\n"
"0.3 0.190 0\n"
"0.4 0.410 0\n"
"0.45 0.60 0\n"
"0.5 0.72 0\n"
"0.6 0.87 0\n"
"0.7 0.94 0\n"
"0.85 0.98 0\n"
"0.88 0.984 0 /\n"
"\n"
"WSF\n"
"0.12 0 \n"
"0.22 0 \n"
"0.55 0.005 \n"
"0.88 0.984 /\n";
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
static constexpr const char* fam2DeckStringGasWater =
Implement WSF/GSF satfunc familiy for CO2STORE
2023-02-21 03:09:36 -06:00
"RUNSPEC\n"
"\n"
"DIMENS\n"
" 10 10 3 /\n"
"\n"
"TABDIMS\n"
"/\n"
"\n"
"WATER\n"
"GAS\n"
"\n"
"DISGAS\n"
"\n"
"FIELD\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*8325 /\n"
"\n"
"\n"
"PORO\n"
" 300*0.15 /\n"
"PROPS\n"
"\n"
"\n"
"SGFN\n"
"0 0 0\n"
"0.001 0 0\n"
"0.02 0 0\n"
"0.05 0.005 0\n"
"0.12 0.025 0\n"
"0.2 0.075 0\n"
"0.25 0.125 0\n"
"0.3 0.190 0\n"
"0.4 0.410 0\n"
"0.45 0.60 0\n"
"0.5 0.72 0\n"
"0.6 0.87 0\n"
"0.7 0.94 0\n"
"0.85 0.98 0\n"
"0.88 0.984 0 /\n"
"\n"
"SWFN\n"
"0.12 0 0\n"
"0.22 0 0\n"
"0.55 0.005 0\n"
"0.88 0.984 0 /\n";
Some LET related tests.
2022-01-28 11:02:31 -06:00
template <class Scalar>
inline Scalar computeLetCurve(const Scalar S, const Scalar L, const Scalar E, const Scalar T)
{
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
if (S <= 0.0) {
return 0.0;
} else if (S >= 1.0) {
return 1.0;
}
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
const Scalar powS = Opm::pow(S,L);
const Scalar pow1mS = Opm::pow(1.0-S,T);
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
return powS/(powS+pow1mS*E);
Some LET related tests.
2022-01-28 11:02:31 -06:00
}
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
template<class Scalar>
struct Fixture {
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
enum { numPhases = 3 };
enum { waterPhaseIdx = 0 };
enum { oilPhaseIdx = 1 };
enum { gasPhaseIdx = 2 };
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
using MaterialTraits = Opm::ThreePhaseMaterialTraits<Scalar,
waterPhaseIdx,
oilPhaseIdx,
gasPhaseIdx>;
using FluidState = Opm::SimpleModularFluidState<Scalar,
/*numPhases=*/3,
/*numComponents=*/3,
void,
/*storePressure=*/false,
/*storeTemperature=*/false,
/*storeComposition=*/false,
/*storeFugacity=*/false,
/*storeSaturation=*/true,
/*storeDensity=*/false,
/*storeViscosity=*/false,
/*storeEnthalpy=*/false>;
using MaterialLawManager = Opm::EclMaterialLawManager<MaterialTraits>;
using MaterialLaw = typename MaterialLawManager::MaterialLaw;
};
fixed: earlier boosts only accepts boost::mpl::list for BOOST_AUTO_TEST_TEMPLATE
2023-11-10 03:12:13 -06:00
using Types = boost::mpl::list<float,double>;
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
BOOST_AUTO_TEST_CASE_TEMPLATE(Fam1Fam2Hysteresis, Scalar, Types)
{
using MaterialLaw = typename Fixture<Scalar>::MaterialLaw;
using MaterialLawManager = typename Fixture<Scalar>::MaterialLawManager;
constexpr int numPhases = Fixture<Scalar>::numPhases;
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
Opm::Parser parser;
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
const auto deck = parser.parseString(fam1DeckString);
const Opm::EclipseState eclState(deck);
const auto& eclGrid = eclState.getInputGrid();
size_t n = eclGrid.getCartesianSize();
MaterialLawManager fam1materialLawManager;
fam1materialLawManager.initFromState(eclState);
fam1materialLawManager.initParamsForElements(eclState, n);
BOOST_CHECK(!fam1materialLawManager.enableEndPointScaling());
BOOST_CHECK(!fam1materialLawManager.enableHysteresis());
const auto fam2Deck = parser.parseString(fam2DeckString);
const Opm::EclipseState fam2EclState(fam2Deck);
MaterialLawManager fam2MaterialLawManager;
fam2MaterialLawManager.initFromState(fam2EclState);
fam2MaterialLawManager.initParamsForElements(fam2EclState, n);
BOOST_CHECK(!fam2MaterialLawManager.enableEndPointScaling());
BOOST_CHECK(!fam2MaterialLawManager.enableHysteresis());
const auto hysterDeck = parser.parseString(hysterDeckString);
const Opm::EclipseState hysterEclState(hysterDeck);
MaterialLawManager hysterMaterialLawManager;
hysterMaterialLawManager.initFromState(hysterEclState);
hysterMaterialLawManager.initParamsForElements(hysterEclState, n);
BOOST_CHECK(!hysterMaterialLawManager.enableEndPointScaling());
BOOST_CHECK(hysterMaterialLawManager.enableHysteresis());
// make sure that the saturation functions for both keyword families are
// identical, and that setting and getting the hysteresis parameters works
for (unsigned elemIdx = 0; elemIdx < n; ++elemIdx) {
for (int i = -10; i < 120; ++i) {
Scalar Sw = Scalar(i) / 100;
for (int j = i; j < 120; ++j) {
Scalar So = Scalar(j) / 100;
Scalar Sg = 1 - Sw - So;
typename Fixture<Scalar>::FluidState fs;
fs.setSaturation(Fixture<Scalar>::waterPhaseIdx, Sw);
fs.setSaturation(Fixture<Scalar>::oilPhaseIdx, So);
fs.setSaturation(Fixture<Scalar>::gasPhaseIdx, Sg);
std::array<Scalar,numPhases> pcFam1 = {0.0, 0.0};
std::array<Scalar,numPhases> pcFam2 = {0.0, 0.0};
MaterialLaw::capillaryPressures(pcFam1,
fam1materialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::capillaryPressures(pcFam2,
fam2MaterialLawManager.materialLawParams(elemIdx),
fs);
std::array<Scalar,numPhases> krFam1 = {0.0, 0.0};
std::array<Scalar,numPhases> krFam2 = {0.0, 0.0};
MaterialLaw::relativePermeabilities(krFam1,
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
fam1materialLawManager.materialLawParams(elemIdx),
fs);
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
MaterialLaw::relativePermeabilities(krFam2,
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
fam2MaterialLawManager.materialLawParams(elemIdx),
fs);
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
BOOST_CHECK_MESSAGE(std::abs(pcFam1[phaseIdx] - pcFam2[phaseIdx]) <= 1e-5,
"Discrepancy between capillary pressure of family 1 and family 2 keywords");
BOOST_CHECK_MESSAGE(std::abs(krFam1[phaseIdx] - krFam2[phaseIdx]) <= 1e-1,
"Discrepancy between relative permeabilities of family 1 and family 2 keywords");
Add support for SOF2 for 2p runs. Make it possible to use SOF2 in combination with either SGFN og SWFN for gas-oil and water-oil case respectivly.
2017-11-23 03:46:55 -06:00
}
Implement WSF/GSF satfunc familiy for CO2STORE
2023-02-21 03:09:36 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
// This should ideally test each of the materials (stone1, stone2, default, two-phase),
// but currently only tests default
const std::array<Scalar,2> pcSwMdc_in = {1.0 / 2.0, 1.0 / 3.0};
const std::array<Scalar,2> krnSwMdc_in = {1.0 / 5.0, 1.0 / 7.0};
hysterMaterialLawManager.setOilWaterHysteresisParams(pcSwMdc_in[0],
krnSwMdc_in[0],
elemIdx);
hysterMaterialLawManager.setGasOilHysteresisParams(pcSwMdc_in[1],
krnSwMdc_in[1],
elemIdx);
std::array<Scalar,2> pcSwMdc_out = {0.0, 0.0};
std::array<Scalar,2> krnSwMdc_out = {0.0, 0.0};
hysterMaterialLawManager.oilWaterHysteresisParams(pcSwMdc_out[0],
krnSwMdc_out[0],
elemIdx);
hysterMaterialLawManager.gasOilHysteresisParams(pcSwMdc_out[1],
krnSwMdc_out[1],
elemIdx);
for (unsigned phasePairIdx = 0; phasePairIdx < 2; ++phasePairIdx) {
BOOST_CHECK_MESSAGE(pcSwMdc_in[phasePairIdx] == pcSwMdc_out[phasePairIdx],
"Hysteresis parameters did not propagate correctly");
BOOST_CHECK_MESSAGE(krnSwMdc_in[phasePairIdx] == krnSwMdc_out[phasePairIdx],
"Hysteresis parameters did not propagate correctly");
Implement WSF/GSF satfunc familiy for CO2STORE
2023-02-21 03:09:36 -06:00
}
}
}
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
}
}
Implement WSF/GSF satfunc familiy for CO2STORE
2023-02-21 03:09:36 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
BOOST_AUTO_TEST_CASE_TEMPLATE(GasOil, Scalar, Types)
{
using MaterialLaw = typename Fixture<Scalar>::MaterialLaw;
using MaterialLawManager = typename Fixture<Scalar>::MaterialLawManager;
constexpr int numPhases = Fixture<Scalar>::numPhases;
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
Opm::Parser parser;
const auto fam1Deck = parser.parseString(fam1DeckStringGasOil);
const Opm::EclipseState fam1EclState(fam1Deck);
const auto& eclGrid = fam1EclState.getInputGrid();
const size_t n = eclGrid.getCartesianSize();
MaterialLawManager fam1materialLawManager;
fam1materialLawManager.initFromState(fam1EclState);
fam1materialLawManager.initParamsForElements(fam1EclState, n);
const auto fam2Deck = parser.parseString(fam2DeckStringGasOil);
const Opm::EclipseState fam2EclState(fam2Deck);
MaterialLawManager fam2MaterialLawManager;
fam2MaterialLawManager.initFromState(fam2EclState);
fam2MaterialLawManager.initParamsForElements(fam2EclState, n);
for (unsigned elemIdx = 0; elemIdx < n; ++elemIdx) {
for (int i = 0; i < 100; ++ i) {
Scalar Sw = 0;
Scalar So = Scalar(i) / 100;
Scalar Sg = 1 - Sw - So;
typename Fixture<Scalar>::FluidState fs;
fs.setSaturation(Fixture<Scalar>::waterPhaseIdx, Sw);
fs.setSaturation(Fixture<Scalar>::oilPhaseIdx, So);
fs.setSaturation(Fixture<Scalar>::gasPhaseIdx, Sg);
std::array<Scalar,numPhases> pcFam1 = {0.0, 0.0};
std::array<Scalar,numPhases> pcFam2 = {0.0, 0.0};
MaterialLaw::capillaryPressures(pcFam1,
fam1materialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::capillaryPressures(pcFam2,
fam2MaterialLawManager.materialLawParams(elemIdx),
fs);
std::array<Scalar,numPhases> krFam1 = {0.0, 0.0};
std::array<Scalar,numPhases> krFam2 = {0.0, 0.0};
MaterialLaw::relativePermeabilities(krFam1,
fam1materialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::relativePermeabilities(krFam2,
fam2MaterialLawManager.materialLawParams(elemIdx),
fs);
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
BOOST_CHECK_MESSAGE(std::abs(pcFam1[phaseIdx] - pcFam2[phaseIdx]) <= 1e-5,
"Discrepancy between capillary pressure of family 1 and family 2 keywords");
BOOST_CHECK_MESSAGE(krFam1[phaseIdx] - krFam2[phaseIdx] <= 1e-1,
"Discrepancy between relative permeabilities of family 1 and family 2 keywords");
}
}
}
}
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
BOOST_AUTO_TEST_CASE_TEMPLATE(GasWater, Scalar, Types)
{
using MaterialLaw = typename Fixture<Scalar>::MaterialLaw;
using MaterialLawManager = typename Fixture<Scalar>::MaterialLawManager;
constexpr int numPhases = Fixture<Scalar>::numPhases;
Some LET related tests.
2022-01-28 11:02:31 -06:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
Opm::Parser parser;
const auto fam2Deck = parser.parseString(fam2DeckStringGasWater);
const Opm::EclipseState fam2EclState(fam2Deck);
const auto& eclGrid = fam2EclState.getInputGrid();
const size_t n = eclGrid.getCartesianSize();
MaterialLawManager fam2materialLawManager;
fam2materialLawManager.initFromState(fam2EclState);
fam2materialLawManager.initParamsForElements(fam2EclState, n);
const auto fam3Deck = parser.parseString(fam3DeckStringGasWater);
const Opm::EclipseState fam3EclState(fam3Deck);
MaterialLawManager fam3materialLawManager;
fam3materialLawManager.initFromState(fam3EclState);
fam3materialLawManager.initParamsForElements(fam3EclState, n);
for (unsigned elemIdx = 0; elemIdx < n; ++elemIdx) {
for (int i = 0; i < 100; ++ i) {
Scalar Sw = 0;
Scalar So = Scalar(i) / 100;
Scalar Sg = 1 - Sw - So;
typename Fixture<Scalar>::FluidState fs;
fs.setSaturation(Fixture<Scalar>::waterPhaseIdx, Sw);
fs.setSaturation(Fixture<Scalar>::oilPhaseIdx, So);
fs.setSaturation(Fixture<Scalar>::gasPhaseIdx, Sg);
std::array<Scalar,numPhases> pcFam2 = {0.0, 0.0};
std::array<Scalar,numPhases> pcFam3 = {0.0, 0.0};
MaterialLaw::capillaryPressures(pcFam2,
fam2materialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::capillaryPressures(pcFam3,
fam3materialLawManager.materialLawParams(elemIdx),
fs);
std::array<Scalar,numPhases> krFam2 = {0.0, 0.0};
std::array<Scalar,numPhases> krFam3 = {0.0, 0.0};
MaterialLaw::relativePermeabilities(krFam2,
fam2materialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::relativePermeabilities(krFam3,
fam3materialLawManager.materialLawParams(elemIdx),
fs);
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
BOOST_CHECK_MESSAGE(std::abs(pcFam2[phaseIdx] - pcFam3[phaseIdx]) <= 1e-5,
"Discrepancy between capillary pressure of family 2 and family 3 keywords");
BOOST_CHECK_MESSAGE(std::abs(krFam2[phaseIdx] - krFam3[phaseIdx]) <= 1e-1,
"Discrepancy between relative permeabilities of family 2 and family 3 keywords");
Some LET related tests.
2022-01-28 11:02:31 -06:00
}
}
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
}
tests: added checks with Scalar=float. All tests compile with Scalar=float, but only some work.
2016-01-14 18:28:07 -06:00
}
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
BOOST_AUTO_TEST_CASE_TEMPLATE(Let, Scalar, Types)
tests: added checks with Scalar=float. All tests compile with Scalar=float, but only some work.
2016-01-14 18:28:07 -06:00
{
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
using MaterialLaw = typename Fixture<Scalar>::MaterialLaw;
using MaterialLawManager = typename Fixture<Scalar>::MaterialLawManager;
constexpr int numPhases = Fixture<Scalar>::numPhases;
constexpr int gasPhaseIdx = Fixture<Scalar>::gasPhaseIdx;
constexpr int oilPhaseIdx = Fixture<Scalar>::oilPhaseIdx;
constexpr int waterPhaseIdx = Fixture<Scalar>::waterPhaseIdx;
initialize MPI for all unit tests This started to cause crashes for me. I don't really have an idea why, but it is not too bad of a thing to do anyway...
2016-04-17 04:28:06 -05:00
test_eclmateriallawmanager: convert to boost::test
2023-05-25 15:53:25 -05:00
Opm::Parser parser;
const auto letDeck = parser.parseString(letDeckString);
const Opm::EclipseState letEclState(letDeck);
const auto& eclGrid = letEclState.getInputGrid();
const size_t n = eclGrid.getCartesianSize();
MaterialLawManager letmaterialLawManager;
letmaterialLawManager.initFromState(letEclState);
letmaterialLawManager.initParamsForElements(letEclState, n);
Scalar Swco = 0.1;
Scalar psi2Pa = 6894.7573;
for (unsigned elemIdx = 0; elemIdx < 1; ++elemIdx) {
for (int i = -10; i < 120; ++i) {
Scalar So = Scalar(i) / 100;
// Oil in gas and conate water
Scalar Sw = Swco;
Scalar Sg = 1 - Sw - So;
typename Fixture<Scalar>::FluidState fs;
fs.setSaturation(waterPhaseIdx, Sw);
fs.setSaturation(oilPhaseIdx, So);
fs.setSaturation(gasPhaseIdx, Sg);
std::array<Scalar,numPhases> pcLET = {0.0, 0.0, 0.0};
MaterialLaw::capillaryPressures(pcLET,
letmaterialLawManager.materialLawParams(elemIdx),
fs);
std::array<Scalar,numPhases> krLET = {0.0, 0.0, 0.0};
MaterialLaw::relativePermeabilities(krLET,
letmaterialLawManager.materialLawParams(elemIdx),
fs);
//SGOFLET
//0.0 0.03 1.8 1.9 1.0 0.95 0.0 0.01 3.5 4.0 1.1 1.0 1.0 1.0 1.0 0.2 0.01
Scalar krg_let = 0.95 * computeLetCurve((1.0 - So - 0.03 - Swco) / (1.0 - 0.03 - 0.01 - Swco), 1.8, 1.9, 1.0);
Scalar krog_let = 1.0 * computeLetCurve((So - 0.01) / (1.0 - 0.03 - 0.01 - Swco), 3.5, 4.0, 1.1);
//S=(So-Sorg)/(1-Sorg-Sgl-Swco), Pc = Pct + (pcir_pc-Pct)*(1-S)^L/[(1-S)^L+E*S^T]
Scalar pcog_let = psi2Pa*(0.01 + (0.2 - 0.01) * computeLetCurve((1.0 - So - Swco) / (1.0 - Swco), 1.0, 1.0, 1.0));
BOOST_CHECK_MESSAGE(std::abs(krLET[gasPhaseIdx] - krg_let) <= 1e-5 &&
std::abs(krLET[oilPhaseIdx] - krog_let) <= 1e-5,
"Inconsistent LET relative permeabilities family 1 gas/oil");
BOOST_CHECK_MESSAGE(pcLET[gasPhaseIdx] - pcog_let <= 1e-2,
"Inconsistent LET capillary pressure family 1 gas/oil");
// Oil in water
Sg = 0.0;
Sw = 1.0 - So - Sg;
fs.setSaturation(waterPhaseIdx, Sw);
fs.setSaturation(oilPhaseIdx, So);
fs.setSaturation(gasPhaseIdx, Sg);
MaterialLaw::capillaryPressures(pcLET,
letmaterialLawManager.materialLawParams(elemIdx),
fs);
MaterialLaw::relativePermeabilities(krLET,
letmaterialLawManager.materialLawParams(elemIdx),
fs);
//SWOFLET
//0.1 0.2 1.5 7.0 1.5 0.5 0.05 0.15 3.5 3.0 1.3 1.0 0.7 17.0 0.95 3.8 0.04
Scalar krw_let = 0.5 * computeLetCurve((Sw - 0.2) / (1.0 - 0.2 - 0.15), 1.5, 7.0, 1.5);
Scalar krow_let = 1.0 * computeLetCurve((1.0 - Sw - 0.15) /(1.0 - 0.2 - 0.15), 3.5, 3.0, 1.3);
// S=(Sw-Swco)/(1-Swco-Sorw), Pc = Pct + (Pcir-Pct)*(1-S)^L/[(1-S)^L+E*S^T]
Scalar pcow_let = -psi2Pa*(0.04 + (3.8 - 0.04) * computeLetCurve(1.0 - (Sw - Swco) / (1.0 - 0.05 - Swco), 0.7, 17.0, 0.95));
BOOST_CHECK_MESSAGE(std::abs(krLET[waterPhaseIdx] - krw_let) <= 1e-5 &&
std::abs(krLET[oilPhaseIdx] - krow_let) <= 1e-5,
"Inconsistent LET relative permeabilities family 1 wat/oil");
BOOST_CHECK_MESSAGE(std::abs(pcLET[waterPhaseIdx] - pcow_let) <= 1e-2,
"Inconsistent LET capillary pressure family 1 wat/oil");
}
}
various improvements - add unit tests for the EclMaterialLawManager and the black-oil PVT classes - rework the API of the blackoil PVT classes so that it follows the initFromDeck() pattern and uses multiplexer classes instead of dynamic polymorphism to select the concrete PVT approach - the former change allowed to use the BlackOil fluid system with arbitrary Evaluations, i.e., it is not limited to a single Evaluation class which is specified as a class template argument anymore. - more fixes for GCC and clang warnings. now it should be reasonably hard to trigger warnings using these compilers. this commit squashes the following patches: - add a unit test for the EclMaterialLawManager - black oil: simplify the WaterPvtInterface - add a unit test for the PVT classes of the black oil fluid system - black oil: improvements to the PVT classes - blackoil PVT: remove dynamic polymorphism
2015-09-22 09:44:28 -05:00
}
Reference in New Issue Copy Permalink