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

equil initializer: use TEMPI to specify the initial temperature

Browse Source 
this also implicitly takes care of RTEMPVD, RTEMP and all the other
ways of specifying the initial reservoir temperature.
This commit is contained in:
Andreas Lauser
2018-02-08 17:59:14 +01:00
parent 3288548c7e
commit 4241986e94
1 changed files with 13 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
ebos/eclequilinitializer.hh
View File

@@ -88,16 +88,20 @@ public:
: simulator_(simulator)
{
const auto& vanguard = simulator.vanguard();
const auto& eclState = vanguard.eclState();
unsigned numElems = vanguard.grid().size(0);
unsigned numCartesianElems = vanguard.cartesianSize();
typedef typename GET_PROP_TYPE(TypeTag, FluidSystem) FluidSystem;
EQUIL::DeckDependent::InitialStateComputer<TypeTag> initialState(materialLawManager,
vanguard.eclState(),
eclState,
vanguard.grid(),
simulator.problem().gravity()[dimWorld - 1]);
const std::vector<double>& tempiData =
eclState.get3DProperties().getDoubleGridProperty("TEMPI").getData();
// copy the result into the array of initial fluid states
initialFluidStates_.resize(numCartesianElems);
for (unsigned int elemIdx = 0; elemIdx < numElems; ++elemIdx) {
@@ -127,10 +131,14 @@ public:
fluidState.setRv(0.0);
// set the temperature
// TODO Get the temperature from the initialState
Scalar T = FluidSystem::surfaceTemperature;
fluidState.setTemperature(T);
// set the temperature.
//
// TODO: setting temperature explicitly while computing static equilibirum
// for everything else is a bit inconsistent, i.e.,
// Opm::initStateEquil() should be extended to provide correct initial
// temperatures
assert(std::isfinite(tempiData[cartesianElemIdx]));
fluidState.setTemperature(tempiData[cartesianElemIdx]);
// set the phase pressures.
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)