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
95 = getParamFromGroup<bool>(problem.paramGroup(),
"KOmega.EnableProductionLimiter",
false);
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;