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();