Skip to content
Snippets Groups Projects
Commit 1bfe0fd9 authored by Katharina Heck's avatar Katharina Heck Committed by Timo Koch
Browse files

[cleanup] fix indent and comments, use better assert in mp adapter

parent b31dbfa2
No related branches found
No related tags found
2 merge requests!1337WIP Fix/dirichlet caching v2,!1195Cleanup/interfacial areas
...@@ -31,7 +31,7 @@ namespace Dumux ...@@ -31,7 +31,7 @@ namespace Dumux
{ {
/*! /*!
* \ingroup Fluidmatrixinteractions * \ingroup Fluidmatrixinteractions
* \brief Relation for a simple effective thermal conductivity * \brief Description of a interfacial area between solid and fluid phase
*/ */
template<class Scalar> template<class Scalar>
class FluidSolidInterfacialAreaShiWang class FluidSolidInterfacialAreaShiWang
......
...@@ -38,7 +38,7 @@ namespace Dumux ...@@ -38,7 +38,7 @@ namespace Dumux
template <class MaterialLaw, int numPhases> template <class MaterialLaw, int numPhases>
class MPAdapter 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> template <class MaterialLaw>
......
...@@ -149,7 +149,7 @@ public: ...@@ -149,7 +149,7 @@ public:
//heat conduction for the solid phases //heat conduction for the solid phases
for(int sPhaseIdx=0; sPhaseIdx<numEnergyEqSolid; ++sPhaseIdx) 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: ...@@ -370,16 +370,16 @@ public:
//we only need to do this for when there is more than 1 fluid phase //we only need to do this for when there is more than 1 fluid phase
if (enableChemicalNonEquilibrium) if (enableChemicalNonEquilibrium)
{ {
// Here comes the catch: We are not doing energy conservation for the whole // Here comes the catch: We are not doing energy conservation for the whole
// system, but rather for each individual phase. // system, but rather for each individual phase.
// -> Therefore the energy fluxes over each phase boundary need be // -> Therefore the energy fluxes over each phase boundary need be
// individually accounted for. // individually accounted for.
// -> Each particle crossing a phase boundary does carry some mass and // -> Each particle crossing a phase boundary does carry some mass and
// thus energy! // thus energy!
// -> Therefore, this contribution needs to be added. // -> Therefore, this contribution needs to be added.
// -> the particle always brings the energy of the originating phase. // -> 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 // -> 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 // * molMass * enthalpy of the component in the *originating* phase
const auto& fluidState = volVars.fluidState(); const auto& fluidState = volVars.fluidState();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment