83 void update(
const ElementContext& elemCtx,
unsigned vertexIdx,
unsigned timeIdx)
87 const auto& problem = elemCtx.problem();
88 const auto& fractureMapper = problem.fractureMapper();
89 unsigned globalVertexIdx = elemCtx.globalSpaceIndex(vertexIdx, timeIdx);
91 Opm::Valgrind::SetUndefined(fractureFluidState_);
92 Opm::Valgrind::SetUndefined(fractureVolume_);
93 Opm::Valgrind::SetUndefined(fracturePorosity_);
94 Opm::Valgrind::SetUndefined(fractureIntrinsicPermeability_);
95 Opm::Valgrind::SetUndefined(fractureRelativePermeabilities_);
98 if (!fractureMapper.isFractureVertex(globalVertexIdx)) {
106 std::min<Scalar>(1.0, this->fluidState_.saturation(wettingPhaseIdx));
107 this->fluidState_.setSaturation(wettingPhaseIdx, SwMatrix);
108 this->fluidState_.setSaturation(nonWettingPhaseIdx, 1 - SwMatrix);
113 problem.fracturePorosity(elemCtx, vertexIdx, timeIdx);
114 fractureIntrinsicPermeability_ =
115 problem.fractureIntrinsicPermeability(elemCtx, vertexIdx, timeIdx);
121 const auto& vertexPos = elemCtx.pos(vertexIdx, timeIdx);
122 for (
unsigned vertex2Idx = 0; vertex2Idx < elemCtx.numDof(0); ++ vertex2Idx) {
123 unsigned globalVertex2Idx = elemCtx.globalSpaceIndex(vertex2Idx, timeIdx);
125 if (vertexIdx == vertex2Idx ||
126 !fractureMapper.isFractureEdge(globalVertexIdx, globalVertex2Idx))
129 Scalar fractureWidth =
130 problem.fractureWidth(elemCtx, vertexIdx, vertex2Idx, timeIdx);
132 auto distVec = elemCtx.pos(vertex2Idx, timeIdx);
133 distVec -= vertexPos;
135 Scalar edgeLength = distVec.two_norm();
143 fractureVolume_ += (fractureWidth / 2) * (edgeLength / 2);
153 fractureFluidState_.assign(this->fluidState_);
157 const auto& fractureMatParams =
158 problem.fractureMaterialLawParams(elemCtx, vertexIdx, timeIdx);
162 Scalar saturations[numPhases];
163 MaterialLaw::saturations(saturations, fractureMatParams,
164 fractureFluidState_);
165 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
166 fractureFluidState_.setSaturation(phaseIdx, saturations[phaseIdx]);
171 std::max<Scalar>(0.0, fractureFluidState_.saturation(wettingPhaseIdx));
172 fractureFluidState_.setSaturation(wettingPhaseIdx, SwFracture);
173 fractureFluidState_.setSaturation(nonWettingPhaseIdx, 1 - SwFracture);
176 MaterialLaw::relativePermeabilities(fractureRelativePermeabilities_,
178 fractureFluidState_);
182 fractureFluidState_.checkDefined();