107 if (params_.saturationOutput_) {
110 if (params_.mobilityOutput_) {
113 if (params_.relativePermeabilityOutput_) {
117 if (params_.porosityOutput_) {
120 if (params_.intrinsicPermeabilityOutput_) {
123 if (params_.volumeFractionOutput_) {
127 if (params_.velocityOutput_) {
128 size_t nDof = this->simulator_.model().numGridDof();
129 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
130 fractureVelocity_[phaseIdx].resize(nDof);
131 for (
unsigned dofIdx = 0; dofIdx < nDof; ++dofIdx) {
132 fractureVelocity_[phaseIdx][dofIdx].resize(dimWorld);
133 fractureVelocity_[phaseIdx][dofIdx] = 0.0;
146 if (!Parameters::Get<Parameters::EnableVtkOutput>()) {
150 const auto& fractureMapper = elemCtx.simulator().vanguard().fractureMapper();
152 for (
unsigned i = 0; i < elemCtx.numPrimaryDof(0); ++i) {
153 unsigned I = elemCtx.globalSpaceIndex(i, 0);
154 if (!fractureMapper.isFractureVertex(I)) {
158 const auto& intQuants = elemCtx.intensiveQuantities(i, 0);
159 const auto& fs = intQuants.fractureFluidState();
161 if (params_.porosityOutput_) {
162 Opm::Valgrind::CheckDefined(intQuants.fracturePorosity());
163 fracturePorosity_[I] = intQuants.fracturePorosity();
165 if (params_.intrinsicPermeabilityOutput_) {
166 const auto& K = intQuants.fractureIntrinsicPermeability();
167 fractureIntrinsicPermeability_[I] = K[0][0];
170 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
171 if (params_.saturationOutput_) {
172 Opm::Valgrind::CheckDefined(fs.saturation(phaseIdx));
173 fractureSaturation_[phaseIdx][I] = fs.saturation(phaseIdx);
175 if (params_.mobilityOutput_) {
176 Opm::Valgrind::CheckDefined(intQuants.fractureMobility(phaseIdx));
177 fractureMobility_[phaseIdx][I] = intQuants.fractureMobility(phaseIdx);
179 if (params_.relativePermeabilityOutput_) {
180 Opm::Valgrind::CheckDefined(intQuants.fractureRelativePermeability(phaseIdx));
181 fractureRelativePermeability_[phaseIdx][I] =
182 intQuants.fractureRelativePermeability(phaseIdx);
184 if (params_.volumeFractionOutput_) {
185 Opm::Valgrind::CheckDefined(intQuants.fractureVolume());
186 fractureVolumeFraction_[I] += intQuants.fractureVolume();
191 if (params_.velocityOutput_) {
193 for (
unsigned scvfIdx = 0; scvfIdx < elemCtx.numInteriorFaces(0); ++ scvfIdx) {
194 const auto& extQuants = elemCtx.extensiveQuantities(scvfIdx, 0);
196 unsigned i = extQuants.interiorIndex();
197 unsigned I = elemCtx.globalSpaceIndex(i, 0);
199 unsigned j = extQuants.exteriorIndex();
200 unsigned J = elemCtx.globalSpaceIndex(j, 0);
202 if (!fractureMapper.isFractureEdge(I, J)) {
206 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
208 std::max<Scalar>(1e-16, std::abs(extQuants.fractureVolumeFlux(phaseIdx)));
209 Opm::Valgrind::CheckDefined(extQuants.extrusionFactor());
210 assert(extQuants.extrusionFactor() > 0);
211 weight *= extQuants.extrusionFactor();
213 Dune::FieldVector<Scalar, dim> v(extQuants.fractureFilterVelocity(phaseIdx));
216 for (
unsigned dimIdx = 0; dimIdx < dimWorld; ++dimIdx) {
217 fractureVelocity_[phaseIdx][I][dimIdx] += v[dimIdx];
218 fractureVelocity_[phaseIdx][J][dimIdx] += v[dimIdx];
221 fractureVelocityWeight_[phaseIdx][I] += weight;
222 fractureVelocityWeight_[phaseIdx][J] += weight;
238 if (params_.saturationOutput_) {
241 if (params_.mobilityOutput_) {
244 if (params_.relativePermeabilityOutput_) {
245 this->
commitPhaseBuffer_(baseWriter,
"fractureRelativePerm_%s", fractureRelativePermeability_);
248 if (params_.porosityOutput_) {
251 if (params_.intrinsicPermeabilityOutput_) {
252 this->
commitScalarBuffer_(baseWriter,
"fractureIntrinsicPerm", fractureIntrinsicPermeability_);
254 if (params_.volumeFractionOutput_) {
256 for (
unsigned I = 0; I < fractureVolumeFraction_.size(); ++I) {
257 fractureVolumeFraction_[I] /= this->simulator_.model().dofTotalVolume(I);
262 if (params_.velocityOutput_) {
263 size_t nDof = this->simulator_.model().numGridDof();
265 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
268 for (
unsigned dofIdx = 0; dofIdx < nDof; ++dofIdx) {
269 fractureVelocity_[phaseIdx][dofIdx] /=
270 std::max<Scalar>(1e-20, fractureVelocityWeight_[phaseIdx][dofIdx]);
274 snprintf(name, 512,
"fractureFilterVelocity_%s", FluidSystem::phaseName(phaseIdx).data());
276 DiscBaseOutputModule::attachVectorDofData_(baseWriter, fractureVelocity_[phaseIdx], name);