Commit 7f0c936e authored by Martin Schneider's avatar Martin Schneider
Browse files

[md][ffpm] Add pm velocity function for tpfa coupling

parent f9ddb3f3
......@@ -63,7 +63,7 @@ class StokesDarcyCouplingDataBoxBase : public StokesDarcyCouplingDataImplementat
static constexpr auto freeFlowIdx = CouplingManager::freeFlowIdx;
static constexpr auto porousMediumIdx = CouplingManager::porousMediumIdx;
using VelocityVector = typename Element<freeFlowIdx>::Geometry::GlobalCoordinate;;
using VelocityVector = typename Element<freeFlowIdx>::Geometry::GlobalCoordinate;
template<std::size_t id> using BoundaryTypes = typename ProblemTraits<Problem<id>>::BoundaryTypes;
using StokesVelocityGradients = StaggeredVelocityGradients<Scalar, GridGeometry<freeFlowIdx>, BoundaryTypes<freeFlowIdx>, Indices<freeFlowIdx>>;
......
......@@ -57,6 +57,8 @@ class StokesDarcyCouplingDataTpfaBase : public StokesDarcyCouplingDataImplementa
static constexpr auto freeFlowIdx = CouplingManager::freeFlowIdx;
static constexpr auto porousMediumIdx = CouplingManager::porousMediumIdx;
using VelocityVector = typename Element<freeFlowIdx>::Geometry::GlobalCoordinate;
using AdvectionType = GetPropType<SubDomainTypeTag<porousMediumIdx>, Properties::AdvectionType>;
using DarcysLaw = DarcysLawImplementation<SubDomainTypeTag<porousMediumIdx>, GridGeometry<porousMediumIdx>::discMethod>;
using ForchheimersLaw = ForchheimersLawImplementation<SubDomainTypeTag<porousMediumIdx>, GridGeometry<porousMediumIdx>::discMethod>;
......@@ -107,6 +109,14 @@ public:
return momentumFlux;
}
/*!
* \brief Returns the averaged velocity vector at the interface of the porous medium according to darcys law
*/
VelocityVector porousMediumVelocity(const Element<freeFlowIdx>& element, const SubControlVolumeFace<freeFlowIdx>& scvf) const
{
DUNE_THROW(Dune::NotImplemented, "The calculation of tangential interface velocities is not implemented for Tpfa");
}
protected:
/*!
* \brief Returns the pressure at the interface
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment