68 const SubControlVolume& scv,
69 const VolumeVariables& volVars)
const
71 CellCenterPrimaryVariables storage = ParentType::computeStorageForCellCenter(problem, scv, volVars);
73 storage[turbulentKineticEnergyEqIdx] = volVars.turbulentKineticEnergy() * volVars.density();
74 storage[dissipationEqIdx] = volVars.dissipationTilde() * volVars.density();
80 const Element &element,
81 const FVElementGeometry& fvGeometry,
82 const ElementVolumeVariables& elemVolVars,
83 const ElementFaceVariables& elemFaceVars,
84 const SubControlVolume &scv)
const
86 CellCenterPrimaryVariables source = ParentType::computeSourceForCellCenter(problem, element, fvGeometry,
87 elemVolVars, elemFaceVars, scv);
89 const auto& volVars = elemVolVars[scv];
92 source[turbulentKineticEnergyEqIdx] += 2.0 * volVars.dynamicEddyViscosity()
93 * volVars.stressTensorScalarProduct();
94 source[dissipationEqIdx] += volVars.cOneEpsilon() * volVars.fOne()
95 * volVars.dissipationTilde() / volVars.turbulentKineticEnergy()
96 * 2.0 * volVars.dynamicEddyViscosity()
97 * volVars.stressTensorScalarProduct();
100 source[turbulentKineticEnergyEqIdx] -= volVars.dissipationTilde() * volVars.density();
101 source[dissipationEqIdx] -= volVars.cTwoEpsilon() * volVars.fTwo() * volVars.density()
102 * volVars.dissipationTilde() * volVars.dissipationTilde()
103 / volVars.turbulentKineticEnergy();
106 source[turbulentKineticEnergyEqIdx] -= volVars.dValue() * volVars.density();
107 source[dissipationEqIdx] += volVars.eValue() * volVars.density();