71                                                           const SubControlVolume& scv,
 
   72                                                           const VolumeVariables& volVars)
 const 
   74        CellCenterPrimaryVariables storage = ParentType::computeStorageForCellCenter(problem, scv, volVars);
 
   76        storage[turbulentKineticEnergyEqIdx] = volVars.turbulentKineticEnergy()*volVars.density();
 
   77        storage[dissipationEqIdx] = volVars.dissipation()*volVars.density();
 
 
   83                                                          const Element &element,
 
   84                                                          const FVElementGeometry& fvGeometry,
 
   85                                                          const ElementVolumeVariables& elemVolVars,
 
   86                                                          const ElementFaceVariables& elemFaceVars,
 
   87                                                          const SubControlVolume &scv)
 const 
   89        CellCenterPrimaryVariables source = ParentType::computeSourceForCellCenter(problem, element, fvGeometry,
 
   90                                                                                   elemVolVars, elemFaceVars, scv);
 
   93        const auto& volVars = elemVolVars[scv];
 
   98            Scalar productionTerm = 2.0 * volVars.dynamicEddyViscosity() * volVars.stressTensorScalarProduct();
 
   99            source[turbulentKineticEnergyEqIdx] += productionTerm;
 
  100            source[dissipationEqIdx] +=  volVars.gammaBSL() / volVars.kinematicEddyViscosity() * productionTerm;
 
  103            source[turbulentKineticEnergyEqIdx] -= volVars.betaStarBSL() * volVars.density() * volVars.turbulentKineticEnergy() * volVars.dissipation();
 
  104            source[dissipationEqIdx] -= volVars.betaBSL() * volVars.density() * volVars.dissipation() * volVars.dissipation();
 
  107            Scalar gradientProduct = 0.0;
 
  108            for (
unsigned int i = 0; i < ModelTraits::dim(); ++i)
 
  109                gradientProduct += volVars.storedTurbulentKineticEnergyGradient()[i]
 
  110                                * volVars.storedDissipationGradient()[i];
 
  112            source[dissipationEqIdx] += 2.0 * volVars.density() * (1.0-volVars.F1()) * volVars.sigmaOmega2() / volVars.dissipation() * gradientProduct;
 
  117            Scalar productionTerm = 2.0 * volVars.dynamicEddyViscosity() * volVars.stressTensorScalarProduct();
 
  118            source[turbulentKineticEnergyEqIdx] += productionTerm;
 
  119            source[dissipationEqIdx] += volVars.gammaSST() / volVars.kinematicEddyViscosity() * productionTerm;
 
  122            source[turbulentKineticEnergyEqIdx] -= volVars.betaStarSST() * volVars.density() * volVars.turbulentKineticEnergy() * volVars.dissipation();
 
  123            source[dissipationEqIdx] -= volVars.betaSST() * volVars.density() * volVars.dissipation() * volVars.dissipation();
 
  126            Scalar gradientProduct = 0.0;
 
  127            for (
unsigned int i = 0; i < ModelTraits::dim(); ++i)
 
  128                gradientProduct += volVars.storedTurbulentKineticEnergyGradient()[i]
 
  129                                * volVars.storedDissipationGradient()[i];
 
  130            source[dissipationEqIdx] += 2.0 * volVars.density() * (1.0-volVars.F1()) * volVars.sigmaOmega2() / volVars.dissipation() * gradientProduct;
 
  133            DUNE_THROW(Dune::NotImplemented, 
"\nThis SST Model is not implemented.\n");