diff --git a/dumux/multidomain/2cnistokes2p2cni/2p2cnicouplinglocalresidual.hh b/dumux/multidomain/2cnistokes2p2cni/2p2cnicouplinglocalresidual.hh index 4487c37390814ccd395adf82854cddf77a51a7cc..f3422f812ddda6061f52a8b3bf0fcba1ce79b612 100644 --- a/dumux/multidomain/2cnistokes2p2cni/2p2cnicouplinglocalresidual.hh +++ b/dumux/multidomain/2cnistokes2p2cni/2p2cnicouplinglocalresidual.hh @@ -51,24 +51,13 @@ class TwoPTwoCNICouplingLocalResidual : public NILocalResidual enum { numEq = GET_PROP_VALUE(TypeTag, NumEq) }; enum { numPhases = GET_PROP_VALUE(TypeTag, NumPhases) }; enum { useMoles = GET_PROP_VALUE(TypeTag, UseMoles) }; - enum { - pressureIdx = Indices::pressureIdx, - temperatureIdx = Indices::temperatureIdx - }; - enum { - nPhaseIdx = Indices::nPhaseIdx - }; + enum { nPhaseIdx = Indices::nPhaseIdx }; + enum { wCompIdx = Indices::wCompIdx, }; enum { massBalanceIdx = GET_PROP_VALUE(TypeTag, ReplaceCompEqIdx), contiWEqIdx = Indices::contiWEqIdx, energyEqIdx = Indices::energyEqIdx }; - enum { - wCompIdx = Indices::wCompIdx, - nCompIdx = Indices::nCompIdx - }; - enum { phaseIdx = nPhaseIdx }; // index of the phase for the phase flux calculation - enum { compIdx = wCompIdx}; // index of the component for the phase flux calculation typedef typename GET_PROP_TYPE(TypeTag, VolumeVariables) VolumeVariables; typedef typename GET_PROP_TYPE(TypeTag, ElementVolumeVariables) ElementVolumeVariables; @@ -123,8 +112,6 @@ public: const VolumeVariables &volVars = this->curVolVars_()[scvIdx]; // set pressure as part of the momentum coupling - static_assert(GET_PROP_VALUE(TypeTag, Formulation) == TwoPTwoCFormulation::pnsw, - "This coupling condition is only implemented for a pnsw formulation."); if (this->bcTypes_(scvIdx).isCouplingDirichlet(massBalanceIdx)) this->residual_[scvIdx][massBalanceIdx] = volVars.pressure(nPhaseIdx); diff --git a/dumux/multidomain/2cstokes2p2c/2p2ccouplinglocalresidual.hh b/dumux/multidomain/2cstokes2p2c/2p2ccouplinglocalresidual.hh index 36a614276036466a7385857893732b6d4f921f57..9d9a480da3b3b8ded2e421413809fc210067ab4b 100644 --- a/dumux/multidomain/2cstokes2p2c/2p2ccouplinglocalresidual.hh +++ b/dumux/multidomain/2cstokes2p2c/2p2ccouplinglocalresidual.hh @@ -48,26 +48,13 @@ class TwoPTwoCCouplingLocalResidual : public TwoPTwoCLocalResidual typedef typename GET_PROP_TYPE(TypeTag, Indices) Indices; enum { dim = GridView::dimension }; - enum { numEq = GET_PROP_VALUE(TypeTag, NumEq) }; - enum { numPhases = GET_PROP_VALUE(TypeTag, NumPhases) }; enum { useMoles = GET_PROP_VALUE(TypeTag, UseMoles) }; - enum { - pressureIdx = Indices::pressureIdx - }; - enum { - wPhaseIdx = Indices::wPhaseIdx, - nPhaseIdx = Indices::nPhaseIdx - }; + enum { nPhaseIdx = Indices::nPhaseIdx }; + enum { wCompIdx = Indices::wCompIdx }; enum { massBalanceIdx = GET_PROP_VALUE(TypeTag, ReplaceCompEqIdx), contiWEqIdx = Indices::contiWEqIdx }; - enum { - wCompIdx = Indices::wCompIdx, - nCompIdx = Indices::nCompIdx - }; - enum { phaseIdx = nPhaseIdx }; // index of the phase for the phase flux calculation - enum { compIdx = wCompIdx}; // index of the component for the phase flux calculation typedef typename GET_PROP_TYPE(TypeTag, VolumeVariables) VolumeVariables; typedef typename GET_PROP_TYPE(TypeTag, ElementVolumeVariables) ElementVolumeVariables; @@ -123,8 +110,6 @@ public: const VolumeVariables &volVars = this->curVolVars_()[scvIdx]; // set pressure as part of the momentum coupling - static_assert(GET_PROP_VALUE(TypeTag, Formulation) == TwoPTwoCFormulation::pnsw, - "This coupling condition is only implemented for a pnsw formulation."); if (this->bcTypes_(scvIdx).isCouplingDirichlet(massBalanceIdx)) this->residual_[scvIdx][massBalanceIdx] = volVars.pressure(nPhaseIdx);