Skip to content
Snippets Groups Projects
Commit fc882257 authored by Simon Emmert's avatar Simon Emmert
Browse files

Merge branch 'fix/fickslaw-maincomponentindex' into 'master'

Use main component index instead of phase index in the box version of Fick's law.

See merge request !1710
parents b4a2a64c fb8b7e38
No related branches found
No related tags found
1 merge request!1710Use main component index instead of phase index in the box version of Fick's law.
...@@ -123,7 +123,7 @@ public: ...@@ -123,7 +123,7 @@ public:
// compute the diffusive flux // compute the diffusive flux
componentFlux[compIdx] = -1.0*rho*vtmv(scvf.unitOuterNormal(), D, gradX)*scvf.area(); componentFlux[compIdx] = -1.0*rho*vtmv(scvf.unitOuterNormal(), D, gradX)*scvf.area();
if (BalanceEqOpts::mainComponentIsBalanced(phaseIdx) && !FluidSystem::isTracerFluidSystem()) if (BalanceEqOpts::mainComponentIsBalanced(phaseIdx) && !FluidSystem::isTracerFluidSystem())
componentFlux[phaseIdx] -= componentFlux[compIdx]; componentFlux[FluidSystem::getMainComponent(phaseIdx)] -= componentFlux[compIdx];
} }
return componentFlux; return componentFlux;
} }
......
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