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