From 65a30f7510471f3fb958144144c8f131f40ee938 Mon Sep 17 00:00:00 2001
From: Katharina Heck <katharina.heck@iws.uni-stuttgart.de>
Date: Mon, 5 Nov 2018 17:07:08 +0100
Subject: [PATCH] [nonequilibrium] change number of transport equations to
numcomponents in case of thermal nonequilibrium but chemical equilibrium.
otherwise models with more components than phases do not work
---
dumux/porousmediumflow/nonequilibrium/model.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/porousmediumflow/nonequilibrium/model.hh b/dumux/porousmediumflow/nonequilibrium/model.hh
index 1013d40312..893f81330c 100644
--- a/dumux/porousmediumflow/nonequilibrium/model.hh
+++ b/dumux/porousmediumflow/nonequilibrium/model.hh
@@ -63,7 +63,7 @@ template<class ET, bool chem, bool therm, int numEF, int numES, NusseltFormulati
struct NonEquilibriumModelTraits : public ET
{
static constexpr int numEq() { return numEnergyEqFluid()+numEnergyEqSolid()+numTransportEq()+ET::numConstraintEq(); }
- static constexpr int numTransportEq() { return chem ? ET::numPhases()*ET::numComponents() : ET::numPhases(); }
+ static constexpr int numTransportEq() { return chem ? ET::numPhases()*ET::numComponents() : ET::numComponents(); }
static constexpr int numEnergyEqFluid() { return therm ? numEF : 0; }
static constexpr int numEnergyEqSolid() { return therm ? numES : 0; }
--
GitLab