Skip to content
Snippets Groups Projects
Commit 86314f1a authored by Dennis Gläser's avatar Dennis Gläser
Browse files

[tracer][localres] make upwind term work for general vol vars

This is necessary for the upwind scheme in facet coupling models, as
the upwind/downwind volume variables type can be from the other domain.
parent 05db944a
No related branches found
No related tags found
1 merge request!1288Feature/md explicit assembler
...@@ -130,7 +130,7 @@ public: ...@@ -130,7 +130,7 @@ public:
for (int compIdx = 0; compIdx < numComponents; ++compIdx) for (int compIdx = 0; compIdx < numComponents; ++compIdx)
{ {
// the physical quantities for which we perform upwinding // the physical quantities for which we perform upwinding
auto upwindTerm = [compIdx](const VolumeVariables& volVars) auto upwindTerm = [compIdx](const auto& volVars)
{ return volVars.molarDensity()*volVars.moleFraction(phaseIdx, compIdx); }; { return volVars.molarDensity()*volVars.moleFraction(phaseIdx, compIdx); };
// advective fluxes // advective fluxes
...@@ -145,7 +145,7 @@ public: ...@@ -145,7 +145,7 @@ public:
for (int compIdx = 0; compIdx < numComponents; ++compIdx) for (int compIdx = 0; compIdx < numComponents; ++compIdx)
{ {
// the physical quantities for which we perform upwinding // the physical quantities for which we perform upwinding
auto upwindTerm = [compIdx](const VolumeVariables& volVars) auto upwindTerm = [compIdx](const auto& volVars)
{ return volVars.density()*volVars.massFraction(phaseIdx, compIdx); }; { return volVars.density()*volVars.massFraction(phaseIdx, compIdx); };
// advective fluxes // advective fluxes
......
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