69                                                           const SubControlVolume& scv,
 
   70                                                           const VolumeVariables& volVars)
 const 
   72        CellCenterPrimaryVariables storage = ParentType::computeStorageForCellCenter(problem, scv, volVars);
 
   74        storage[turbulentKineticEnergyEqIdx] = volVars.turbulentKineticEnergy()*volVars.density();
 
   75        storage[dissipationEqIdx] = volVars.dissipation()*volVars.density();
 
 
   81                                                          const Element &element,
 
   82                                                          const FVElementGeometry& fvGeometry,
 
   83                                                          const ElementVolumeVariables& elemVolVars,
 
   84                                                          const ElementFaceVariables& elemFaceVars,
 
   85                                                          const SubControlVolume &scv)
 const 
   87        CellCenterPrimaryVariables source = ParentType::computeSourceForCellCenter(problem, element, fvGeometry,
 
   88                                                                                   elemVolVars, elemFaceVars, scv);
 
   91        const auto& volVars = elemVolVars[scv];
 
   94        static const auto enableKOmegaProductionLimiter
 
   96        Scalar productionTerm = 2.0 * volVars.dynamicEddyViscosity() * volVars.stressTensorScalarProduct();
 
   97        if (enableKOmegaProductionLimiter)
 
   99            Scalar productionAlternative = 20.0 * volVars.density() * volVars.betaK() * volVars.turbulentKineticEnergy() * volVars.dissipation();
 
  100            productionTerm = min(productionTerm, productionAlternative);
 
  102        source[turbulentKineticEnergyEqIdx] += productionTerm;
 
  103        source[dissipationEqIdx] += volVars.alpha() * volVars.dissipation() / volVars.turbulentKineticEnergy() * productionTerm;
 
  106        source[turbulentKineticEnergyEqIdx] -= volVars.betaK() * volVars.density() * volVars.turbulentKineticEnergy() * volVars.dissipation();
 
  107        source[dissipationEqIdx] -= volVars.betaOmega() * volVars.density() * volVars.dissipation() * volVars.dissipation();
 
  110        Scalar gradientProduct = 0.0;
 
  111        for (
unsigned int i = 0; i < ModelTraits::dim(); ++i)
 
  112            gradientProduct += volVars.storedTurbulentKineticEnergyGradient()[i]
 
  113                               * volVars.storedDissipationGradient()[i];
 
  114        if (gradientProduct > 0.0)
 
  115            source[dissipationEqIdx] += 0.125 *volVars.density() / volVars.dissipation() * gradientProduct;