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

Merge pull request #3184 from bska/write-rft-seg-dyndata

Browse Source 
Finish Outputting Per Segment Dynamic State to RFT File
This commit is contained in:
Bård Skaflestad 2022-11-28 14:26:16 +01:00 committed by GitHub
commit 028ea4c281
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 430 additions and 110 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

151
src/opm/output/eclipse/WriteRFT.cpp
View File

@ -457,6 +457,126 @@ namespace {
// -----------------------------------------------------------------------
class PLTSegmentPhaseVelocity : public PLTPhaseQuantity
{
public:
explicit PLTSegmentPhaseVelocity(const std::size_t nseg = 0);
void addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol);
private:
[[nodiscard]] Opm::UnitSystem::measure oilUnit() const override
{ return Opm::UnitSystem::measure::pipeflow_velocity; }
[[nodiscard]] Opm::UnitSystem::measure gasUnit() const override
{ return Opm::UnitSystem::measure::pipeflow_velocity; }
[[nodiscard]] Opm::UnitSystem::measure waterUnit() const override
{ return Opm::UnitSystem::measure::pipeflow_velocity; }
};
PLTSegmentPhaseVelocity::PLTSegmentPhaseVelocity(const std::size_t nseg)
: PLTPhaseQuantity{nseg}
{}
void PLTSegmentPhaseVelocity::addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol)
{
using Ix = ::Opm::data::SegmentPhaseQuantity::Item;
auto velocityValue = [&segSol](const Ix i)
{
return segSol.velocity.has(i) ? -segSol.velocity.get(i) : 0.0;
};
this->addOil (usys, velocityValue(Ix::Oil));
this->addGas (usys, velocityValue(Ix::Gas));
this->addWater(usys, velocityValue(Ix::Water));
}
// -----------------------------------------------------------------------
class PLTSegmentPhaseHoldupFraction : public PLTPhaseQuantity
{
public:
explicit PLTSegmentPhaseHoldupFraction(const std::size_t nseg = 0);
void addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol);
private:
[[nodiscard]] Opm::UnitSystem::measure oilUnit() const override
{ return Opm::UnitSystem::measure::identity; }
[[nodiscard]] Opm::UnitSystem::measure gasUnit() const override
{ return Opm::UnitSystem::measure::identity; }
[[nodiscard]] Opm::UnitSystem::measure waterUnit() const override
{ return Opm::UnitSystem::measure::identity; }
};
PLTSegmentPhaseHoldupFraction::PLTSegmentPhaseHoldupFraction(const std::size_t nseg)
: PLTPhaseQuantity{nseg}
{}
void PLTSegmentPhaseHoldupFraction::addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol)
{
using Ix = ::Opm::data::SegmentPhaseQuantity::Item;
auto holdupValue = [&segSol](const Ix i)
{
return segSol.holdup.has(i) ? segSol.holdup.get(i) : 0.0;
};
this->addOil (usys, holdupValue(Ix::Oil));
this->addGas (usys, holdupValue(Ix::Gas));
this->addWater(usys, holdupValue(Ix::Water));
}
// -----------------------------------------------------------------------
class PLTSegmentPhaseViscosity : public PLTPhaseQuantity
{
public:
explicit PLTSegmentPhaseViscosity(const std::size_t nseg = 0);
void addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol);
private:
[[nodiscard]] Opm::UnitSystem::measure oilUnit() const override
{ return Opm::UnitSystem::measure::viscosity; }
[[nodiscard]] Opm::UnitSystem::measure gasUnit() const override
{ return Opm::UnitSystem::measure::viscosity; }
[[nodiscard]] Opm::UnitSystem::measure waterUnit() const override
{ return Opm::UnitSystem::measure::viscosity; }
};
PLTSegmentPhaseViscosity::PLTSegmentPhaseViscosity(const std::size_t nseg)
: PLTPhaseQuantity{nseg}
{}
void PLTSegmentPhaseViscosity::addSegment(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol)
{
using Ix = ::Opm::data::SegmentPhaseQuantity::Item;
auto viscosityValue = [&segSol](const Ix i)
{
return segSol.viscosity.has(i) ? segSol.viscosity.get(i) : 0.0;
};
this->addOil (usys, viscosityValue(Ix::Oil));
this->addGas (usys, viscosityValue(Ix::Gas));
this->addWater(usys, viscosityValue(Ix::Water));
}
// -----------------------------------------------------------------------
class PLTRecord
{
public:
@ -1030,6 +1150,9 @@ namespace {
private:
PLTFlowRate rate_{};
PLTSegmentPhaseVelocity velocity_{};
PLTSegmentPhaseHoldupFraction holdup_fraction_{};
PLTSegmentPhaseViscosity viscosity_{};
std::vector<int> neighbour_id_{};
std::vector<int> branch_id_{};
@ -1064,7 +1187,6 @@ namespace {
const ::Opm::Segment& segment);
void recordDynamicState(const ::Opm::UnitSystem& usys,
const ::Opm::Segment& segment,
const ::Opm::data::Segment& segSol);
void recordAutoICDTypeProperties(const ::Opm::UnitSystem& usys,
@ -1078,7 +1200,10 @@ namespace {
};
SegmentRecord::SegmentRecord(const std::size_t nseg)
: rate_{ nseg }
: rate_ { nseg }
, velocity_ { nseg }
, holdup_fraction_{ nseg }
, viscosity_ { nseg }
{
if (nseg == std::size_t{0}) {
return;
@ -1113,6 +1238,18 @@ namespace {
rftFile.write("SEGWRAT", this->rate_.water());
rftFile.write("SEGGRAT", this->rate_.gas());
rftFile.write("SEGOVEL", this->velocity_.oil());
rftFile.write("SEGWVEL", this->velocity_.water());
rftFile.write("SEGGVEL", this->velocity_.gas());
rftFile.write("SEGOHF", this->holdup_fraction_.oil());
rftFile.write("SEGWHF", this->holdup_fraction_.water());
rftFile.write("SEGGHF", this->holdup_fraction_.gas());
rftFile.write("SEGOVIS", this->viscosity_.oil());
rftFile.write("SEGWVIS", this->viscosity_.water());
rftFile.write("SEGGVIS", this->viscosity_.gas());
rftFile.write("SEGSSTR", this->strength_);
rftFile.write("SEGSFOPN", this->icd_setting_);
rftFile.write("SEGBRNO", this->branch_id_);
@ -1181,7 +1318,7 @@ namespace {
this->recordPhysicalLocation(usys, segments, segment);
this->recordSegmentConnectivity(segment);
this->recordSegmentProperties(usys, segment);
this->recordDynamicState(usys, segment, segSol);
this->recordDynamicState(usys, segSol);
}
void SegmentRecord::recordPhysicalLocation(const ::Opm::UnitSystem& usys,
@ -1244,15 +1381,17 @@ namespace {
}
}
void SegmentRecord::recordDynamicState(const ::Opm::UnitSystem& usys,
[[maybe_unused]] const ::Opm::Segment& segment,
const ::Opm::data::Segment& segSol)
void SegmentRecord::recordDynamicState(const ::Opm::UnitSystem& usys,
const ::Opm::data::Segment& segSol)
{
using M = ::Opm::UnitSystem::measure;
using SegPress = ::Opm::data::SegmentPressures::Value;
this->pressure_.push_back(usys.from_si(M::pressure, segSol.pressures[SegPress::Pressure]));
this->rate_.addConnection(usys, segSol.rates);
this->velocity_.addSegment(usys, segSol);
this->holdup_fraction_.addSegment(usys, segSol);
this->viscosity_.addSegment(usys, segSol);
}
void SegmentRecord::recordAutoICDTypeProperties(const ::Opm::UnitSystem& usys,

389
tests/test_RFT.cpp
View File

@ -58,6 +58,7 @@
#include <array>
#include <cstddef>
#include <ctime>
#include <functional>
#include <iomanip>
#include <iterator>
#include <map>
@ -446,6 +447,51 @@ namespace {
return this->value(segNum, this->grat_);
}
float ovel(const int segNum) const
{
return this->value(segNum, this->ovel_);
}
float wvel(const int segNum) const
{
return this->value(segNum, this->wvel_);
}
float gvel(const int segNum) const
{
return this->value(segNum, this->gvel_);
}
float hf_o(const int segNum) const
{
return this->value(segNum, this->hf_o_);
}
float hf_w(const int segNum) const
{
return this->value(segNum, this->hf_w_);
}
float hf_g(const int segNum) const
{
return this->value(segNum, this->hf_g_);
}
float ovis(const int segNum) const
{
return this->value(segNum, this->ovis_);
}
float wvis(const int segNum) const
{
return this->value(segNum, this->wvis_);
}
float gvis(const int segNum) const
{
return this->value(segNum, this->gvis_);
}
float icd_strength(const int segNum) const
{
return this->value(segNum, this->icd_strength_);
@ -487,6 +533,15 @@ namespace {
std::vector<float> orat_{};
std::vector<float> wrat_{};
std::vector<float> grat_{};
std::vector<float> ovel_{};
std::vector<float> wvel_{};
std::vector<float> gvel_{};
std::vector<float> hf_o_{};
std::vector<float> hf_w_{};
std::vector<float> hf_g_{};
std::vector<float> ovis_{};
std::vector<float> wvis_{};
std::vector<float> gvis_{};
std::vector<float> icd_strength_{};
std::vector<float> icd_setting_{};
@ -516,6 +571,15 @@ namespace {
BOOST_REQUIRE(rft.hasArray("SEGORAT", well, date));
BOOST_REQUIRE(rft.hasArray("SEGWRAT", well, date));
BOOST_REQUIRE(rft.hasArray("SEGGRAT", well, date));
BOOST_REQUIRE(rft.hasArray("SEGOVEL", well, date));
BOOST_REQUIRE(rft.hasArray("SEGWVEL", well, date));
BOOST_REQUIRE(rft.hasArray("SEGGVEL", well, date));
BOOST_REQUIRE(rft.hasArray("SEGOHF", well, date));
BOOST_REQUIRE(rft.hasArray("SEGWHF", well, date));
BOOST_REQUIRE(rft.hasArray("SEGGHF", well, date));
BOOST_REQUIRE(rft.hasArray("SEGOVIS", well, date));
BOOST_REQUIRE(rft.hasArray("SEGWVIS", well, date));
BOOST_REQUIRE(rft.hasArray("SEGGVIS", well, date));
BOOST_REQUIRE(rft.hasArray("SEGSSTR", well, date));
BOOST_REQUIRE(rft.hasArray("SEGSFOPN", well, date));
BOOST_REQUIRE(rft.hasArray("SEGBRNO", well, date));
@ -537,6 +601,18 @@ namespace {
this->wrat_ = rft.getRft<float>("SEGWRAT", well, date);
this->grat_ = rft.getRft<float>("SEGGRAT", well, date);
this->ovel_ = rft.getRft<float>("SEGOVEL", well, date);
this->wvel_ = rft.getRft<float>("SEGWVEL", well, date);
this->gvel_ = rft.getRft<float>("SEGGVEL", well, date);
this->hf_o_ = rft.getRft<float>("SEGOHF", well, date);
this->hf_w_ = rft.getRft<float>("SEGWHF", well, date);
this->hf_g_ = rft.getRft<float>("SEGGHF", well, date);
this->ovis_ = rft.getRft<float>("SEGOVIS", well, date);
this->wvis_ = rft.getRft<float>("SEGWVIS", well, date);
this->gvis_ = rft.getRft<float>("SEGGVIS", well, date);
this->icd_strength_ = rft.getRft<float>("SEGSSTR", well, date);
this->icd_setting_ = rft.getRft<float>("SEGSFOPN", well, date);
@ -4190,6 +4266,54 @@ END
return xcon;
}
Opm::data::SegmentPhaseQuantity phaseVelocity(const std::size_t segNum)
{
const auto metres_per_second = ::Opm::UnitSystem::newMETRIC()
.to_si(::Opm::UnitSystem::measure::pipeflow_velocity, 1.0);
const auto vel = - (12.0 - 1.0*(segNum - 1))*metres_per_second;
const auto v_oil = vel;
const auto v_wat = vel;
const auto v_gas = (segNum < 4) ? vel : 0.0; // No free gas in segments 4..11.
return Opm::data::SegmentPhaseQuantity{}
.set(Opm::data::SegmentPhaseQuantity::Item::Oil, v_oil)
.set(Opm::data::SegmentPhaseQuantity::Item::Gas, v_gas)
.set(Opm::data::SegmentPhaseQuantity::Item::Water, v_wat);
}
Opm::data::SegmentPhaseQuantity holdupFractions(const std::size_t segNum)
{
const auto resv_oil = 200.0 - 5*(segNum - 1);
const auto resv_wat = 100.0 - 2*(segNum - 1);
const auto resv_gas = (segNum < 4) // No free gas in segments 4..11
? 5000.0 - 100*(segNum - 1)
: 0.0;
const auto resv_tot = resv_oil + resv_wat + resv_gas;
return Opm::data::SegmentPhaseQuantity{}
.set(Opm::data::SegmentPhaseQuantity::Item::Oil, resv_oil / resv_tot)
.set(Opm::data::SegmentPhaseQuantity::Item::Gas, resv_gas / resv_tot)
.set(Opm::data::SegmentPhaseQuantity::Item::Water, resv_wat / resv_tot);
}
Opm::data::SegmentPhaseQuantity phaseViscosity(const std::size_t segNum)
{
const auto cP = ::Opm::UnitSystem::newMETRIC()
.to_si(::Opm::UnitSystem::measure::viscosity, 1.0);
const auto mu_oil = (0.25 + 0.01*(segNum - 1))*cP;
const auto mu_gas = (0.25 + 0.005*(segNum - 1))*cP;
const auto mu_wat = 0.29*cP;
return Opm::data::SegmentPhaseQuantity{}
.set(Opm::data::SegmentPhaseQuantity::Item::Oil, mu_oil)
.set(Opm::data::SegmentPhaseQuantity::Item::Gas, mu_gas)
.set(Opm::data::SegmentPhaseQuantity::Item::Water, mu_wat);
}
Opm::data::Segment segSol_P1(const std::size_t segNum)
{
const auto m3_d = ::Opm::UnitSystem::newMETRIC()
@ -4206,6 +4330,9 @@ END
xs.pressures[Opm::data::SegmentPressures::Value::Pressure] = (135.7 + 9.0*(segNum - 1))*barsa;
xs.velocity = phaseVelocity(segNum);
xs.holdup = holdupFractions(segNum);
xs.viscosity = phaseViscosity(segNum);
xs.segNumber = segNum;
return xs;
@ -4231,36 +4358,42 @@ END
return xw;
}
SegmentResults readWriteSegmentRFTData(std::function<Opm::Deck()> dataSet = &segmentDataSet,
const std::string& caseName = "TESTSEG")
{
using RftDate = ::Opm::EclIO::ERft::RftDate;
const auto rset = RSet { caseName };
const auto model = Setup{ dataSet() };
{
auto rftFile = ::Opm::EclIO::OutputStream::RFT {
rset, ::Opm::EclIO::OutputStream::Formatted { false },
::Opm::EclIO::OutputStream::RFT::OpenExisting{ false }
};
const auto reportStep = 1;
const auto elapsed = model.sched.seconds(reportStep);
const auto& grid = model.es.getInputGrid();
::Opm::RftIO::write(reportStep, elapsed, model.es.getUnits(),
grid, model.sched, wellSol(grid), rftFile);
}
const auto rft = ::Opm::EclIO::ERft {
::Opm::EclIO::OutputStream::outputFileName(rset, "RFT")
};
return SegmentResults {
rft, "P1", RftDate{ 2000, 1, 2 }
};
}
} // Anonymous namespace
BOOST_AUTO_TEST_CASE(Static_Data)
{
using RftDate = ::Opm::EclIO::ERft::RftDate;
const auto rset = RSet { "TESTSEG" };
const auto model = Setup{ segmentDataSet() };
{
auto rftFile = ::Opm::EclIO::OutputStream::RFT {
rset, ::Opm::EclIO::OutputStream::Formatted { false },
::Opm::EclIO::OutputStream::RFT::OpenExisting{ false }
};
const auto reportStep = 1;
const auto elapsed = model.sched.seconds(reportStep);
const auto& grid = model.es.getInputGrid();
::Opm::RftIO::write(reportStep, elapsed, model.es.getUnits(),
grid, model.sched, wellSol(grid), rftFile);
}
const auto rft = ::Opm::EclIO::ERft {
::Opm::EclIO::OutputStream::outputFileName(rset, "RFT")
};
const auto xSEG = SegmentResults {
rft, "P1", RftDate{ 2000, 1, 2 }
};
const auto xSEG = readWriteSegmentRFTData();
BOOST_REQUIRE_EQUAL(xSEG.numSegments(), std::size_t{11});
BOOST_REQUIRE_EQUAL(xSEG.numBranches(), std::size_t{1});
@ -4391,32 +4524,7 @@ BOOST_AUTO_TEST_CASE(Static_Data)
BOOST_AUTO_TEST_CASE(Segment_Pressure)
{
using RftDate = ::Opm::EclIO::ERft::RftDate;
const auto rset = RSet { "TESTSEG" };
const auto model = Setup{ segmentDataSet() };
{
auto rftFile = ::Opm::EclIO::OutputStream::RFT {
rset, ::Opm::EclIO::OutputStream::Formatted { false },
::Opm::EclIO::OutputStream::RFT::OpenExisting{ false }
};
const auto reportStep = 1;
const auto elapsed = model.sched.seconds(reportStep);
const auto& grid = model.es.getInputGrid();
::Opm::RftIO::write(reportStep, elapsed, model.es.getUnits(),
grid, model.sched, wellSol(grid), rftFile);
}
const auto rft = ::Opm::EclIO::ERft {
::Opm::EclIO::OutputStream::outputFileName(rset, "RFT")
};
const auto xSEG = SegmentResults {
rft, "P1", RftDate{ 2000, 1, 2 }
};
const auto xSEG = readWriteSegmentRFTData();
BOOST_CHECK_CLOSE(xSEG.pressure( 1), 135.7f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.pressure( 2), 144.7f, 1.0e-5f);
@ -4433,32 +4541,7 @@ BOOST_AUTO_TEST_CASE(Segment_Pressure)
BOOST_AUTO_TEST_CASE(Segment_Phase_Rates)
{
using RftDate = ::Opm::EclIO::ERft::RftDate;
const auto rset = RSet { "TESTSEG" };
const auto model = Setup{ segmentDataSet() };
{
auto rftFile = ::Opm::EclIO::OutputStream::RFT {
rset, ::Opm::EclIO::OutputStream::Formatted { false },
::Opm::EclIO::OutputStream::RFT::OpenExisting{ false }
};
const auto reportStep = 1;
const auto elapsed = model.sched.seconds(reportStep);
const auto& grid = model.es.getInputGrid();
::Opm::RftIO::write(reportStep, elapsed, model.es.getUnits(),
grid, model.sched, wellSol(grid), rftFile);
}
const auto rft = ::Opm::EclIO::ERft {
::Opm::EclIO::OutputStream::outputFileName(rset, "RFT")
};
const auto xSEG = SegmentResults {
rft, "P1", RftDate{ 2000, 1, 2 }
};
const auto xSEG = readWriteSegmentRFTData();
BOOST_CHECK_CLOSE(xSEG.orat( 1), 123.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.orat( 2), 113.0f, 1.0e-5f);
@ -4497,34 +4580,132 @@ BOOST_AUTO_TEST_CASE(Segment_Phase_Rates)
BOOST_CHECK_CLOSE(xSEG.wrat(11), 11.10f, 1.0e-5f);
}
BOOST_AUTO_TEST_CASE(Segment_Phase_Velocity)
{
const auto xSEG = readWriteSegmentRFTData();
BOOST_CHECK_CLOSE(xSEG.ovel( 1), 12.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 2), 11.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 3), 10.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 4), 9.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 5), 8.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 6), 7.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 7), 6.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 8), 5.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel( 9), 4.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel(10), 3.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovel(11), 2.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 1), 12.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 2), 11.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 3), 10.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 4), 9.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 5), 8.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 6), 7.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 7), 6.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 8), 5.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel( 9), 4.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel(10), 3.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvel(11), 2.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 1), 12.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 2), 11.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 3), 10.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 4), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 5), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 6), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 7), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 8), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel( 9), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel(10), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvel(11), 0.0f, 1.0e-5f);
}
BOOST_AUTO_TEST_CASE(Segment_Holdup_Fractions)
{
const auto xSEG = readWriteSegmentRFTData();
BOOST_CHECK_CLOSE(xSEG.hf_o( 1), 3.773585e-2f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 2), 3.755055e-2f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 3), 3.735745e-2f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 4), 6.630824e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 5), 6.617647e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 6), 6.603774e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 7), 6.589147e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 8), 6.573705e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o( 9), 6.557377e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o(10), 6.540084e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_o(11), 6.521739e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 1), 1.886792e-2f, 3.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 2), 1.887156e-2f, 2.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 3), 1.887534e-2f, 2.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 4), 3.369176e-1f, 2.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 5), 3.382353e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 6), 3.396226e-1f, 2.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 7), 3.410853e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 8), 3.426295e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w( 9), 3.442623e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w(10), 3.459916e-1f, 2.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_w(11), 3.478261e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 1), 9.433962e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 2), 9.435779e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 3), 9.437672e-1f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 4), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 5), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 6), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 7), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 8), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g( 9), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g(10), 0.0f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.hf_g(11), 0.0f, 1.0e-5f);
}
BOOST_AUTO_TEST_CASE(Segment_Phase_Viscosity)
{
const auto xSEG = readWriteSegmentRFTData();
BOOST_CHECK_CLOSE(xSEG.ovis( 1), 0.25f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 2), 0.26f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 3), 0.27f, 1.2e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 4), 0.28f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 5), 0.29f, 1.1e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 6), 0.30f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 7), 0.31f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 8), 0.32f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis( 9), 0.33f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis(10), 0.34f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.ovis(11), 0.35f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 1), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 2), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 3), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 4), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 5), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 6), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 7), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 8), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis( 9), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis(10), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.wvis(11), 0.29f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 1), 0.250f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 2), 0.255f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 3), 0.260f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 4), 0.265f, 1.2e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 5), 0.270f, 1.2e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 6), 0.275f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 7), 0.280f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 8), 0.285f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis( 9), 0.290f, 1.0e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis(10), 0.295f, 1.1e-5f);
BOOST_CHECK_CLOSE(xSEG.gvis(11), 0.300f, 1.0e-5f);
}
BOOST_AUTO_TEST_CASE(Valve)
{
using RftDate = ::Opm::EclIO::ERft::RftDate;
const auto rset = RSet { "TESTVALVE" };
const auto model = Setup{ valveDataSet() };
{
auto rftFile = ::Opm::EclIO::OutputStream::RFT {
rset, ::Opm::EclIO::OutputStream::Formatted { false },
::Opm::EclIO::OutputStream::RFT::OpenExisting{ false }
};
const auto reportStep = 1;
const auto elapsed = model.sched.seconds(reportStep);
const auto& grid = model.es.getInputGrid();
::Opm::RftIO::write(reportStep, elapsed, model.es.getUnits(),
grid, model.sched, wellSol(grid), rftFile);
}
const auto rft = ::Opm::EclIO::ERft {
::Opm::EclIO::OutputStream::outputFileName(rset, "RFT")
};
const auto xSEG = SegmentResults {
rft, "P1", RftDate{ 2000, 1, 2 }
};
const auto xSEG = readWriteSegmentRFTData(&valveDataSet, "TESTVALVE");
auto dfltArea = [](const float diam) -> float {
return 3.14159'26535'89793'23846'26433'83279'50288f * (diam / 2) * (diam / 2);