//no diffusion of the main component, this is a hack to use normal fick's law which computes both diffusions (main and component). We only add the part from the component here
//no diffusion of the main component, this is a hack to use normal fick's law which computes both diffusions (main and component). We only add the part from the component here
if(phaseIdx==compIdx)
continue;
//we need the upwind enthapy. Even better would be the componentEnthalpy
//we need the upwind enthapy. Even better would be the componentEnthalpy
autoenthalpy=0.0;
if(diffusiveFluxes[compIdx]>0)
enthalpy+=insideEnthalpy;
...
...
@@ -304,7 +305,7 @@ public:
}
}
/*!
* \brief Calculate the source term of the equation
* \brief Calculates the source term of the equation.
*
* \param source The source term which ought to be simulated
* \param element An element which contains part of the control volume