84 static NumEqVector
flux(
const Problem& problem,
85 const Element& element,
86 const FVElementGeometry& fvGeometry,
87 const ElementVolumeVariables& elemVolVars,
88 const SubControlVolumeFace &scvf)
90 NumEqVector flux(0.0);
95 if (scvf.boundary() && problem.boundaryTypes(element, scvf).isOutflow(Indices::conti0EqIdx + 1))
98 const int phaseIdx = 0;
100 const auto& insideScv = fvGeometry.scv(scvf.insideScvIdx());
101 const auto& insideVolVars = elemVolVars[scvf.insideScvIdx()];
102 const auto& outsideVolVars = elemVolVars[scvf.outsideScvIdx()];
104 const Scalar insideDistance = (insideScv.dofPosition() - scvf.ipGlobal()).two_norm();
105 const Scalar insideDensity =
massOrMolarDensity(insideVolVars, referenceSystem, phaseIdx);
107 for (
int compIdx = 0; compIdx < numComponents; ++compIdx)
109 if (compIdx == FluidSystem::getMainComponent(phaseIdx))
112 const Scalar massOrMoleFractionInside =
massOrMoleFraction(insideVolVars, referenceSystem, phaseIdx, compIdx);
113 const Scalar massOrMoleFractionOutside =
massOrMoleFraction(outsideVolVars, referenceSystem, phaseIdx, compIdx);
114 const Scalar insideDiffCoeff = getEffectiveDiffusionCoefficient_(insideVolVars, phaseIdx, compIdx) * insideVolVars.extrusionFactor();
119 const Scalar outsideDensity =
massOrMolarDensity(outsideVolVars, referenceSystem, phaseIdx);
120 const Scalar avgDensity = 0.5*(insideDensity + outsideDensity);
121 flux[compIdx] = avgDensity * insideDiffCoeff
122 * (massOrMoleFractionInside - massOrMoleFractionOutside) / insideDistance;
126 const auto& outsideScv = fvGeometry.scv(scvf.outsideScvIdx());
127 const Scalar outsideDiffCoeff = getEffectiveDiffusionCoefficient_(outsideVolVars, phaseIdx, compIdx)
128 * outsideVolVars.extrusionFactor();
129 const Scalar outsideDistance = (outsideScv.dofPosition() - scvf.ipGlobal()).two_norm();
130 const Scalar outsideDensity =
massOrMolarDensity(outsideVolVars, referenceSystem, phaseIdx);
132 const Scalar avgDensity = 0.5*(insideDensity + outsideDensity);
133 const Scalar avgDiffCoeff =
harmonicMean(insideDiffCoeff, outsideDiffCoeff, insideDistance, outsideDistance);
135 flux[compIdx] = avgDensity * avgDiffCoeff
136 * (massOrMoleFractionInside - massOrMoleFractionOutside) / (insideDistance + outsideDistance);
141 const Scalar cumulativeFlux = std::accumulate(flux.begin(), flux.end(), 0.0);
142 flux[FluidSystem::getMainComponent(0)] = -cumulativeFlux;
144 flux *= Extrusion::area(fvGeometry, scvf);