Skip to content
Snippets Groups Projects
Commit 0b92e472 authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[mpfa][fickslaw] pass elemVolVars to interior boundary method()

parent b6e6fbf9
No related branches found
No related tags found
Hide whitespace changes
Inline Side-by-side
dumux/discretization/cellcentered/mpfa/fickslaw.hh
+ 2
− 1
View file @ 0b92e472
...@@ -219,7 +219,7 @@ public: ...@@ -219,7 +219,7 @@ public:
return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx); return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
// Handle interior boundaries // Handle interior boundaries
flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, fluxVarsCache, getX, phaseIdx, compIdx); flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, elemVolVars, fluxVarsCache, getX, phaseIdx, compIdx);
// return overall resulting flux // return overall resulting flux
return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx); return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
...@@ -314,6 +314,7 @@ public: ...@@ -314,6 +314,7 @@ public:
template<typename GetXFunction> template<typename GetXFunction>
static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry, static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars,
const FluxVariablesCache& fluxVarsCache, const FluxVariablesCache& fluxVarsCache,
const GetXFunction& getX, const GetXFunction& getX,
unsigned int phaseIdx, unsigned int compIdx) unsigned int phaseIdx, unsigned int compIdx)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment