From e89ace02d0ada62b3bd85cba5a74f10b3a2535c8 Mon Sep 17 00:00:00 2001 From: andy Date: Mon, 17 Dec 2012 12:35:22 +0000 Subject: [PATCH] ENH: VanDriest delta - reset cutoff after usage --- .../compressible/LES/vanDriestDelta/vanDriestDelta.C | 4 +++- .../incompressible/LES/vanDriestDelta/vanDriestDelta.C | 4 +++- 2 files changed, 6 insertions(+), 2 deletions(-) diff --git a/src/turbulenceModels/compressible/LES/vanDriestDelta/vanDriestDelta.C b/src/turbulenceModels/compressible/LES/vanDriestDelta/vanDriestDelta.C index f6d85d5464..8c9137846d 100644 --- a/src/turbulenceModels/compressible/LES/vanDriestDelta/vanDriestDelta.C +++ b/src/turbulenceModels/compressible/LES/vanDriestDelta/vanDriestDelta.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -82,8 +82,10 @@ void vanDriestDelta::calcDelta() } } + scalar cutOff = wallPointYPlus::yPlusCutOff; wallPointYPlus::yPlusCutOff = 500; wallDistData y(mesh_, ystar); + wallPointYPlus::yPlusCutOff = cutOff; delta_ = min ( diff --git a/src/turbulenceModels/incompressible/LES/vanDriestDelta/vanDriestDelta.C b/src/turbulenceModels/incompressible/LES/vanDriestDelta/vanDriestDelta.C index 22aa5cb076..35cf23f686 100644 --- a/src/turbulenceModels/incompressible/LES/vanDriestDelta/vanDriestDelta.C +++ b/src/turbulenceModels/incompressible/LES/vanDriestDelta/vanDriestDelta.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -81,8 +81,10 @@ void vanDriestDelta::calcDelta() } } + scalar cutOff = wallPointYPlus::yPlusCutOff; wallPointYPlus::yPlusCutOff = 500; wallDistData y(mesh_, ystar); + wallPointYPlus::yPlusCutOff = cutOff; delta_ = min (