OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Adapt to inj prefix and cleanup

Browse Source
This commit is contained in:
Tor Harald Sandve 2024-08-01 11:20:00 +02:00
parent 5f8916c0c3
commit f678c082de
3 changed files with 10 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -2587,7 +2587,7 @@ namespace Opm
if (this->isInjector() && cq_s[activeCompIdx] > 0.0){
// only handles single phase injection now
assert(this->well_ecl_.injectorType() != InjectorType::MULTI);
fs.setTemperature(this->well_ecl_.temperature());
fs.setTemperature(this->well_ecl_.inj_temperature());
typedef typename std::decay<decltype(fs)>::type::Scalar FsScalar;
typename FluidSystem::template ParameterCache<FsScalar> paramCache;
const unsigned pvtRegionIdx = intQuants.pvtRegionIndex();

13
opm/simulators/wells/WellState.cpp
View File

@ -208,10 +208,7 @@ void WellState<Scalar>::initSingleInjector(const Well& well,
const SummaryState& summary_state)
{
const auto& pu = this->phase_usage_;
Scalar temp = temperature_first_connection;
if (well.hasTemperature()) {
temp = well.temperature();
}
const Scalar temp = well.hasInjTemperature() ? well.inj_temperature() : temperature_first_connection;
auto& ws = this->wells_.add(well.name(), SingleWellState<Scalar>{well.name(),
well_info,
false,
@ -239,15 +236,17 @@ void WellState<Scalar>::initSingleWell(const std::vector<Scalar>& cellPressures,
const SummaryState& summary_state)
{
Scalar pressure_first_connection = -1;
Scalar temperature_first_connection = -1;
if (!well_perf_data.empty()) {
pressure_first_connection = cellPressures[well_perf_data[0].cell_index];
temperature_first_connection = cellTemperatures[well_perf_data[0].cell_index];
}
pressure_first_connection = well_info.broadcastFirstPerforationValue(pressure_first_connection);
temperature_first_connection = well_info.broadcastFirstPerforationValue(temperature_first_connection);
if (well.isInjector()) {
Scalar temperature_first_connection = -1;
if (!well_perf_data.empty()) {
temperature_first_connection = cellTemperatures[well_perf_data[0].cell_index];
}
temperature_first_connection = well_info.broadcastFirstPerforationValue(temperature_first_connection);
this->initSingleInjector(well, well_info, pressure_first_connection, temperature_first_connection,
well_perf_data, summary_state);
} else {

4
tests/test_wellstate.cpp
View File

@ -158,9 +158,11 @@ namespace {
std::vector<double>(setup.grid.c_grid()->number_of_cells,
100.0*Opm::unit::barsa);
const auto& unit_sytem = setup.es.getDeckUnitSystem();
const double temp = unit_sytem.to_si(Opm::UnitSystem::measure::temperature, 25);
const auto ctemp =
std::vector<double>(setup.grid.c_grid()->number_of_cells,
25+Opm::unit::degCelsiusOffset);
temp);
auto wells = setup.sched.getWells(timeStep);
pinfos.resize(wells.size());