27#ifndef OPM_VTK_DISCRETE_FRACTURE_MODULE_HPP
28#define OPM_VTK_DISCRETE_FRACTURE_MODULE_HPP
30#include <dune/common/fvector.hh>
32#include <opm/material/common/Valgrind.hpp>
60template <
class TypeTag>
75 static const int vtkFormat = getPropValue<TypeTag, Properties::VtkOutputFormat>();
78 enum { dim = GridView::dimension };
79 enum { dimWorld = GridView::dimensionworld };
80 enum { numPhases = getPropValue<TypeTag, Properties::NumPhases>() };
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();
162 Opm::Valgrind::CheckDefined(intQuants.fracturePorosity());
163 fracturePorosity_[I] = intQuants.fracturePorosity();
166 const auto& K = intQuants.fractureIntrinsicPermeability();
167 fractureIntrinsicPermeability_[I] = K[0][0];
170 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
172 Opm::Valgrind::CheckDefined(fs.saturation(phaseIdx));
173 fractureSaturation_[phaseIdx][I] = fs.saturation(phaseIdx);
176 Opm::Valgrind::CheckDefined(intQuants.fractureMobility(phaseIdx));
177 fractureMobility_[phaseIdx][I] = intQuants.fractureMobility(phaseIdx);
180 Opm::Valgrind::CheckDefined(intQuants.fractureRelativePermeability(phaseIdx));
181 fractureRelativePermeability_[phaseIdx][I] =
182 intQuants.fractureRelativePermeability(phaseIdx);
185 Opm::Valgrind::CheckDefined(intQuants.fractureVolume());
186 fractureVolumeFraction_[I] += intQuants.fractureVolume();
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;
245 this->
commitPhaseBuffer_(baseWriter,
"fractureRelativePerm_%s", fractureRelativePermeability_);
252 this->
commitScalarBuffer_(baseWriter,
"fractureIntrinsicPerm", fractureIntrinsicPermeability_);
256 for (
unsigned I = 0; I < fractureVolumeFraction_.size(); ++I) {
257 fractureVolumeFraction_[I] /= this->
simulator_.model().dofTotalVolume(I);
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);
283 PhaseBuffer fractureSaturation_{};
284 PhaseBuffer fractureMobility_{};
285 PhaseBuffer fractureRelativePermeability_{};
287 ScalarBuffer fracturePorosity_{};
288 ScalarBuffer fractureVolumeFraction_{};
289 ScalarBuffer fractureIntrinsicPermeability_{};
291 PhaseVectorBuffer fractureVelocity_{};
292 PhaseBuffer fractureVelocityWeight_{};
294 PhaseVectorBuffer potentialGradient_{};
295 PhaseBuffer potentialWeight_{};
The base class for writer modules.
Definition: baseoutputmodule.hh:67
BaseOutputWriter::ScalarBuffer ScalarBuffer
Definition: baseoutputmodule.hh:85
void commitPhaseBuffer_(BaseOutputWriter &baseWriter, const char *pattern, PhaseBuffer &buffer, BufferType bufferType=DofBuffer)
Add a phase-specific buffer to the result file.
Definition: baseoutputmodule.hh:345
void resizePhaseBuffer_(PhaseBuffer &buffer, BufferType bufferType=DofBuffer)
Allocate the space for a buffer storing a phase-specific quantity.
Definition: baseoutputmodule.hh:201
const Simulator & simulator_
Definition: baseoutputmodule.hh:434
std::array< VectorBuffer, numPhases > PhaseVectorBuffer
Definition: baseoutputmodule.hh:94
std::array< ScalarBuffer, numPhases > PhaseBuffer
Definition: baseoutputmodule.hh:90
void commitScalarBuffer_(BaseOutputWriter &baseWriter, const char *name, ScalarBuffer &buffer, BufferType bufferType=DofBuffer)
Add a buffer containing scalar quantities to the result file.
Definition: baseoutputmodule.hh:244
void resizeScalarBuffer_(ScalarBuffer &buffer, BufferType bufferType=DofBuffer)
Allocate the space for a buffer storing a scalar quantity.
Definition: baseoutputmodule.hh:156
The base class for all output writers.
Definition: baseoutputwriter.hh:44
VTK output module for quantities which make sense for all models which deal with discrete fractures i...
Definition: vtkdiscretefracturemodule.hpp:62
void allocBuffers()
Allocate memory for the scalar fields we would like to write to the VTK file.
Definition: vtkdiscretefracturemodule.hpp:105
static void registerParameters()
Register all run-time parameters for the multi-phase VTK output module.
Definition: vtkdiscretefracturemodule.hpp:96
void commitBuffers(BaseOutputWriter &baseWriter)
Add all buffers to the VTK output writer.
Definition: vtkdiscretefracturemodule.hpp:231
void processElement(const ElementContext &elemCtx)
Modify the internal buffers according to the intensive quantities relevant for an element.
Definition: vtkdiscretefracturemodule.hpp:144
VtkDiscreteFractureModule(const Simulator &simulator)
Definition: vtkdiscretefracturemodule.hpp:87
Simplifies writing multi-file VTK datasets.
Definition: vtkmultiwriter.hh:66
Declare the properties used by the infrastructure code of the finite volume discretizations.
Definition: blackoilboundaryratevector.hh:37
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
This file provides the infrastructure to retrieve run-time parameters.
The Opm property system, traits with inheritance.
Struct holding the parameters for VtkDiscreteFractureModule.
Definition: vtkdiscretefractureparams.hpp:49
void read()
Reads the parameter values from the parameter system.
bool velocityOutput_
Definition: vtkdiscretefractureparams.hpp:62
bool relativePermeabilityOutput_
Definition: vtkdiscretefractureparams.hpp:58
bool porosityOutput_
Definition: vtkdiscretefractureparams.hpp:59
static void registerParameters()
Registers the parameters in parameter system.
bool intrinsicPermeabilityOutput_
Definition: vtkdiscretefractureparams.hpp:60
bool volumeFractionOutput_
Definition: vtkdiscretefractureparams.hpp:61
bool saturationOutput_
Definition: vtkdiscretefractureparams.hpp:56
bool mobilityOutput_
Definition: vtkdiscretefractureparams.hpp:57