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

using the same types from the co2-compositional

Browse Source 
it compiles now. Not sure why the running is terminated due to NaN.
This commit is contained in:
Kai Bao 2021-11-11 23:23:39 +01:00 committed by Trine Mykkeltvedt
parent c4a708e500
commit c6064aa3e9
1 changed files with 8 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
tests/test_chiflash.cpp
View File

@ -46,15 +46,13 @@ void testChiFlash()
using FluidSystem = Opm::TwoPhaseThreeComponentFluidSystem<Scalar>;
constexpr auto numComponents = FluidSystem::numComponents;
typedef Dune::FieldVector<Scalar, numComponents> ComponentVector;
typedef Opm::CompositionalFluidState<Scalar, FluidSystem> FluidState;
using Evaluation = Opm::DenseAd::Evaluation<double, numComponents>;
typedef Dune::FieldVector<Evaluation, numComponents> ComponentVector;
typedef Opm::CompositionalFluidState<Evaluation, FluidSystem> FluidState;
// From co2-compositional branch, it uses
// typename FluidSystem::template ParameterCache<Evaluation> paramCache;
typedef typename FluidSystem::template ParameterCache<Scalar> ParameterCache;
using ComponentVector = Dune::FieldVector<Scalar, numComponents>;
using Evaluation = Opm::DenseAd::Evaluation<double, numComponents>;
typename FluidSystem::template ParameterCache<Evaluation> paramCache;
// typedef typename FluidSystem::template ParameterCache<Scalar> ParameterCache;
FluidState fs;
// TODO: no capillary pressure for now
@ -73,7 +71,7 @@ void testChiFlash()
fs.setTemperature(303);
ParameterCache paramCache;
// ParameterCache paramCache;
paramCache.updatePhase(fs, FluidSystem::oilPhaseIdx);
paramCache.updatePhase(fs, FluidSystem::gasPhaseIdx);
fs.setDensity(FluidSystem::oilPhaseIdx, FluidSystem::density(fs, paramCache, FluidSystem::oilPhaseIdx));
@ -86,10 +84,10 @@ void testChiFlash()
// Set initial K and L
for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
const Scalar Ktmp = fs.wilsonK_(compIdx);
const Evaluation Ktmp = fs.wilsonK_(compIdx);
fs.setKvalue(compIdx, Ktmp);
}
const Scalar Ltmp = -1.0;
const Evaluation Ltmp = -1.0;
fs.setLvalue(Ltmp);
const int spatialIdx = 0;