66 enum { numPhases = FluidSystem::numPhases };
67 enum { oilPhaseIdx = FluidSystem::oilPhaseIdx };
68 enum { gasPhaseIdx = FluidSystem::gasPhaseIdx };
69 enum { waterPhaseIdx = FluidSystem::waterPhaseIdx };
71 enum { numComponents = FluidSystem::numComponents };
72 enum { oilCompIdx = FluidSystem::oilCompIdx };
73 enum { gasCompIdx = FluidSystem::gasCompIdx };
74 enum { waterCompIdx = FluidSystem::waterCompIdx };
76 enum { dimWorld = GridView::dimensionworld };
96 enableSaltPrecipitation,
101 template <
class EclMaterialLawManager>
103 EclMaterialLawManager& materialLawManager)
104 : simulator_(simulator)
106 const auto& vanguard = simulator.vanguard();
107 const auto& eclState = vanguard.eclState();
109 unsigned numElems = vanguard.grid().size(0);
117 CartesianIndexMapper>;
123 vanguard.cartesianMapper(),
124 simulator.problem().gravity()[dimWorld - 1],
125 simulator.problem().numPressurePointsEquil());
128 initialFluidStates_.resize(
numElems);
129 for (
unsigned int elemIdx = 0; elemIdx <
numElems; ++elemIdx) {
130 auto& fluidState = initialFluidStates_[elemIdx];
133 unsigned regionIdx = simulator_.problem().pvtRegionIndex(elemIdx);
138 if (FluidSystem::phaseIsActive(
phaseIdx))
140 else if (Indices::numPhases == 3)
141 fluidState.setSaturation(
phaseIdx, 0.0);
144 if (FluidSystem::enableDissolvedGas())
146 else if (Indices::gasEnabled && Indices::oilEnabled)
147 fluidState.setRs(0.0);
149 if (FluidSystem::enableVaporizedOil())
151 else if (Indices::gasEnabled && Indices::oilEnabled)
152 fluidState.setRv(0.0);
154 if (FluidSystem::enableVaporizedWater())
158 if (enableTemperature || enableEnergy)
159 fluidState.setTemperature(
initialState.temperature()[elemIdx]);
163 if (!FluidSystem::phaseIsActive(
phaseIdx))
167 const auto&
b = FluidSystem::inverseFormationVolumeFactor(fluidState,
phaseIdx,
regionIdx);
171 fluidState.setDensity(
phaseIdx, rho);
180 if constexpr (enableBrine)
181 fluidState.setSaltConcentration(
initialState.saltConcentration()[elemIdx]);
184 if constexpr (enableSaltPrecipitation)
185 fluidState.setSaltSaturation(
initialState.saltSaturation()[elemIdx]);
197 return initialFluidStates_[elemIdx];
203 std::vector<ScalarFluidState> initialFluidStates_;
const ScalarFluidState & initialFluidState(unsigned elemIdx) const
Return the initial thermodynamic state which should be used as the initial condition.
Definition EquilInitializer.hpp:195
constexpr auto getPropValue()
get the value data member of a property
Definition propertysystem.hh:242
typename Properties::Detail::GetPropImpl< TypeTag, Property >::type::type GetPropType
get the type alias defined in the property (equivalent to old macro GET_PROP_TYPE(....
Definition propertysystem.hh:235