From 64aee1982312374f76f1d30e5c2efe7e7de56c71 Mon Sep 17 00:00:00 2001 From: Kutalmis Bercin Date: Wed, 25 Oct 2023 11:51:58 +0100 Subject: [PATCH] BUG: omegaWallFunction: add missing corner weight (fixes #3003) --- .../omegaWallFunction/omegaWallFunctionFvPatchScalarField.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C index 7339e1bd88..7c5a596642 100644 --- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C +++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C @@ -6,7 +6,7 @@ \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2011-2016, 2019 OpenFOAM Foundation - Copyright (C) 2017-2022 OpenCFD Ltd. + Copyright (C) 2017-2023 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -249,7 +249,7 @@ void Foam::omegaWallFunctionFvPatchScalarField::calculate case blenderType::MAX: { // (PH:Eq. 27) - omega0[celli] += max(omegaVis, omegaLog); + omega0[celli] += w*max(omegaVis, omegaLog); break; }