Skip to content
Snippets Groups Projects
Commit 958abd24 authored by Kutalmış Berçin's avatar Kutalmış Berçin
Browse files

BUG: applyBoundaryLayer: disable application on disconnected cells (fixes #1513)

parent fcd7423f
No related branches found
No related tags found
No related merge requests found
......@@ -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));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment