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

ebos: use the TEMPI keyword to set the temperatures of the cells

Browse Source 
this allows to use a "pseudo thermal" simulator, i.e. each cell may
exhibit a different temperature, but energy is not conserved over
time...
This commit is contained in:
Andreas Lauser 2015-02-04 11:19:05 +01:00
parent a320dd0fa7
commit d83cd2b797
1 changed files with 15 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
applications/ebos/eclproblem.hh
View File

@ -76,9 +76,6 @@ namespace Opm {
namespace Properties {
NEW_TYPE_TAG(EclBaseProblem, INHERITS_FROM(EclGridManager, EclOutputBlackOil));
// The temperature inside the reservoir
NEW_PROP_TAG(Temperature);
// Write all solutions for visualization, not just the ones for the
// report steps...
NEW_PROP_TAG(EnableWriteAllSolutions);
@ -131,9 +128,6 @@ SET_BOOL_PROP(EclBaseProblem, EnablePartialRelinearization, false);
// only write the solutions for the report steps to disk
SET_BOOL_PROP(EclBaseProblem, EnableWriteAllSolutions, false);
// set the defaults for some problem specific properties
SET_SCALAR_PROP(EclBaseProblem, Temperature, 293.15);
// The default for the end time of the simulation [s]
//
// By default, stop it after the universe will probably have stopped
@ -232,8 +226,6 @@ public:
Ewoms::EclOutputBlackOilModule<TypeTag>::registerParameters();
EWOMS_REGISTER_PARAM(TypeTag, Scalar, Temperature,
"The temperature [K] in the reservoir");
EWOMS_REGISTER_PARAM(TypeTag, bool, EnableWriteAllSolutions,
"Write all solutions to disk instead of only the ones for the "
"report steps");
@ -266,8 +258,6 @@ public:
auto& simulator = this->simulator();
temperature_ = EWOMS_GET_PARAM(TypeTag, Scalar, Temperature);
// invert the direction of the gravity vector for ECL problems
// (z coodinates represent depth, not height.)
this->gravity_[dim - 1] *= -1;
@ -556,15 +546,15 @@ public:
/*!
* \copydoc FvBaseMultiPhaseProblem::temperature
*
* The black-oil model assumes constant temperature to define its
* parameters. Although temperature is thus not really used by the
* model, it gets written to the VTK output. Who nows, maybe we
* will need it one day?
*/
template <class Context>
Scalar temperature(const Context &context, int spaceIdx, int timeIdx) const
{ return temperature_; }
{
// use the temporally constant temperature, i.e. use the initial temperature of
// the DOF
int globalDofIdx = context.globalSpaceIndex(spaceIdx, timeIdx);
return initialFluidStates_[globalDofIdx].temperature(/*phaseIdx=*/0);
}
// \}
@ -990,6 +980,9 @@ private:
deck->getKeyword("PRESSURE")->getSIDoubleData();
const std::vector<double> &rsData =
deck->getKeyword("RS")->getSIDoubleData();
// initial reservoir temperature
const std::vector<double> &tempiData =
eclState->getDoubleGridProperty("TEMPI")->getData();
// make sure that the size of the data arrays is correct
#ifndef NDEBUG
@ -1010,9 +1003,13 @@ private:
assert(cartesianDofIdx <= numCartesianCells);
//////
// set temperatures
// set temperature
//////
dofFluidState.setTemperature(temperature_);
Scalar temperature = tempiData[cartesianDofIdx];
if (!std::isfinite(temperature) || temperature <= 0)
temperature = FluidSystem::surfaceTemperature;
dofFluidState.setTemperature(temperature);
//////
// set saturations
@ -1098,8 +1095,6 @@ private:
std::vector<BlackOilFluidState> initialFluidStates_;
Scalar temperature_;
EclWellManager<TypeTag> wellManager_;
EclDeckUnits<TypeTag> deckUnits_;