Skip to content
Snippets Groups Projects
Commit cb60ba16 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[freeflow] Clean-up diffusion laws

parent 1d42e7f1
No related branches found
No related tags found
1 merge request!928Freeflow/lowrekepsilonnc
...@@ -82,10 +82,10 @@ public: ...@@ -82,10 +82,10 @@ public:
//! We don't cache anything for this law //! We don't cache anything for this law
using Cache = FluxVariablesCaching::EmptyDiffusionCache; using Cache = FluxVariablesCaching::EmptyDiffusionCache;
static CellCenterPrimaryVariables diffusiveFluxForCellCenter(const Problem& problem, static CellCenterPrimaryVariables flux(const Problem& problem,
const FVElementGeometry& fvGeometry, const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars, const ElementVolumeVariables& elemVolVars,
const SubControlVolumeFace &scvf) const SubControlVolumeFace &scvf)
{ {
CellCenterPrimaryVariables flux(0.0); CellCenterPrimaryVariables flux(0.0);
...@@ -109,7 +109,7 @@ public: ...@@ -109,7 +109,7 @@ public:
return flux; return flux;
} }
const Scalar tij = transmissibility_(problem, fvGeometry, elemVolVars, scvf, compIdx); const Scalar tij = transmissibility_(fvGeometry, elemVolVars, scvf, compIdx);
const Scalar insideMoleFraction = insideVolVars.moleFraction(compIdx); const Scalar insideMoleFraction = insideVolVars.moleFraction(compIdx);
const Scalar outsideMolarDensity = scvf.boundary() ? insideVolVars.molarDensity() : outsideVolVars.molarDensity(); const Scalar outsideMolarDensity = scvf.boundary() ? insideVolVars.molarDensity() : outsideVolVars.molarDensity();
...@@ -148,8 +148,7 @@ public: ...@@ -148,8 +148,7 @@ public:
return flux; return flux;
} }
static Scalar transmissibility_(const Problem& problem, static Scalar transmissibility_(const FVElementGeometry& fvGeometry,
const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars, const ElementVolumeVariables& elemVolVars,
const SubControlVolumeFace& scvf, const SubControlVolumeFace& scvf,
const int compIdx) const int compIdx)
......
...@@ -85,10 +85,10 @@ public: ...@@ -85,10 +85,10 @@ public:
using Cache = FluxVariablesCaching::EmptyDiffusionCache; using Cache = FluxVariablesCaching::EmptyDiffusionCache;
using CacheFiller = FluxVariablesCaching::EmptyCacheFiller; using CacheFiller = FluxVariablesCaching::EmptyCacheFiller;
static CellCenterPrimaryVariables diffusiveFluxForCellCenter(const Problem& problem, static CellCenterPrimaryVariables flux(const Problem& problem,
const FVElementGeometry& fvGeometry, const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars, const ElementVolumeVariables& elemVolVars,
const SubControlVolumeFace& scvf) const SubControlVolumeFace& scvf)
{ {
//this is to calculate the maxwellStefan diffusion in a multicomponent system. //this is to calculate the maxwellStefan diffusion in a multicomponent system.
//see: Multicomponent Mass Transfer. R. Taylor u. R. Krishna. J. Wiley & Sons, New York 1993 //see: Multicomponent Mass Transfer. R. Taylor u. R. Krishna. J. Wiley & Sons, New York 1993
......
...@@ -106,7 +106,7 @@ public: ...@@ -106,7 +106,7 @@ public:
flux[Indices::replaceCompEqIdx] = std::accumulate(flux.begin(), flux.end(), 0.0); flux[Indices::replaceCompEqIdx] = std::accumulate(flux.begin(), flux.end(), 0.0);
} }
flux += MolecularDiffusionType::diffusiveFluxForCellCenter(problem, fvGeometry, elemVolVars, scvf); flux += MolecularDiffusionType::flux(problem, fvGeometry, elemVolVars, scvf);
return flux; return flux;
} }
}; };
......
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