diff --git a/applications/solvers/combustion/reactingFoam/chemistry.H b/applications/solvers/combustion/reactingFoam/chemistry.H index 8bfa872d7b395dd0907d4fdb6355322d97853731..d7faf86b0c56719b3be3d75a91b692fae6df92b6 100644 --- a/applications/solvers/combustion/reactingFoam/chemistry.H +++ b/applications/solvers/combustion/reactingFoam/chemistry.H @@ -11,7 +11,7 @@ if (turbulentReaction) { volScalarField tk = - Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon()); + Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon()); volScalarField tc = chemistry.tc(); // Chalmers PaSR model diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C index 09e07ef57b0fd353151deb5c6e788015249ca389..96d14423eca0f612986fcc706ff9a02999e9d7aa 100644 --- a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C +++ b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C @@ -268,7 +268,7 @@ kOmegaSST::kOmegaSST / max ( a1_*omega_, - F2()*sqrt(magSqr(symm(fvc::grad(U_)))) + F2()*sqrt(2.0*magSqr(symm(fvc::grad(U_)))) ) ); mut_.correctBoundaryConditions(); @@ -365,7 +365,7 @@ void kOmegaSST::correct() // Re-calculate viscosity mut_ = a1_*rho_*k_ - /max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_))))); + /max(a1_*omega_, F2()*sqrt(2.0*magSqr(symm(fvc::grad(U_))))); mut_.correctBoundaryConditions(); // Re-calculate thermal diffusivity @@ -446,7 +446,7 @@ void kOmegaSST::correct() // Re-calculate viscosity - mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(S2)); + mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(2.0*S2)); mut_.correctBoundaryConditions(); // Re-calculate thermal diffusivity diff --git a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C index 5ebacb3c2fcde4ffc4b45128d8c2499dd58239cb..6dbebeb0d76e03a888cae873eead8f3240fded84 100644 --- a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C +++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C @@ -246,7 +246,7 @@ kOmegaSST::kOmegaSST / max ( a1_*omega_, - F2()*mag(symm(fvc::grad(U_))) + F2()*sqrt(2.0)*mag(symm(fvc::grad(U_))) ) ); nut_.correctBoundaryConditions(); @@ -400,7 +400,7 @@ void kOmegaSST::correct() // Re-calculate viscosity - nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(S2)); + nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(2*S2)); nut_.correctBoundaryConditions(); }