56 GetPropType<TypeTag, Properties::GridView>,
57 GetPropType<TypeTag, Properties::ElementMapper>,
58 GetPropType<TypeTag, Properties::Scalar>>
70 enum { numPhases = FluidSystem::numPhases };
74 :
BaseType(simulator.vanguard().cartesianIndexMapper(),
75 simulator.vanguard().gridView(),
76 simulator.model().elementMapper(),
77 simulator.vanguard().eclState())
78 , simulator_(simulator)
88 if (this->enableThresholdPressure_ && !this->thpresDefault_.empty()) {
89 this->computeDefaultThresholdPressures_();
90 this->applyExplicitThresholdPressures_();
96 void computeDefaultThresholdPressures_()
98 const auto& vanguard = simulator_.vanguard();
99 const auto& gridView = vanguard.gridView();
104 ElementContext
elemCtx(simulator_);
107 const auto& stencil =
elemCtx.stencil(0);
112 unsigned i =
face.interiorIndex();
113 unsigned j =
face.exteriorIndex();
126 const Evaluation& trans = simulator_.problem().transmissibility(
elemCtx, i, j);
127 Scalar faceArea =
face.area();
128 if (std::abs(faceArea*
getValue(trans)) < 1
e-18)
155 for (
unsigned i = 0; i < this->thpresDefault_.size(); ++i)
156 this->thpresDefault_[i] = gridView.comm().max(this->thpresDefault_[i]);
158 this->logPressures();
161 const Simulator& simulator_;
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