diff --git a/src/turbulenceModels/RAS/compressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/compressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C index a156b44da3990bf1f640a9cc50c898f637c36389..5b0ee67630fa909764d9965fb3f9d49ed2c01af3 100644 --- a/src/turbulenceModels/RAS/compressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C +++ b/src/turbulenceModels/RAS/compressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C @@ -162,10 +162,10 @@ void omegaWallFunctionFvPatchScalarField::updateCoeffs() label faceCellI = patch().faceCells()[faceI]; scalar yPlus = - Cmu25*y[faceCellI]*sqrt(k[faceCellI]) + Cmu25*y[faceI]*sqrt(k[faceCellI]) /(muw[faceI]/rhow[faceI]); - omega[faceCellI] = sqrt(k[faceCellI])/(Cmu25*kappa*y[faceCellI]); + omega[faceCellI] = sqrt(k[faceCellI])/(Cmu25*kappa*y[faceI]); if (yPlus > yPlusLam) { @@ -173,7 +173,11 @@ void omegaWallFunctionFvPatchScalarField::updateCoeffs() (mutw[faceI] + muw[faceI]) *magGradUw[faceI] *Cmu25*sqrt(k[faceCellI]) - /(kappa*y[faceCellI]); + /(kappa*y[faceI]); + } + else + { + G[faceCellI] = 0.0; } } diff --git a/src/turbulenceModels/RAS/incompressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/turbulenceModels/RAS/incompressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C index 63bdc2e7d4e80080425351d2d2d08aef8fe84dcc..2fed92c105a4f1cbe30c77501fde7d3e990161bc 100644 --- a/src/turbulenceModels/RAS/incompressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C +++ b/src/turbulenceModels/RAS/incompressible/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C @@ -158,7 +158,7 @@ void omegaWallFunctionFvPatchScalarField::updateCoeffs() { label faceCellI = patch().faceCells()[faceI]; - scalar yPlus = Cmu25*y[faceCellI]*sqrt(k[faceCellI])/nuw[faceI]; + scalar yPlus = Cmu25*y[faceI]*sqrt(k[faceCellI])/nuw[faceI]; omega[faceCellI] = sqrt(k[faceCellI])/(Cmu25*kappa*y[faceCellI]); @@ -170,6 +170,10 @@ void omegaWallFunctionFvPatchScalarField::updateCoeffs() *Cmu25*sqrt(k[faceCellI]) /(kappa*y[faceCellI]); } + else + { + G[faceCellI] = 0.0; + } } // TODO: perform averaging for cells sharing more than one boundary face