68                                                           const SubControlVolume& scv,
 
   69                                                           const VolumeVariables& volVars)
 const 
   71        CellCenterPrimaryVariables storage = ParentType::computeStorageForCellCenter(problem, scv, volVars);
 
   72        storage[viscosityTildeEqIdx] = volVars.viscosityTilde() * volVars.density();
 
 
   77                                                          const Element &element,
 
   78                                                          const FVElementGeometry& fvGeometry,
 
   79                                                          const ElementVolumeVariables& elemVolVars,
 
   80                                                          const ElementFaceVariables& elemFaceVars,
 
   81                                                          const SubControlVolume &scv)
 const 
   83        CellCenterPrimaryVariables source = ParentType::computeSourceForCellCenter(problem, element, fvGeometry,
 
   84                                                                                   elemVolVars, elemFaceVars, scv);
 
   86        const auto& volVars = elemVolVars[scv];
 
   88        source[viscosityTildeEqIdx] += volVars.cb1() * (1.0 - volVars.ft2())
 
   89                                       * volVars.stressTensorScalarProductTilde()
 
   90                                       * volVars.viscosityTilde() * volVars.density();
 
   92        source[viscosityTildeEqIdx] -= (volVars.cw1() * volVars.fW()
 
   93                                        - volVars.cb1() * volVars.ft2() / problem.karmanConstant() / problem.karmanConstant())
 
   94                                       * volVars.viscosityTilde() * volVars.viscosityTilde()
 
   95                                       / volVars.wallDistance() / volVars.wallDistance() * volVars.density();;
 
   97        for (
unsigned int axisIdx = 0; axisIdx < ModelTraits::dim(); ++axisIdx)
 
   99            source[viscosityTildeEqIdx] += volVars.cb2() / volVars.sigma()
 
  100                                           * volVars.storedViscosityTildeGradient()[axisIdx]
 
  101                                           * volVars.storedViscosityTildeGradient()[axisIdx]