diff --git a/dumux/freeflow/navierstokes/momentum/fluxvariables.hh b/dumux/freeflow/navierstokes/momentum/fluxvariables.hh index 815dcd88e02ea84f7d9ac5288164947e70dfab97..765620d2fe55aee919111fe35cfd917e865e48b3 100644 --- a/dumux/freeflow/navierstokes/momentum/fluxvariables.hh +++ b/dumux/freeflow/navierstokes/momentum/fluxvariables.hh @@ -428,6 +428,8 @@ public: { if (const auto& scv = fvGeometry.scv(scvf.insideScvIdx()); this->elemBcTypes()[scvf.localIndex()].isDirichlet(scv.dofAxis())) return problem.dirichlet(this->element(), scvf)[scv.dofAxis()] * this->problem().density(this->element(), scv); + else + DUNE_THROW(Dune::InvalidStateException, "Neither Dirichlet nor Neumann BC set at " << scvf.ipGlobal()); } const bool selfIsUpstream = scvf.directionSign() == sign(transportingVelocity);