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

[mpfa][fourierslaw] pass elemVolVars to interiorBoundaryContribution() method

parent 3256fafd
No related branches found
No related tags found
2 merge requests!617[WIP] Next,!380mixeddimension-facetcoupling
...@@ -166,13 +166,14 @@ public: ...@@ -166,13 +166,14 @@ public:
return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux(); return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux();
// Handle interior boundaries // Handle interior boundaries
flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, fluxVarsCache); flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, elemVolVars, fluxVarsCache);
// return overall resulting flux // return overall resulting flux
return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux(); return useTpfaBoundary ? flux : flux + fluxVarsCache.heatNeumannFlux();
} }
static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry, static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars,
const FluxVariablesCache& fluxVarsCache) const FluxVariablesCache& fluxVarsCache)
{ {
// obtain the transmissibilites associated with all pressures // obtain the transmissibilites associated with all pressures
......
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