From 6e829fb32ed451cabd01b3ad15ba19efbd60530c Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 17 Feb 2016 15:16:57 -0800
Subject: [PATCH] ENH: renaming epsilon field in applyBoundaryLayer
---
.../preProcessing/applyBoundaryLayer/applyBoundaryLayer.C | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
index 4473b8391de..3ec45ddfc91 100644
--- a/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
+++ b/applications/utilities/preProcessing/applyBoundaryLayer/applyBoundaryLayer.C
@@ -104,7 +104,7 @@ Foam::tmp<Foam::volScalarField> calcK
// Do not correct BC
// - operation may use inconsistent fields wrt these local manipulations
- // k.correctBoundaryConditions();
+ //k.correctBoundaryConditions();
correctProcessorPatches(k);
Info<< "Writing k\n" << endl;
@@ -132,6 +132,7 @@ Foam::tmp<Foam::volScalarField> calcEpsilon
scalar ce0 = ::pow(Cmu, 0.75)/kappa;
epsilon = (1 - mask)*epsilon + mask*ce0*k*sqrt(k)/min(y, ybl);
epsilon.max(SMALL);
+ epsilon.rename("epsilon");
// Do not correct BC
// - operation may use inconsistent fields wrt these local manipulations
@@ -297,6 +298,7 @@ void calcIncompressible
// Calculate nut - reference nut is calculated by the turbulence model
// on its construction
volScalarField& nut = tnut();
+
volScalarField S(mag(dev(symm(fvc::grad(U)))));
nut = (1 - mask)*nut + mask*sqr(kappa*min(y, ybl))*::sqrt(2)*S;
--
GitLab