Skip to content
Snippets Groups Projects
Commit 66366b9e authored by Mathis Kelm's avatar Mathis Kelm
Browse files

[pm] Fix velocity reconstruction from neumann boundary conditions

parent 7b42aa57
No related branches found
No related tags found
1 merge request!3989[pm] Fix velocity reconstruction from neumann boundary conditions
Checking pipeline status
Hide whitespace changes
Inline Side-by-side
dumux/porousmediumflow/velocity.hh
+ 1
− 1
View file @ 66366b9e
...@@ -335,7 +335,7 @@ public: ...@@ -335,7 +335,7 @@ public:
// We assume that the density at the face equals the one at the cell center and reconstruct a volume flux from the Neumann mass flux. // We assume that the density at the face equals the one at the cell center and reconstruct a volume flux from the Neumann mass flux.
const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()]; const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
const auto eqIdx = VolumeVariables::Indices::conti0EqIdx + phaseIdx; const auto eqIdx = VolumeVariables::Indices::conti0EqIdx + phaseIdx;
scvfFluxes[fIdx] += neumannFlux[eqIdx] / insideVolVars.density(phaseIdx) * scvf.area() * insideVolVars.extrusionFactor(); scvfFluxes[fIdx] += neumannFlux[eqIdx] / insideVolVars.density(phaseIdx) * scvf.area();
} }
else else
{ {
......
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