From 958abd24839bcee1c13f510b2863c53e95b55c64 Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Thu, 2 Sep 2021 16:20:45 +0100
Subject: [PATCH] BUG: applyBoundaryLayer: disable application on disconnected
cells (fixes #1513)
---
.../preProcessing/applyBoundaryLayer/applyBoundaryLayer.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
index b7dd15df497..bd44fb76728 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
@@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation
- Copyright (C) 2015-2019 OpenCFD Ltd.
+ Copyright (C) 2015-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -333,10 +333,10 @@ int main(int argc, char *argv[])
// u/U0 = (y/ybl)^(1/7)
// assumes U0 is the same as the current cell velocity
Info<< "Setting boundary layer velocity" << nl << endl;
- scalar yblv = ybl.value();
+ const scalar yblv = ybl.value();
forAll(U, celli)
{
- if (y[celli] <= yblv)
+ if ((y[celli] > 0) && (y[celli] <= yblv))
{
mask[celli] = 1;
U[celli] *= ::pow(y[celli]/yblv, (1.0/7.0));
--
GitLab