diff --git a/dumux/material/fluidmatrixinteractions/1pia/fluidsolidinterfacialareashiwang.hh b/dumux/material/fluidmatrixinteractions/1pia/fluidsolidinterfacialareashiwang.hh index 6b0da034c25fe9a209c533a4585894bd56208dfe..ece673473bb93dd149988325de1c16bc8a217c02 100644 --- a/dumux/material/fluidmatrixinteractions/1pia/fluidsolidinterfacialareashiwang.hh +++ b/dumux/material/fluidmatrixinteractions/1pia/fluidsolidinterfacialareashiwang.hh @@ -31,7 +31,7 @@ namespace Dumux { /*! * \ingroup Fluidmatrixinteractions - * \brief Relation for a simple effective thermal conductivity + * \brief Description of a interfacial area between solid and fluid phase */ template<class Scalar> class FluidSolidInterfacialAreaShiWang diff --git a/dumux/material/fluidmatrixinteractions/mp/mpadapter.hh b/dumux/material/fluidmatrixinteractions/mp/mpadapter.hh index a0d51aaece39a8a4039d534f37c36d77e482f24d..698615cb5ea4386c568037b92b63873d765393a8 100644 --- a/dumux/material/fluidmatrixinteractions/mp/mpadapter.hh +++ b/dumux/material/fluidmatrixinteractions/mp/mpadapter.hh @@ -38,7 +38,7 @@ namespace Dumux template <class MaterialLaw, int numPhases> class MPAdapter { - static_assert(numPhases == 2, "only adapter for 2 phases is implemented"); + static_assert(AlwaysFalse<MaterialLaw>::value, "Adapter not implemented for the specified number of phases"); }; template <class MaterialLaw> diff --git a/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh b/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh index 2539b3e985aa5dfb26ed035d95e77516633c6683..01db8cbafb2827afdd4cacc76bbfa0a7ae824a28 100644 --- a/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh +++ b/dumux/porousmediumflow/nonequilibrium/thermal/localresidual.hh @@ -149,7 +149,7 @@ public: //heat conduction for the solid phases for(int sPhaseIdx=0; sPhaseIdx<numEnergyEqSolid; ++sPhaseIdx) { - flux[energyEqSolidIdx+sPhaseIdx] += fluxVars.heatConductionFlux(numPhases + sPhaseIdx); + flux[energyEqSolidIdx+sPhaseIdx] += fluxVars.heatConductionFlux(numPhases + sPhaseIdx); } } @@ -370,16 +370,16 @@ public: //we only need to do this for when there is more than 1 fluid phase if (enableChemicalNonEquilibrium) { - // Here comes the catch: We are not doing energy conservation for the whole - // system, but rather for each individual phase. - // -> Therefore the energy fluxes over each phase boundary need be - // individually accounted for. - // -> Each particle crossing a phase boundary does carry some mass and - // thus energy! - // -> Therefore, this contribution needs to be added. - // -> the particle always brings the energy of the originating phase. - // -> Energy advectivly transported into a phase = the moles of a component that go into a phase - // * molMass * enthalpy of the component in the *originating* phase + // Here comes the catch: We are not doing energy conservation for the whole + // system, but rather for each individual phase. + // -> Therefore the energy fluxes over each phase boundary need be + // individually accounted for. + // -> Each particle crossing a phase boundary does carry some mass and + // thus energy! + // -> Therefore, this contribution needs to be added. + // -> the particle always brings the energy of the originating phase. + // -> Energy advectivly transported into a phase = the moles of a component that go into a phase + // * molMass * enthalpy of the component in the *originating* phase const auto& fluidState = volVars.fluidState();